Tự học ôn thi HSG Quốc gia Vật lý - Động lực học giọt mưa và Hệ có khối lượng biến thiên

TỦ SÁCH CHUYÊN KHẢO VẬT LÝ CHUYÊN & OLYMPIAD

ĐỘNG LỰC HỌC GIỌT MƯA RƠI TRONG ĐÁM MÂY: PHƯƠNG TRÌNH KHỐI LƯỢNG BIẾN THIÊN, LỰC CẢN KHÍ ĐỘNG & GIỚI HẠN PHÂN RÃ

Chuyên đề: Cơ học hệ có khối lượng biến thiên & Cơ học chất lưu • Thang điểm: 20,0 điểm (Phát triển từ chuyên khảo Kvant của P. Korkin)

Tổng quan mô hình & Dữ kiện vật lý chung

Trong giáo trình cơ học lý thuyết cổ điển của A. Sommerfeld (Mechanics) cũng như các bài toán kinh điển trên tạp chí Kvant (A. Stasenko, 1982 & 2019), chuyển động của một giọt nước hình cầu rơi xuyên qua lớp sương mù tĩnh trong khi liên tục sáp nhập các vi giọt nhỏ trên đường đi từng dẫn tới một kết quả giải tích nổi tiếng: nếu bỏ qua lực cản không khí và kích thước ban đầu rất nhỏ, giọt nước rơi nhanh dần đều với gia tốc không đổi bằng $g/7$.

Tuy nhiên, để mô tả sát thực tế một giọt mưa thật rơi từ đám mây tích vũ xuống mặt đất, ta bắt buộc phải xét đồng thời sự cạnh tranh của ba cơ chế vật lý: trọng lực, lực cản khí động và ma sát nhớt của không khí, cùng quá trình tăng khối lượng kèm tổn hao xung lượng do sáp nhập các vi giọt sương. Hơn nữa, sự chênh lệch áp suất khí động với áp suất căng bề mặt Laplace sẽ quyết định kích thước và vận tốc tối đa trước khi giọt mưa bị biến dạng và vỡ vụn.

Các hằng số vật lý sử dụng xuyên suốt bài toán:

• Khối lượng riêng của nước lỏng: $\rho = 1000\text{ kg/m}^3$; gia tốc trọng trường: $g = 9{,}8\text{ m/s}^2$.
• Khối lượng riêng của không khí: $\rho_0 = 1{,}2\text{ kg/m}^3$; độ nhớt động lực học của không khí: $\eta = 2{,}0 \cdot 10^{-5}\text{ Pa}\cdot\text{s}$.
• Hệ số cản khí động của giọt gần hình cầu: $C = 0{,}50$; sức căng bề mặt của nước: $\sigma = 7{,}5 \cdot 10^{-2}\text{ N/m}$.

Hình 1. Mô hình giọt mưa rơi thẳng đứng và hấp thụ các vi giọt nước lơ lửng trên quãng đường $\Delta x = v\,\Delta t$.

Chặng 1: Động lực học giọt mưa đơn lẻ & Ngưỡng phân rã thủy động

5,0 điểm

1.1 (2,0 điểm). Xét một giọt mưa hình cầu bán kính $r$, khối lượng $m = \frac{4}{3}\pi r^3 \rho$ đang rơi thẳng đứng với vận tốc $v$. Trong khoảng thời gian ngắn $\Delta t$, giọt mưa sáp nhập thêm khối lượng $\Delta m$ từ các vi giọt đứng yên. Lực cản của không khí gồm hai thành phần: lực cản áp suất khí động ở mặt trước $F_1 = \frac{1}{2}C\rho_0 \pi r^2 v^2$ và lực ma sát nhớt Stokes $F_2 = 6\pi \eta r v$.
Viết phương trình biến thiên xung lượng tổng quát của giọt mưa. Từ đó, xét trường hợp giọt mưa đơn lẻ (tạm bỏ qua sự sáp nhập vi giọt, $\Delta m = 0$), hãy chứng minh rằng giọt mưa sẽ tiến tới vận tốc rơi giới hạn $v_{\text{gh}}$ xác định bởi:

$$ v_{\text{gh}}(r) = \sqrt{\left(\frac{6\eta}{C\rho_0 r}\right)^2 + \frac{8\rho g r}{3C\rho_0}} - \frac{6\eta}{C\rho_0 r} \tag{1} $$

1.2 (1,5 điểm). Đồ thị Hình 2 dưới đây biểu diễn sự phụ thuộc của vận tốc giới hạn $v_{\text{gh}}$ vào bán kính $r$ theo hai mô hình: có xét ma sát nhớt (đường liền nét màu xanh) và bỏ qua ma sát nhớt $\eta \to 0$ (đường nét đứt màu hồng), đối chiếu với số liệu thực nghiệm của B. J. Mason (Physics of Clouds). Hãy tính độ lớn của lực cản khí động $F_1$ và lực cản nhớt $F_2$ đối với một giọt mưa có bán kính $r_0 = 0{,}5\text{ mm}$ đang rơi với vận tốc $v = 3{,}0\text{ m/s}$ để giải thích vì sao ở miền bán kính $r \ge 0{,}5\text{ mm}$, hai đường cong gần như sát nhau và ta có thể bỏ qua thành phần ma sát nhớt $F_2$.

Hình 2. Sự phụ thuộc của vận tốc giới hạn $v_{\text{gh}}$ vào bán kính giọt mưa $r$ theo lý thuyết và thực nghiệm.

1.3 (1,5 điểm). Theo công thức (1), vận tốc giới hạn tăng theo bán kính $r$. Tuy nhiên, trong tự nhiên không tồn tại những giọt mưa khổng lồ rơi với vận tốc hàng chục mét trên giây. Nguyên nhân là khi áp suất động của luồng khí phía trước $\Delta p_{\text{động}} = \frac{1}{2}C\rho_0 v^2$ chạm ngưỡng áp suất Laplace do sức căng bề mặt $\Delta p_{\sigma} = \frac{2\sigma}{r}$, giọt nước mất tính ổn định cầu, bị dẹt ra và vỡ vụn thành nhiều giọt nhỏ.
Thiết lập hệ thức liên hệ giữa vận tốc tới hạn phân rã $v_{\text{cr}}$ và bán kính tới hạn $r_{\text{cr}}$. Áp dụng tính giá trị của vận tốc tới hạn $v_{\text{cr}}$ (theo đơn vị $\text{m/s}$) đối với một giọt mưa lớn có bán kính $r_{\text{cr}} = 5{,}0\text{ mm}$.

🎯 Kiểm tra kết quả Chặng 1 (Ý 1.3): Nhập vận tốc tới hạn phân rã $v_{\text{cr}}$ (đơn vị $\text{m/s}$) của giọt mưa có bán kính $r_{\text{cr}} = 5{,}0\text{ mm}$:

💡 Gợi ý 1: Câu hỏi định hướng tư duy (Tự suy nghĩ trước)

• Câu hỏi 1: Trong khoảng thời gian $\Delta t$, hệ gồm giọt mưa (khối lượng $m$, vận tốc $v$) và phần vi giọt sắp sáp nhập (khối lượng $\Delta m$, vận tốc ban đầu bằng $0$) chịu tác dụng của những ngoại lực nào theo phương thẳng đứng? Độ biến thiên động lượng của toàn hệ bằng bao nhiêu khi bỏ qua vô cùng bé bậc hai $\Delta m \cdot \Delta v$?
• Câu hỏi 2: Khi giọt mưa đơn lẻ ($\Delta m = 0$) đạt trạng thái rơi ổn định ($\Delta v = 0$), phương trình cân bằng lực lúc này là một phương trình bậc mấy đối với vận tốc $v$? Vì sao ta chỉ chọn nghiệm dương?
• Câu hỏi 3: Sức căng bề mặt $\sigma$ giữ cho giọt nước có dạng hình cầu nhờ áp suất phụ Laplace hướng vào tâm. Khi luồng không khí thổi ngược vào đáy giọt với áp suất động $\frac{1}{2}C\rho_0 v^2$, điều gì xảy ra tại điểm đình trệ ở đáy giọt khi hai áp suất này cân bằng nhau?

🔍 Gợi ý 2: Hướng dẫn tháo gỡ & Phương pháp giải tích

1. Phân tích lực và xung lượng: Chọn trục $Ox$ thẳng đứng hướng xuống. Trước va chạm (thời điểm $t$), tổng động lượng của giọt mưa và phần nước $\Delta m$ đứng yên là $p(t) = mv + \Delta m \cdot 0 = mv$. Sau thời gian $\Delta t$, khối lượng hợp nhất là $m + \Delta m$ chuyển động với vận tốc $v + \Delta v$, nên động lượng là:

$$ p(t+\Delta t) = (m+\Delta m)(v+\Delta v) \approx mv + m\,\Delta v + v\,\Delta m $$

Định lý biến thiên động lượng cho hệ trong thời gian $\Delta t$:

$$ m\,\Delta v = mg\,\Delta t - F_1\,\Delta t - F_2\,\Delta t - v\,\Delta m \tag{2} $$

Hình 3 (Giải tích). Cân bằng lực trên giọt mưa (trái) và sự biến dạng đáy giọt khi áp suất khí động tiệm cận áp suất Laplace (phải).

2. Giải phương trình bậc hai tìm $v_{\text{gh}}$: Khi $\Delta m = 0$ và $\Delta v = 0$, thay $m = \frac{4}{3}\pi r^3 \rho$ vào phương trình cân bằng $mg = F_1 + F_2$, rồi chia hai vế cho $\frac{1}{2}C\rho_0 \pi r^2$, ta thu được phương trình bậc hai chuẩn tắc theo $v$.

3. Ngưỡng phân rã: Từ điều kiện $\frac{1}{2}C\rho_0 v_{\text{cr}}^2 \approx \frac{2\sigma}{r_{\text{cr}}}$, rút ra $v_{\text{cr}} = \sqrt{\frac{4\sigma}{C\rho_0 r_{\text{cr}}}}$. Thay số với đơn vị SI chuẩn ($r_{\text{cr}} = 5{,}0 \cdot 10^{-3}\text{ m}$).

✅ Bài giải chi tiết & Đáp số chuẩn

Ý 1.1. Thiết lập phương trình chuyển động và vận tốc giới hạn:

Xét hệ gồm giọt mưa có khối lượng $m$ đang chuyển động với vận tốc $v$ và phần vi giọt nước có khối lượng $\Delta m$ đang đứng yên trong không khí sẽ sáp nhập vào giọt mưa trong khoảng thời gian $\Delta t$. Áp dụng định luật II Newton dạng xung lượng - động lượng theo phương thẳng đứng hướng xuống:

$$ m\,\Delta v + v\,\Delta m = \left(mg - F_1 - F_2\right)\Delta t \implies m\,\Delta v = mg\,\Delta t - v\,\Delta m - \frac{1}{2}C\rho_0 \pi r^2 v^2\,\Delta t - 6\pi\eta r v\,\Delta t $$

Khi xét một giọt mưa đơn lẻ rơi không sáp nhập ($\Delta m = 0$), giọt tăng tốc cho tới khi đạt vận tốc giới hạn $v_{\text{gh}}$ (khi đó $\Delta v = 0$). Phương trình cân bằng lực có dạng:

$$ \frac{4}{3}\pi r^3 \rho g = \frac{1}{2}C\rho_0 \pi r^2 v_{\text{gh}}^2 + 6\pi\eta r v_{\text{gh}} $$

Chia cả hai vế cho $\frac{1}{2}C\rho_0 \pi r^2$, ta đưa về phương trình bậc hai đối với $v_{\text{gh}}$:

$$ v_{\text{gh}}^2 + \frac{12\eta}{C\rho_0 r}v_{\text{gh}} - \frac{8\rho g r}{3C\rho_0} = 0 $$

Do $v_{\text{gh}} \ge 0$, chọn nghiệm dương của phương trình bậc hai, ta thu được biểu thức cần chứng minh:

$$ v_{\text{gh}}(r) = \sqrt{\left(\frac{6\eta}{C\rho_0 r}\right)^2 + \frac{8\rho g r}{3C\rho_0}} - \frac{6\eta}{C\rho_0 r} $$

Ý 1.2. Đánh giá tỷ lệ lực cản khí động và lực ma sát nhớt:

Thay các số liệu $r_0 = 0{,}5\text{ mm} = 5{,}0 \cdot 10^{-4}\text{ m}$, $v = 3{,}0\text{ m/s}$, $C = 0{,}50$, $\rho_0 = 1{,}2\text{ kg/m}^3$ và $\eta = 2{,}0 \cdot 10^{-5}\text{ Pa}\cdot\text{s}$ vào biểu thức của hai lực cản:

$$ F_1 = \frac{1}{2}C\rho_0 \pi r_0^2 v^2 = \frac{1}{2} \times 0{,}50 \times 1{,}2 \times \pi \times \left(5{,}0 \cdot 10^{-4}\right)^2 \times 3{,}0^2 \approx 2{,}12 \cdot 10^{-6}\text{ N} $$

(Lưu ý: Ngay cả khi tính theo tiết diện hiệu dụng hoặc hệ số cản thực tế ở số Reynolds trung gian khiến $F_1$ dao động trong khoảng $(2{,}1\text{–}4{,}2)\cdot 10^{-6}\text{ N}$, hãy so sánh với lực cản nhớt Stokes $F_2$):

$$ F_2 = 6\pi \eta r_0 v = 6\pi \times 2{,}0 \cdot 10^{-5} \times 5{,}0 \cdot 10^{-4} \times 3{,}0 \approx 5{,}65 \cdot 10^{-7}\text{ N} $$

Ta thấy rõ $F_1 \gg F_2$ (lực cản áp suất khí động lớn gấp khoảng $4$ đến $7{,}5$ lần lực ma sát nhớt ngay tại $r_0 = 0{,}5\text{ mm}$, và khi $r$ tăng lên tới $1\text{–}3\text{ mm}$ thì $F_1/F_2 \propto r v$ còn lớn hơn hàng chục lần). Nhờ vậy, đối với các hạt mưa kích thước điển hình ($r \ge 0{,}5\text{ mm}$), việc bỏ qua ma sát nhớt $F_2$ chỉ làm vận tốc giới hạn bị đánh giá cao hơn một lượng rất nhỏ. Khi đặt $\eta = 0$ vào công thức (1), với $r = 0{,}3\text{ mm};\, 1{,}0\text{ mm};\, 3{,}0\text{ mm}$, lý thuyết cho $v_{\text{gh}} \approx 3{,}6\text{ m/s};\, 6{,}6\text{ m/s};\, 11{,}4\text{ m/s}$, rất sát với giá trị thực nghiệm của B. J. Mason ($2{,}5\text{ m/s};\, 6{,}5\text{ m/s};\, 10\text{ m/s}$).

Ý 1.3. Điều kiện bất ổn định bề mặt và vận tốc tới hạn phân rã:

Khi giọt mưa rơi nhanh, áp suất khí động ép lên mặt dưới của giọt là $\Delta p_{\text{động}} = \frac{1}{2}C\rho_0 v^2$, trong khi áp suất căng bề mặt giữ giọt ở dạng cầu là $\Delta p_{\sigma} = \frac{2\sigma}{r}$. Giọt bắt đầu biến dạng mạnh và vỡ vụn khi hai đại lượng này ngang ngửa nhau:

$$ \frac{1}{2}C\rho_0 v_{\text{cr}}^2 \approx \frac{2\sigma}{r_{\text{cr}}} \implies v_{\text{cr}}^2 \, r_{\text{cr}} \approx \frac{4\sigma}{C\rho_0} \tag{3} $$

Thay số với $r_{\text{cr}} = 5{,}0\text{ mm} = 5{,}0 \cdot 10^{-3}\text{ m}$, $\sigma = 7{,}5 \cdot 10^{-2}\text{ N/m}$, $C = 0{,}50$, $\rho_0 = 1{,}2\text{ kg/m}^3$:

$$ v_{\text{cr}} = \sqrt{\frac{4\sigma}{C\rho_0 r_{\text{cr}}}} = \sqrt{\frac{4 \times 7{,}5 \cdot 10^{-2}}{0{,}50 \times 1{,}2 \times 5{,}0 \cdot 10^{-3}}} = \sqrt{100} = 10\text{ m/s} $$

Đáp số chuẩn Chặng 1 (Ý 1.3): Vận tốc tới hạn phân rã của giọt mưa bán kính $5{,}0\text{ mm}$ là $v_{\text{cr}} = 10\text{ m/s}$. Các giọt lớn hơn nếu hình thành sẽ bị dẹt, tăng lực cản rồi tách thành các giọt con chuyển động chậm lại theo đường cong Hình 2.

🚀 Thử thách mở rộng Chặng 1: Số Weber và Bán kính cực đại tuyệt đối của hạt mưa

Bài toán mở rộng: Trong cơ học chất lưu, tỷ số không thứ nguyên $\text{We} = \frac{\rho_0 v^2 r}{\sigma}$ được gọi là số Weber. Kết hợp điều kiện phân rã (3) với vận tốc rơi giới hạn khi bỏ qua độ nhớt $v_{\text{gh}}^2 = \frac{8\rho g r}{3C\rho_0}$, hãy tìm biểu thức giải tích xác định bán kính cực đại tuyệt đối $r_{\max}$ mà một giọt mưa rơi tự nhiên trong khí quyển có thể duy trì mà không bị vỡ. Tính giá trị số của $r_{\max}$.

Hướng dẫn: Thế $v_{\text{cr}}^2 = \frac{8\rho g r_{\max}}{3C\rho_0}$ vào hệ thức $v_{\text{cr}}^2 r_{\max} = \frac{4\sigma}{C\rho_0}$, ta triệt tiêu được cả $C$ lẫn $\rho_0$: $$ \frac{8\rho g r_{\max}^2}{3C\rho_0} = \frac{4\sigma}{C\rho_0} \implies r_{\max} = \sqrt{\frac{3\sigma}{2\rho g}} \approx 3{,}4\text{ mm}. $$ Điều này giải thích hoàn hảo vì sao hạt mưa ngoài đời thực hiếm khi có bán kính vượt quá $3{,}5\text{–}4{,}5\text{ mm}$!

Chặng 2: Quy luật phát triển bán kính & Tích phân phương trình vi phân sáp nhập

5,0 điểm

2.1 (1,5 điểm). Bây giờ ta xét quá trình giọt mưa (bán kính ban đầu $r_0$ tại $x = 0$) rơi xuyên qua đám mây chứa các vi giọt nước đứng yên có nồng độ thể tích là $\alpha$ (tức thể tích nước lỏng chứa trong một đơn vị thể tích không khí). Giả sử mọi vi giọt nằm trong hình trụ quét bởi tiết diện ngang $\pi r^2$ của giọt mưa đều sáp nhập hoàn toàn vào nó.
Hãy viết biểu thức độ tăng khối lượng $\Delta m$ trong thời gian $\Delta t$ (tương ứng quãng đường dịch chuyển $\Delta x = v\,\Delta t$), từ đó chứng minh bán kính giọt mưa tăng tuyến tính theo quãng đường rơi $x$:

$$ r(x) = r_0 + \frac{\alpha x}{4} \tag{4} $$

2.2 (1,5 điểm). Bỏ qua thành phần lực cản nhớt $F_2$, kết hợp phương trình xung lượng (2) với biểu thức $\Delta m$ và $r(x)$, hãy chứng minh phương trình vi phân cho vận tốc $v$ theo thời gian có dạng:

$$ dv = g\,dt - \frac{3}{4}\left(\alpha + \frac{C\rho_0}{2\rho}\right)\frac{v^2\,dt}{r_0 + \alpha x / 4} \tag{5} $$

Nhận xét về sự đồng dạng toán học giữa số hạng phản lực do sáp nhập vi giọt và số hạng lực cản khí động trong phương trình (5).

2.3 (2,0 điểm). Để giải phương trình phi tuyến (5), ta chuyển biến độc lập từ thời gian $t$ sang quãng đường $x$ (nhờ $dx = v\,dt$) và đặt hàm thế vận tốc $y(x) = \frac{v^2}{2g}$.
a) Chứng minh phương trình (5) quy về phương trình vi phân cấp một:

$$ \frac{dy}{dx} = 1 - \frac{a\,y}{x + b}, \quad \text{với } a = 6 + \frac{3C\rho_0}{\alpha \rho}, \quad b = \frac{4r_0}{\alpha} \tag{6} $$

b) Giải phương trình (6) bằng phép đổi biến $y(x) = z(x)\cdot(x+b)$ với điều kiện đầu $v(0) = 0$, từ đó chứng minh biểu thức tổng quát của vận tốc giọt mưa theo quãng đường rơi $x$ là:

$$ v^2(x) = \frac{8g r_0\left(1 + \frac{\alpha x}{4r_0}\right)}{7\alpha + \frac{3C\rho_0}{\rho}} \left[ 1 - \frac{1}{\left(1 + \frac{\alpha x}{4r_0}\right)^{7 + \frac{3C\rho_0}{\alpha\rho}}} \right] \tag{7} $$

c) Áp dụng công thức (7) cho bài toán cổ điển của A. Sommerfeld: giả sử bỏ qua lực cản không khí ($C \to 0$) và giọt rơi đủ xa trong lớp sương sao cho $r(x) \gg r_0$ (tức $\frac{\alpha x}{4r_0} \gg 1$), khi đó giọt chuyển động nhanh dần đều với gia tốc $a_{\text{S}} = \frac{g}{k}$. Xác định hằng số nguyên $k$.

🎯 Kiểm tra kết quả Chặng 2 (Ý 2.3c): Nhập giá trị của hệ số nguyên $k$ trong biểu thức gia tốc Sommerfeld $a_{\text{S}} = g/k$:

💡 Gợi ý 1: Câu hỏi định hướng tư duy (Tự suy nghĩ trước)

• Câu hỏi 1: Khi giọt mưa bán kính $r$ đi được đoạn $\Delta x$, thể tích vùng không gian hình trụ mà nó quét qua là bao nhiêu? Mặt khác, khi bán kính giọt tăng thêm $\Delta r$, thể tích của bản thân giọt cầu tăng thêm một lớp vỏ cầu mỏng có thể tích bằng bao nhiêu?
• Câu hỏi 2: Tại sao khi nhân hai vế của phương trình (5) với $v$, vế trái xuất hiện đại lượng $v\,dv$ còn vế phải xuất hiện $v\,dt$? Đại lượng $v\,dv$ liên hệ thế nào với vi phân $dy$ của hàm $y = \frac{v^2}{2g}$?
• Câu hỏi 3: Khi thay $y = z(x+b)$ vào phương trình $y' = 1 - \frac{ay}{x+b}$, đạo hàm $y'$ theo quy tắc tích bằng gì? Làm thế nào để tách biến $z$ và $x$ về hai vế riêng biệt?

🔍 Gợi ý 2: Hướng dẫn tháo gỡ & Phương pháp giải tích

1. Bảo toàn khối lượng khi sáp nhập: Khối lượng nước thu gom trong ống trụ quét $\Delta x = v\,\Delta t$ là $\Delta m = \alpha \rho \pi r^2 v\,\Delta t = \alpha \rho \pi r^2 \Delta x$. Mặt khác, vì giọt luôn có dạng cầu $m = \frac{4}{3}\pi r^3 \rho$, vi phân khối lượng của giọt là $\Delta m = 4\pi r^2 \rho\,\Delta r$. Cho hai biểu thức bằng nhau và giản ước $\pi r^2 \rho$, ta tìm được ngay quan hệ giữa $\Delta r$ và $\Delta x$.

2. Đổi biến sang tọa độ $x$: Nhân hai vế của (5) với $\frac{v}{g}$ và thay $v\,dt = dx$, $y = \frac{v^2}{2g} \implies dy = \frac{v\,dv}{g}$:

$$ dy = dx - \frac{3}{4}\left(\alpha + \frac{C\rho_0}{2\rho}\right)\frac{2y\,dx}{r_0 + \alpha x/4} = dx - \frac{\left(6 + \frac{3C\rho_0}{\alpha\rho}\right)y}{x + \frac{4r_0}{\alpha}}\,dx $$

3. Tích phân phương trình tách biến: Với $y = z(x+b) \implies y' = z'(x+b) + z$. Thay vào $y' = 1 - az$, ta được:

$$ (x+b)\frac{dz}{dx} = 1 - (a+1)z \implies \frac{dz}{1 - (a+1)z} = \frac{dx}{x+b} $$

✅ Bài giải chi tiết & Đáp số chuẩn

Ý 2.1. Chứng minh quy luật tăng tuyến tính của bán kính theo quãng đường:

Trong khoảng thời gian $\Delta t$, giọt mưa quét một hình trụ có đáy $\pi r^2$ và chiều cao $\Delta x = v\,\Delta t$. Thể tích nước lỏng nằm trong các vi giọt bị hấp thụ là $\Delta V_{\text{nước}} = \alpha \pi r^2 \Delta x$, suy ra độ tăng khối lượng:

$$ \Delta m = \alpha \rho \pi r^2 v\,\Delta t = \alpha \rho \pi r^2 \Delta x $$

Mặt khác, độ tăng khối lượng của giọt cầu bán kính $r$ khi bán kính tăng thêm $\Delta r$ là $\Delta m = \rho \cdot 4\pi r^2 \Delta r$. Đồng nhất hai biểu thức:

$$ 4\pi r^2 \rho\,\Delta r = \alpha \rho \pi r^2 \Delta x \implies dr = \frac{\alpha}{4}\,dx $$

Tích phân từ $x = 0$ (bán kính $r_0$) tới tọa độ $x$, ta thu được:

$$ r(x) = r_0 + \frac{\alpha x}{4} $$

Ý 2.2. Thiết lập phương trình vi phân cho vận tốc:

Bỏ qua lực cản nhớt $F_2$ trong phương trình (2), ta có:

$$ m\,dv = mg\,dt - v\,dm - \frac{1}{2}C\rho_0 \pi r^2 v^2\,dt $$

Chia cả hai vế cho $m = \frac{4}{3}\pi r^3 \rho$ và thay $dm = \alpha \rho \pi r^2 v\,dt$:

$$ dv = g\,dt - \frac{\alpha \rho \pi r^2 v^2\,dt}{\frac{4}{3}\pi r^3 \rho} - \frac{\frac{1}{2}C\rho_0 \pi r^2 v^2\,dt}{\frac{4}{3}\pi r^3 \rho} = g\,dt - \frac{3}{4}\left(\alpha + \frac{C\rho_0}{2\rho}\right)\frac{v^2\,dt}{r(x)} $$

Thay $r(x) = r_0 + \frac{\alpha x}{4}$, ta nhận được phương trình (5). Nhận xét vật lý: Cả lực hãm phản lực do sáp nhập vi giọt ($v\frac{dm}{dt} \propto r^2 v^2$) lẫn lực cản áp suất khí động ($F_1 \propto r^2 v^2$) đều tỉ lệ thuận với diện tích mặt cắt ngang $\pi r^2$ và bình phương vận tốc $v^2$, nên chúng gộp lại thành một hệ số cản hiệu dụng duy nhất!

Ý 2.3. Tích phân phương trình vi phân tổng quát:

a) Nhân hai vế của (5) với $\frac{v}{g}$ và sử dụng $v\,dt = dx$, $dy = d\left(\frac{v^2}{2g}\right) = \frac{v\,dv}{g}$:

$$ \frac{dy}{dx} = 1 - \frac{3}{4}\left(\alpha + \frac{C\rho_0}{2\rho}\right)\frac{2y}{\frac{\alpha}{4}\left(x + \frac{4r_0}{\alpha}\right)} = 1 - \frac{\left(6 + \frac{3C\rho_0}{\alpha\rho}\right)y}{x + \frac{4r_0}{\alpha}} = 1 - \frac{a\,y}{x+b} $$

trong đó $a = 6 + \frac{3C\rho_0}{\alpha\rho}$ và $b = \frac{4r_0}{\alpha}$.

b) Đặt $y = z(x+b) \implies y' = z'(x+b) + z$. Phương trình trở thành:

$$ z'(x+b) + z = 1 - az \implies \frac{dz}{1 - (a+1)z} = \frac{dx}{x+b} $$

Lấy tích phân bất định hai vế:

$$ -\frac{1}{a+1}\ln\left|1 - (a+1)z\right| = \ln(x+b) + \text{const} \implies z(x) = \frac{1}{a+1} - \frac{A}{(x+b)^{a+1}} $$

Do đó $y(x) = z(x+b) = \frac{x+b}{a+1} - \frac{A}{(x+b)^a}$. Tại $x = 0$, giọt bắt đầu rơi không vận tốc đầu nên $y(0) = 0$:

$$ 0 = \frac{b}{a+1} - \frac{A}{b^a} \implies A = \frac{b^{a+1}}{a+1} \implies y(x) = \frac{x+b}{a+1}\left[1 - \left(\frac{b}{x+b}\right)^{a+1}\right] $$

Thay $v^2 = 2gy$, $x+b = \frac{4r_0}{\alpha}\left(1 + \frac{\alpha x}{4r_0}\right)$ và $a+1 = 7 + \frac{3C\rho_0}{\alpha\rho}$ vào biểu thức trên, ta thu được nghiệm giải tích tổng quát (7):

$$ v^2(x) = \frac{8g r_0\left(1 + \frac{\alpha x}{4r_0}\right)}{7\alpha + \frac{3C\rho_0}{\rho}} \left[ 1 - \frac{1}{\left(1 + \frac{\alpha x}{4r_0}\right)^{7 + \frac{3C\rho_0}{\alpha\rho}}} \right] $$

c) Trong giới hạn cổ điển của Sommerfeld ($C = 0$ và $\frac{\alpha x}{4r_0} \gg 1$), số hạng lũy thừa nghịch đảo $\left(1 + \frac{\alpha x}{4r_0}\right)^{-7} \to 0$ và $1 + \frac{\alpha x}{4r_0} \approx \frac{\alpha x}{4r_0}$. Công thức (7) lập tức rút gọn thành:

$$ v^2(x) \approx \frac{8g r_0 \left(\frac{\alpha x}{4r_0}\right)}{7\alpha} = \frac{2gx}{7} = 2\left(\frac{g}{7}\right)x \implies a_{\text{S}} = \frac{g}{7} $$

Đáp số chuẩn Chặng 2 (Ý 2.3c): Hệ số $k = 7$ (Gia tốc rơi ổn định của giọt nước trong lớp sương tĩnh không lực cản khí động là $a_{\text{S}} = g/7$).

🚀 Thử thách mở rộng Chặng 2: Bài toán tổng quát trong không gian $n$ chiều

Bài toán mở rộng: Giả sử khối lượng của giọt tỉ lệ với lũy thừa bậc $n$ của kích thước đặc trưng ($m \propto r^n$, trong không gian 3 chiều thì $n=3$). Nếu bỏ qua lực cản khí động ($C = 0$) và xét giới hạn $r \gg r_0$, hãy chứng minh gia tốc tiệm cận của giọt theo mô hình Sommerfeld là $a = \frac{g}{2n+1}$.

Hướng dẫn: Khi $m \propto r^n \propto x^n$, phương trình xung lượng $\frac{d(mv)}{dt} = mg$ với $v = \sqrt{2ax}$ dẫn tới $\frac{d}{dx}\left(x^n \sqrt{2ax}\right)\frac{dx}{dt} = x^n g \implies \left(n + \frac{1}{2}\right)x^{n-1/2}\cdot 2a x^{1/2} = x^n g \implies (2n+1)a = g \implies a = \frac{g}{2n+1}$. Với $n = 3$, ta thu lại ngay $a = \frac{g}{7}$ mà không cần giải toàn bộ phương trình vi phân!

Chặng 3: Các giới hạn vật lý đặc biệt & Quãng đường tăng tốc đặc trưng

5,0 điểm

3.1 (2,5 điểm). Từ biểu thức tổng quát (7), hãy thực hiện phép chuyển qua giới hạn giải tích để khảo sát hai trường hợp riêng:
a) Rơi tự do trong chân không ($C \to 0$ và $\alpha \to 0$): Dùng khai triển nhị thức Newton bậc nhất $(1+u)^N \approx 1 + Nu$ khi $u \to 0$, chứng minh (7) trở về định luật rơi tự do Galileo $v^2 = 2gx$.
b) Giọt mưa đơn lẻ rơi trong không khí có lực cản nhưng không sáp nhập ($\alpha \to 0$ và $C = 0{,}50$): Dùng giới hạn cơ bản $\lim_{\alpha \to 0}\left(1 + \frac{\alpha x}{4r_0}\right)^{-\frac{3C\rho_0}{\alpha\rho}} = e^{-\frac{3C\rho_0 x}{4\rho r_0}}$, chứng minh quy luật vận tốc theo quãng đường có dạng:

$$ v^2(x) = \frac{8\rho g r_0}{3C\rho_0}\left(1 - e^{-\frac{3C\rho_0 x}{4\rho r_0}}\right) \tag{8} $$

3.2 (2,5 điểm). Từ công thức (8), vận tốc của giọt mưa đơn lẻ tiệm cận rất nhanh tới vận tốc giới hạn theo quy luật hàm mũ. Người ta định nghĩa quãng đường tăng tốc đặc trưng $L$ là độ dài mà tại đó số mũ trong hàm $e^{-x/L}$ bằng $-1$:

$$ L = \frac{4\rho r_0}{3C\rho_0} \tag{9} $$

Hãy tính giá trị của $L$ (theo đơn vị $\text{m}$, lấy đến 1 chữ số sau dấu phẩy) đối với một giọt mưa có bán kính $r_0 = 0{,}5\text{ mm}$. Nhận xét về tỷ lệ giữa quãng đường tăng tốc $L$ này so với độ cao thông thường của các đám mây (hàng kilômét).

🎯 Kiểm tra kết quả Chặng 3 (Ý 3.2): Nhập quãng đường tăng tốc đặc trưng $L$ (đơn vị $\text{m}$) của giọt mưa đơn lẻ bán kính $r_0 = 0{,}5\text{ mm}$:

💡 Gợi ý 1: Câu hỏi định hướng tư duy (Tự suy nghĩ trước)

• Câu hỏi 1: Trong công thức (7), nếu ta đặt $u = \frac{\alpha x}{4r_0}$ rất nhỏ ($u \ll 1$) và đồng thời cho $C \to 0$, biểu thức $1 - (1+u)^{-7}$ xấp xỉ bằng bao nhiêu khi giữ lại vô cùng bé bậc nhất theo $u$?
• Câu hỏi 2: Khi cho $\alpha \to 0$ nhưng giữ $C > 0$ cố định, cơ số $\left(1 + \frac{\alpha x}{4r_0}\right)$ tiến tới $1$, còn số mũ $7 + \frac{3C\rho_0}{\alpha\rho}$ tiến tới $\infty$. Đây là dạng vô định nào trong giải tích và đưa về số $e$ bằng cách nào?
• Câu hỏi 3: Khi tính $L = \frac{4\rho r_0}{3C\rho_0}$, em cần đổi $r_0 = 0{,}5\text{ mm}$ sang đơn vị $\text{m}$ như thế nào để không bị lệch bậc $10^3$?

🔍 Gợi ý 2: Hướng dẫn tháo gỡ & Phương pháp giải tích

1. Kỹ thuật lấy giới hạn $\alpha \to 0$: Ta viết lại lũy thừa trong ngoặc vuông của công thức (7) dưới dạng chuẩn của giới hạn Euler:

$$ \lim_{\alpha \to 0}\left[\left(1 + \frac{\alpha x}{4r_0}\right)^{\frac{4r_0}{\alpha x}}\right]^{-\frac{\alpha x}{4r_0}\left(7 + \frac{3C\rho_0}{\alpha\rho}\right)} = e^{-\frac{3C\rho_0 x}{4\rho r_0}} $$

Đồng thời ở thừa số phía trước của (7), khi $\alpha \to 0$:

$$ \lim_{\alpha \to 0}\frac{8g r_0\left(1 + \frac{\alpha x}{4r_0}\right)}{7\alpha + \frac{3C\rho_0}{\rho}} = \frac{8\rho g r_0}{3C\rho_0} $$

2. Thay số tính $L$: Với $\rho = 1000\text{ kg/m}^3$, $r_0 = 5{,}0 \cdot 10^{-4}\text{ m}$, $C = 0{,}50$, $\rho_0 = 1{,}2\text{ kg/m}^3$, thay trực tiếp vào (9).

✅ Bài giải chi tiết & Đáp số chuẩn

Ý 3.1a. Giới hạn rơi tự do trong chân không ($C \to 0, \alpha \to 0$):

Cho $C = 0$ trước trong công thức (7), ta có:

$$ v^2(x) = \frac{8g r_0\left(1 + \frac{\alpha x}{4r_0}\right)}{7\alpha}\left[1 - \left(1 + \frac{\alpha x}{4r_0}\right)^{-7}\right] $$

Khi $\alpha \to 0$, đặt $u = \frac{\alpha x}{4r_0} \to 0$. Áp dụng xấp xỉ Taylor $(1+u)^{-7} \approx 1 - 7u$:

$$ v^2(x) \approx \frac{8g r_0(1 + 0)}{7\alpha} \cdot \left(7 \cdot \frac{\alpha x}{4r_0}\right) = 2gx $$

Biểu thức hoàn toàn trùng khớp với công thức rơi tự do trong trường trọng lực đều khi không có lực cản môi trường.

Ý 3.1b. Giới hạn giọt mưa đơn lẻ rơi trong không khí ($\alpha \to 0, C > 0$):

Giữ $C > 0$ và cho $\alpha \to 0$ trong công thức (7):
- Thừa số thứ nhất tiến tới: $\lim_{\alpha \to 0} \frac{8g r_0(1 + \alpha x / 4r_0)}{7\alpha + 3C\rho_0/\rho} = \frac{8\rho g r_0}{3C\rho_0}$.
- Số hạng lũy thừa trong ngoặc vuông có dạng vô định $1^{\infty}$. Lấy logarit tự nhiên của $\Phi(\alpha) = \left(1 + \frac{\alpha x}{4r_0}\right)^{-\left(7 + \frac{3C\rho_0}{\alpha\rho}\right)}$:

$$ \ln \Phi(\alpha) = -\left(7 + \frac{3C\rho_0}{\alpha\rho}\right)\ln\left(1 + \frac{\alpha x}{4r_0}\right) \underset{\alpha \to 0}{\sim} -\frac{3C\rho_0}{\alpha\rho} \cdot \frac{\alpha x}{4r_0} = -\frac{3C\rho_0 x}{4\rho r_0} $$

Suy ra $\lim_{\alpha \to 0}\Phi(\alpha) = e^{-\frac{3C\rho_0 x}{4\rho r_0}}$. Kết hợp lại, ta thu được công thức (8):

$$ v^2(x) = \frac{8\rho g r_0}{3C\rho_0}\left(1 - e^{-\frac{3C\rho_0 x}{4\rho r_0}}\right) $$

Khi $x \to \infty$, $e^{-x/L} \to 0$ và $v(x) \to \sqrt{\frac{8\rho g r_0}{3C\rho_0}}$, hoàn toàn thống nhất với công thức (1) ở Chặng 1 khi bỏ qua độ nhớt ($\eta = 0$).

Ý 3.2. Tính quãng đường tăng tốc đặc trưng $L$:

Thay số với $\rho = 1000\text{ kg/m}^3$, $r_0 = 0{,}5\text{ mm} = 5{,}0 \cdot 10^{-4}\text{ m}$, $C = 0{,}50$ và $\rho_0 = 1{,}2\text{ kg/m}^3$:

$$ L = \frac{4\rho r_0}{3C\rho_0} = \frac{4 \times 1000 \times 5{,}0 \cdot 10^{-4}}{3 \times 0{,}50 \times 1{,}2} = \frac{2{,}0}{1{,}8} \approx 11{,}1\text{ m} $$

Nhận xét vật lý: Toàn bộ giai đoạn tăng tốc của một giọt mưa đơn lẻ chỉ diễn ra trong khoảng vài mét đến vài chục mét đầu tiên ($L \approx 11\text{ m}$), rất nhỏ so với độ cao hàng kilômét của đám mây!

Đáp số chuẩn Chặng 3 (Ý 3.2): Quãng đường tăng tốc đặc trưng là $L \approx 11{,}1\text{ m}$ (hoặc lấy tròn $11\text{ m}$).

🚀 Thử thách mở rộng Chặng 3: Thời gian đặc trưng để giọt mưa đạt vận tốc giới hạn

Bài toán mở rộng: Từ phương trình vi phân theo thời gian của giọt đơn lẻ $\frac{dv}{dt} = g\left(1 - \frac{v^2}{v_{\text{gh}}^2}\right)$, hãy tích phân trực tiếp để tìm hàm vận tốc theo thời gian $v(t)$ và tính thời gian $\tau_{0{,}99}$ để giọt mưa bán kính $r_0 = 0{,}5\text{ mm}$ đạt được $99\%$ vận tốc giới hạn của nó.

Hướng dẫn: Tích phân $\int_0^v \frac{dv}{1 - (v/v_{\text{gh}})^2} = \int_0^t g\,dt \implies v(t) = v_{\text{gh}}\tanh\left(\frac{gt}{v_{\text{gh}}}\right)$. Với $r_0 = 0{,}5\text{ mm}$, $v_{\text{gh}} \approx 4{,}67\text{ m/s}$. Cho $\tanh\left(\frac{g\tau}{v_{\text{gh}}}\right) = 0{,}99 \implies \frac{g\tau}{v_{\text{gh}}} \approx 2{,}65 \implies \tau_{0{,}99} \approx 1{,}26\text{ s}$. Chỉ sau hơn 1 giây rơi, giọt mưa đơn lẻ đã gần như đạt vận tốc cực đại!

Chặng 4: Vật lý đám mây thực tế, Gia tốc trong mây dày & Mật độ vi giọt

5,0 điểm

4.1 (1,5 điểm). Trong khí tượng học, lượng hơi nước ngưng tụ trong một đơn vị thể tích không khí được gọi là độ nước của mây $w$ (đơn vị $\text{g/m}^3$ hoặc $\text{kg/m}^3$), biến thiên từ vài phần trăm $\text{g/m}^3$ đến giá trị cực đại khoảng $w_{\max} = 4{,}0\text{ g/m}^3$ ở các đám mây tích vũ dày.
Viết hệ thức liên hệ giữa nồng độ thể tích $\alpha$ với $w$ và $\rho$. Từ đó, tính giá trị của $\alpha_{\max}$ ứng với $w_{\max} = 4{,}0\text{ g/m}^3$ và so sánh hai số hạng ở mẫu số của công thức (7): $7\alpha_{\max}$ và $\frac{3C\rho_0}{\rho}$. Kết luận xem hiệu ứng nào đóng vai trò chi phối lực hãm trong đám mây thực tế.

4.2 (2,0 điểm). Xét một giọt mưa ban đầu có bán kính $r_0 = 0{,}5\text{ mm}$ rơi xuyên qua một đám mây rất dày có bề dày $H = 4{,}0\text{ km}$ và nồng độ thể tích lớn $\alpha = 4{,}0 \cdot 10^{-6}$.
a) Chứng minh rằng ở phần dưới của đám mây (khi $\frac{\alpha x}{4r_0} \gg 1$), công thức (7) một lần nữa rút gọn về quy luật chuyển động nhanh dần đều với gia tốc hiệu dụng $a_{\text{hd}}$ không đổi:

$$ v^2(x) \approx \frac{2gx}{7 + \frac{3C\rho_0}{\alpha \rho}} \implies a_{\text{hd}} = \frac{g}{7 + \frac{3C\rho_0}{\alpha \rho}} \approx \frac{\alpha \rho}{3C\rho_0}\,g \tag{10} $$

Giải thích bản chất vật lý vì sao dù có lực cản không khí rất lớn, sự tăng bán kính do sáp nhập lại giúp giọt mưa duy trì gia tốc không đổi $a_{\text{hd}}$, và chứng minh rằng vận tốc theo (10) thực chất vẫn tuân theo đường cong vận tốc giới hạn tức thời $v^2(x) \approx \frac{8\rho g r(x)}{3C\rho_0}$.
b) Tính bán kính tổng cộng $r(H)$ (theo đơn vị $\text{mm}$) và vận tốc lý thuyết $v(H)$ của giọt mưa khi đi hết bề dày đám mây $H = 4{,}0\text{ km}$ (giả sử giọt chưa vỡ). Đối chiếu với tiêu chuẩn phân rã (3) ở Chặng 1, hãy cho biết giọt mưa có giữ được nguyên vẹn tới đáy mây hay không.

4.3 (1,5 điểm). Giả sử các vi giọt sương cấu tạo nên đám mây có bán kính trung bình $r_{\text{s}} \approx 10\text{ }\mu\text{m} = 1{,}0 \cdot 10^{-5}\text{ m}$ và nồng độ thể tích trong mây là $\alpha = 4{,}0 \cdot 10^{-6}$. Hãy ước lượng mật độ số lượng vi giọt $n$ (số vi giọt trong $1\text{ m}^3$ không khí) và khoảng cách trung bình $d \approx n^{-1/3}$ giữa các vi giọt lân cận.

🎯 Kiểm tra kết quả Chặng 4 (Ý 4.2b): Nhập bán kính lý thuyết $r(H)$ (đơn vị $\text{mm}$) của giọt mưa sau khi rơi qua đám mây dày $H = 4{,}0\text{ km}$ (với $r_0 = 0{,}5\text{ mm}$, $\alpha = 4{,}0 \cdot 10^{-6}$):

💡 Gợi ý 1: Câu hỏi định hướng tư duy (Tự suy nghĩ trước)

• Câu hỏi 1: Trong $1\text{ m}^3$ không khí của đám mây có chứa khối lượng nước ngưng tụ là $w$ (tính bằng $\text{kg/m}^3$). Vì khối lượng riêng của nước lỏng là $\rho = 1000\text{ kg/m}^3$, thể tích thực sự mà phần nước lỏng này chiếm chỗ là bao nhiêu?
• Câu hỏi 2: Trọng lực tác dụng lên giọt tỉ lệ với lũy thừa bậc mấy của bán kính $r$, trong khi lực cản khí động tỉ lệ với lũy thừa bậc mấy của $r$? Khi $r$ tăng tuyến tính theo quãng đường $x$, tỷ số $\frac{mg}{F_1}$ phụ thuộc vào $r$ và $v$ như thế nào?
• Câu hỏi 3: Mỗi vi giọt sương có thể tích $V_1 = \frac{4}{3}\pi r_{\text{s}}^3$. Nếu tổng thể tích của $n$ vi giọt trong $1\text{ m}^3$ không khí bằng $\alpha$, ta tính $n$ và khoảng cách trung bình $d = n^{-1/3}$ như thế nào?

🔍 Gợi ý 2: Hướng dẫn tháo gỡ & Phương pháp giải tích

1. So sánh hiệu ứng cản khí động và hiệu ứng sáp nhập: Ta có $w = \alpha \rho \implies \alpha = \frac{w}{\rho}$. Với $w_{\max} = 4{,}0\text{ g/m}^3 = 4{,}0 \cdot 10^{-3}\text{ kg/m}^3$, ta được $\alpha_{\max} = 4{,}0 \cdot 10^{-6}$. Khi đó $7\alpha_{\max} = 2{,}8 \cdot 10^{-5}$, trong khi $\frac{3C\rho_0}{\rho} = \frac{3 \times 0{,}50 \times 1{,}2}{1000} = 1{,}8 \cdot 10^{-3} \gg 7\alpha_{\max}$.

Hình 4 (Giải tích). So sánh sự tiến triển vận tốc $v(x)$ của giọt đơn lẻ và giọt sáp nhập mạnh trong mây dày $H = 4\text{ km}$.

2. Tính bán kính và vận tốc ở đáy mây: Áp dụng công thức (4) với $r_0 = 0{,}5\text{ mm}$, $\alpha = 4{,}0 \cdot 10^{-6}$ và $H = 4000\text{ m}$:

$$ \Delta r = \frac{\alpha H}{4} = \frac{4{,}0 \cdot 10^{-6} \times 4000}{4} = 4{,}0 \cdot 10^{-3}\text{ m} = 4{,}0\text{ mm} \implies r(H) = r_0 + \Delta r $$

✅ Bài giải chi tiết & Đáp số chuẩn

Ý 4.1. Ước lượng nồng độ thể tích $\alpha$ và cơ chế chi phối:

Độ nước của đám mây $w$ (khối lượng nước ngưng tụ trong một đơn vị thể tích không khí) liên hệ với nồng độ thể tích $\alpha$ bởi:

$$ w = \alpha \rho \implies \alpha = \frac{w}{\rho} $$

Do $w$ biến thiên từ cỡ $10^{-2}\text{ g/m}^3$ đến $4{,}0\text{ g/m}^3$, nên $\alpha$ biến thiên từ $10^{-8}$ đến $\alpha_{\max} = \frac{4{,}0 \cdot 10^{-3}\text{ kg/m}^3}{1000\text{ kg/m}^3} = 4{,}0 \cdot 10^{-6}$.
So sánh hai số hạng ở mẫu số của công thức (7) ngay tại $\alpha_{\max} = 4{,}0 \cdot 10^{-6}$:

$$ 7\alpha_{\max} = 2{,}8 \cdot 10^{-5} \ll \frac{3C\rho_0}{\rho} = \frac{3 \times 0{,}50 \times 1{,}2}{1000} = 1{,}8 \cdot 10^{-3} $$

Tỷ số $\frac{\alpha \rho}{C\rho_0} \approx 6{,}7 \cdot 10^{-3} \ll 1$. Điều này chứng tỏ trong mọi đám mây thực tế, lực cản khí động của không khí đóng vai trò chi phối lớn hơn rất nhiều so với lực hãm phản lực do sáp nhập vi giọt.

Ý 4.2. Chuyển động nhanh dần đều trong đám mây dày và giới hạn phân rã:

a) Với đám mây dày $H = 4{,}0\text{ km} = 4000\text{ m}$ và $\alpha = 4{,}0 \cdot 10^{-6}$, tỷ số độ tăng bán kính là $\frac{\alpha H}{4r_0} = \frac{4{,}0 \cdot 10^{-6} \times 4000}{4 \times 0{,}5 \cdot 10^{-3}} = 8 \gg 1$. Khi $\frac{\alpha x}{4r_0} \gg 1$, số hạng lũy thừa nghịch đảo trong (7) tắt dần về $0$ và $r_0\left(1 + \frac{\alpha x}{4r_0}\right) \approx \frac{\alpha x}{4}$. Công thức (7) rút gọn thành:

$$ v^2(x) \approx \frac{2gx}{7 + \frac{3C\rho_0}{\alpha \rho}} \implies a_{\text{hd}} = \frac{g}{7 + \frac{3C\rho_0}{\alpha \rho}} \approx \frac{4{,}0 \cdot 10^{-6} \times 1000}{3 \times 0{,}50 \times 1{,}2}\,g \approx 0{,}0022g $$

Giải thích vật lý: Vì sao có lực cản không khí mà giọt vẫn chuyển động nhanh dần đều? Bởi vì trọng lực tăng theo bậc ba của bán kính ($mg \propto r^3$), trong khi lực cản khí động chỉ tăng theo bậc hai của bán kính ($F_1 \propto r^2 v^2$). Do bán kính $r(x) \approx \frac{\alpha x}{4}$ tăng tuyến tính theo quãng đường $x$, nên tại mỗi độ sâu $x$, vận tốc giọt mưa luôn bám sát vận tốc giới hạn tức thời ứng với bán kính $r(x)$ tại điểm đó:

$$ v^2(x) \approx \frac{2g(\alpha x / 4)}{\frac{7}{4}\alpha + \frac{3C\rho_0}{4\rho}} \approx \frac{2g r(x)}{\frac{3C\rho_0}{4\rho}} = \frac{8\rho g r(x)}{3C\rho_0} $$

Vì $r(x) \propto x$, nên $v^2(x) \propto r(x) \propto x$, tạo ra một chuyển động nhanh dần đều rất "êm" với gia tốc $a_{\text{hd}} \approx 0{,}0022g$, không hề vượt qua giới hạn của công thức (1)!

b) Khi đi hết bề dày đám mây $H = 4{,}0\text{ km}$, bán kính giọt tăng thêm $\Delta r = \frac{\alpha H}{4} = 4{,}0\text{ mm}$, đạt giá trị tổng cộng:

$$ r(H) = r_0 + \frac{\alpha H}{4} = 0{,}5\text{ mm} + 4{,}0\text{ mm} = 4{,}5\text{ mm} $$

Vận tốc lý thuyết theo (10) tại đáy mây là $v(H) = \sqrt{2 a_{\text{hd}} H} \approx \sqrt{2 \times 0{,}0022 \times 9{,}8 \times 4000} \approx 13{,}1\text{ m/s}$. Tuy nhiên, theo tiêu chuẩn ổn định bề mặt (3), ứng với bán kính $r = 4{,}5\text{ mm}$, vận tốc tới hạn mà giọt bắt đầu mất dạng cầu và vỡ vụn là:

$$ v_{\text{cr}}(4{,}5\text{ mm}) = \sqrt{\frac{4\sigma}{C\rho_0 r}} = \sqrt{\frac{4 \times 7{,}5 \cdot 10^{-2}}{0{,}50 \times 1{,}2 \times 4{,}5 \cdot 10^{-3}}} \approx 10{,}5\text{ m/s} \quad (\text{hoặc } 10{,}2\text{ m/s nếu } \sigma = 0{,}070\text{ N/m}) $$

Do đó, giọt mưa không thể đạt tới $13\text{ m/s}$ mà sẽ bị biến dạng và phân rã thành các giọt nhỏ hơn ngay khi vận tốc vượt ngưỡng $\sim 10{,}2\text{–}10{,}5\text{ m/s}$.

Ý 4.3. Ước lượng mật độ số lượng vi giọt $n$ và khoảng cách trung bình $d$:

Thể tích của một vi giọt sương bán kính $r_{\text{s}} \approx 10\text{ }\mu\text{m} = 10^{-5}\text{ m}$ là $V_1 = \frac{4}{3}\pi r_{\text{s}}^3 \approx 4{,}19 \cdot 10^{-15}\text{ m}^3$. Với nồng độ thể tích $\alpha \le 4{,}0 \cdot 10^{-6}$, mật độ số lượng vi giọt trong mây là:

$$ n = \frac{\alpha}{\frac{4}{3}\pi r_{\text{s}}^3} \sim 10^7\text{–}10^8\text{ m}^{-3} \quad (\text{ứng với } r_{\text{s}} \text{ từ } 10\text{ }\mu\text{m đến } 40\text{ }\mu\text{m}) $$

Khoảng cách trung bình giữa các vi giọt lân cận là:

$$ d \approx n^{-1/3} \sim \left(10^7\text{–}10^8\text{ m}^{-3}\right)^{-1/3} \approx 2\text{–}5\text{ mm} \text{ (đến cỡ vài xentimét khi } \alpha \text{ nhỏ)} $$

Đáp số chuẩn Chặng 4 (Ý 4.2b): Bán kính lý thuyết của giọt mưa ở cuối đám mây dày $4{,}0\text{ km}$ là $r(H) = 4{,}5\text{ mm}$.

🚀 Thử thách mở rộng Chặng 4: Ảnh hưởng của dòng khí thăng (Updraft) trong đám mây dông

Bài toán mở rộng: Trong các đám mây tích vũ thực tế, luôn có một luồng không khí ấm bốc lên thẳng đứng với vận tốc $u_0$ (dòng thăng, thường cỡ $2\text{–}8\text{ m/s}$). Giả sử các vi giọt sương rất nhỏ nên bị cuốn đi cùng dòng khí với vận tốc hướng lên là $u_0$. Hãy viết lại điều kiện để một giọt mưa bán kính $r$ có thể bắt đầu rơi xuống đối với mặt đất, và giải thích vì sao mưa rào từ các đám mây có dòng thăng mạnh luôn gồm các hạt mưa rất to!

Hướng dẫn: Vận tốc của giọt mưa đối với luồng khí là $v_{\text{tđ}} = v_{\text{đất}} + u_0$. Để giọt mưa rơi xuống mặt đất ($v_{\text{đất}} > 0$), vận tốc rơi giới hạn của nó đối với không khí phải thắng được vận tốc dòng thăng: $v_{\text{gh}}(r) > u_0$. Thay $v_{\text{gh}}(r) \approx \sqrt{\frac{8\rho g r}{3C\rho_0}} > u_0 \implies r > \frac{3C\rho_0 u_0^2}{8\rho g}$. Với dòng thăng $u_0 = 6\text{ m/s}$, chỉ những giọt đã sáp nhập đạt bán kính $r > 0{,}83\text{ mm}$ mới rơi xuống được mặt đất!

📌 Tổng kết mô hình vật lý của giọt mưa:

Mô hình giải tích trên đã kết nối trọn vẹn ba cơ chế cốt lõi (trọng trường, lực cản khí động học chất lưu và động lực học khối lượng biến thiên). Bảng dưới đây tóm tắt các chế độ chuyển động đặc trưng của hạt mưa:

Chế độ vật lý	Điều kiện tham số	Quy luật vận tốc $v^2(x)$	Đặc trưng động lực học
Rơi tự do (Galileo)	$C \to 0,\; \alpha \to 0$	$v^2 = 2gx$	Gia tốc không đổi $a = g$
Sương tĩnh (Sommerfeld)	$C = 0,\; \alpha x \gg 4r_0$	$v^2 = \dfrac{2}{7}gx$	Nhanh dần đều với $a = \dfrac{g}{7}$
Giọt đơn lẻ trong khí quyển	$\alpha \to 0,\; C = 0{,}5$	$v_{\text{gh}}^2\left(1 - e^{-x/L}\right)$	Bão hòa nhanh sau $L \approx 11\text{ m}$
Sáp nhập trong mây dày	$\alpha x \gg 4r_0,\; C = 0{,}5$	$v^2 \approx \dfrac{2gx}{7 + 3C\rho_0/(\alpha\rho)}$	Tăng tốc êm ($a_{\text{hd}} \approx 0{,}002g$) tới ngưỡng vỡ $v_{\text{cr}} \approx 10\text{ m/s}$

Ghi chú mở rộng: Ngoài ba cơ chế chính trên, vật lý đám mây hiện đại còn xét thêm quá trình ngưng tụ hơi nước từ tâm mầm, tương tác tĩnh điện giữa các giọt mang điện trái dấu, sự bay hơi của giọt khi rơi qua lớp không khí khô bên dưới đám mây và sự kết tinh băng ở tầng đối lưu trên cao.