Bốc hơi là quá trình chuyển pha nhiệt động lực học trong đó các phân tử chất lỏng chuyển sang trạng thái khí tại bề mặt thoáng tự do, xảy ra ở mọi nhiệt độ dưới điểm sôi xác định của khối chất lỏng đó. Đây là một cơ chế cơ bản trong chu trình tuần hoàn vật chất và truyền nhiệt, chi phối sự trao đổi ẩm giữa bề mặt hành tinh và khí quyển cũng như các phản ứng công nghệ trong đời sống. Bài viết này trình bày chi tiết về bản chất vi mô, các quy luật nhiệt động lực học, các công thức định lượng kinh điển, phương pháp quan trắc tiêu chuẩn và phạm vi ứng dụng thực tiễn của hiện tượng bốc hơi.
Khái niệm và bản chất vật lý của bốc hơi
Bốc hơi (evaporation) là một dạng hóa hơi (vaporization) bề mặt đặc thù. Khác với hiện tượng sôi hình thành bọt khí trong lòng chất lỏng khi áp suất hơi bão hòa cân bằng với áp suất ngoài, bốc hơi chỉ giới hạn tại lớp phân tử mỏng giáp ranh giữa pha lỏng và pha khí. Trong một khối chất lỏng ở trạng thái cân bằng nhiệt, động năng của các phân tử không phân bố đồng đều mà tuân theo phân bố vận tốc Maxwell-Boltzmann. Hầu hết các phân tử có động năng trung bình tương ứng với nhiệt độ của khối chất lỏng, nhưng luôn tồn tại một tỷ lệ phân tử ở vùng đuôi phân bố mang động năng chuyển động nhiệt vượt trội.
Khi các phân tử giàu năng lượng này di chuyển tới gần bề mặt thoáng với vector vận tốc hướng ra ngoài, chúng có đủ năng lượng để thắng được lực hút liên phân tử (lực van der Waals hoặc liên kết hydro) của các phân tử lân cận trong pha lỏng. Quá trình vượt qua rào cản năng lượng này cho phép phân tử thoát ly khỏi bề mặt lỏng và hòa nhập vào pha hơi. Vì các phân tử mang động năng cao nhất liên tục thoát đi, động năng trung bình của các phân tử còn lại trong pha lỏng bị suy giảm, dẫn đến hiện tượng hạ nhiệt độ khối chất lỏng nếu không có nguồn nhiệt bổ sung từ môi trường bên ngoài.
Cơ chế vi mô và nhiệt động lực học pha
Xét theo quan điểm nhiệt động học cổ điển, quá trình bốc hơi đòi hỏi chất lỏng phải hấp thụ một lượng năng lượng nhiệt tương ứng với ẩn nhiệt hóa hơi (nhiệt ẩn chuyển pha, latent heat of vaporization). Ẩn nhiệt này cung cấp công để phá vỡ các liên kết kết tụ giữa các phân tử trong pha lỏng và cung cấp công dãn nở thể tích chống lại áp suất khí quyển khi chuyển từ thể lỏng đậm đặc sang thể hơi loãng.
Theo tiêu chuẩn quốc tế IAPWS ban hành qua công trình của Wagner và Pruß (2002), đối với nước tinh khiết ở nhiệt độ sôi tiêu chuẩn 100 °C (tương đương 373.15 K), entanpi hóa hơi đo được là 2257 kJ/kg (tương ứng với giá trị mol xấp xỉ 40.66 kJ/mol). Ở nhiệt độ môi trường thấp hơn, khoảng cách liên kết và cường độ liên kết hydro giữa các phân tử nước lớn hơn, đòi hỏi năng lượng hóa hơi cao hơn. Cụ thể, theo số liệu thực nghiệm chuẩn của Wagner và Pruß (2002), tại 20 °C (tương đương 293.15 K), nhiệt hóa hơi của nước đạt 2454 kJ/kg, và tại 0 °C giá trị này tăng lên mức 2501 kJ/kg.
Tại ranh giới tiếp xúc giữa pha lỏng và pha khí, quá trình bay hơi luôn diễn ra đồng thời với quá trình ngưng tụ ngược lại. Nếu không gian phía trên chất lỏng bị giới hạn trong bình kín, số lượng phân tử hơi quay trở lại pha lỏng sẽ tăng dần cho đến khi tốc độ ngưng tụ bằng đúng tốc độ bốc hơi. Khi đó, hệ thống đạt trạng thái cân bằng động nhiệt động lực học, và áp suất riêng phần của hơi ở trạng thái này được gọi là áp suất hơi bão hòa.
Các yếu tố chi phối tốc độ bốc hơi
Tốc độ bốc hơi từ một bề mặt thoáng phụ thuộc đồng thời vào các thông số nội tại của chất lỏng và điều kiện môi trường khí quyển bao quanh:
- Nhiệt độ bề mặt chất lỏng: Khi nhiệt độ tăng, động năng trung bình của các phân tử tăng lên, làm tăng phân số phân tử có động năng vượt rào cản năng lượng thoát ly, đồng thời làm tăng áp suất hơi bão hòa tại bề mặt tiếp xúc.
- Độ ẩm tương đối của không khí: Mức chênh lệch giữa áp suất hơi bão hòa tại nhiệt độ bề mặt và áp suất riêng phần của hơi nước trong không khí là động lực chủ đạo cho quá trình khuếch tán phân tử. Không khí càng khô, tốc độ bay hơi càng cao; khi độ ẩm đạt ngưỡng bão hòa, bốc hơi ròng sẽ triệt tiêu.
- Diện tích bề mặt tiếp xúc: Do bốc hơi là hiện tượng bề mặt thuần túy, tốc độ bốc hơi trên toàn khối chất lỏng tỷ lệ thuận với diện tích bề mặt thoáng tiếp xúc trực tiếp với pha khí.
- Chế độ chuyển động của không khí (gió): Sự lưu thông của không khí giúp quét sạch lớp biên bão hòa ẩm tích tụ ngay sát bề mặt chất lỏng, duy trì građien áp suất hơi cao giữa bề mặt và tầng khí quyển tự do.
- Bản chất hóa học của chất lỏng: Các chất lỏng có lực tương tác liên phân tử yếu (như ether, ethanol, acetone) có áp suất hơi bão hòa cao và nhiệt hóa hơi thấp, dẫn đến tốc độ bốc hơi nhanh hơn đáng kể so với nước ở cùng điều kiện nhiệt độ.
Mô hình toán học và công thức định lượng
Để tính toán tốc độ truyền khối trong quá trình bốc hơi, giới khoa học đã phát triển các mô hình từ quy mô động học phân tử đến mô hình cân bằng năng lượng khí tượng học.
Định luật Dalton và phương trình Hertz-Knudsen
Ở cấp độ cơ bản, John Dalton thiết lập định luật bán thực nghiệm cho thấy tốc độ bốc hơi tỷ lệ thuận với độ hụt áp suất hơi. Biểu thức của tốc độ bốc hơi theo định luật Dalton được mô tả dạng:
trong đó là tốc độ bốc hơi trên một đơn vị diện tích bề mặt, là áp suất hơi bão hòa tại nhiệt độ bề mặt chất lỏng, là áp suất riêng phần của hơi nước trong không khí xung quanh, và là hệ số phụ thuộc vào vận tốc gió và độ nhám bề mặt.
Ở quy mô vi mô trong điều kiện chân không hoặc áp suất thấp, thông lượng bốc hơi ròng được mô tả bằng phương trình kinh điển Hertz-Knudsen dựa trên thuyết động học phân tử:
trong đó là thông lượng mol ròng đi qua bề mặt chất lỏng (tính bằng mol trên mét vuông giây), là hệ số bốc hơi, là hệ số ngưng tụ, là áp suất bão hòa tương ứng với nhiệt độ bề mặt lỏng, là áp suất riêng phần của pha hơi, là khối lượng mol của chất, là hằng số khí lý tưởng và là nhiệt độ nhiệt động lực học tuyệt đối tính bằng Kelvin.
Tuy nhiên, công trình tổng quan của Persad và Ward (2016) trên tạp chí Chemical Reviews đã chỉ ra rằng phương trình Hertz-Knudsen cổ điển có nhược điểm lớn khi coi các hệ số là tham số tùy biến. Khảo sát thực nghiệm trên 127 bộ dữ liệu thực nghiệm bốc hơi dừng của nước và ethanol cho thấy hệ số bốc hơi từng dao động qua 3 bậc độ lớn trong các tài liệu cũ do bỏ qua hiện tượng bước nhảy nhiệt độ tại mặt phân giới và sự ghép pha nhiệt động lực học.
Phương trình Penman cho bốc hơi tự nhiên
Trong khí tượng thủy văn, phương trình do Penman (1948) đề xuất là cột mốc đột phá khi tích hợp cả định luật bảo toàn năng lượng bề mặt và động lực học khuếch tán khí quyển. Dạng tổng quát của phương trình Penman được biểu diễn:
trong đó là lượng bốc hơi tương đương từ bề mặt nước thoáng mở, là độ dốc của đường cong áp suất hơi bão hòa theo nhiệt độ, là bức xạ thuần tại bề mặt tính theo đương lượng bốc hơi, là hằng số tâm trắc của khí quyển, và là thông số khí động học đặc trưng cho khả năng vận chuyển hơi nước của gió.
Động học bốc hơi của giọt chất lỏng bề mặt
Đối với các giọt chất lỏng nằm trên bề mặt chất rắn (sessile droplet), quá trình bốc hơi chịu sự kiểm soát chặt chẽ của hiện tượng thấm ướt và góc tiếp xúc. Theo Cazabat và Guéna (2010), quá trình bốc hơi giọt vĩ mô trên bề mặt rắn diễn ra qua hai chế độ hình học giới hạn đặc trưng:
| Chế độ bốc hơi | Đặc điểm hình học | Cơ chế động học | Hiện tượng điển hình |
|---|---|---|---|
| Ghim bán kính tiếp xúc (CCR) | Bán kính tiếp xúc không đổi, góc tiếp xúc giảm dần theo thời gian bốc hơi | Dòng chảy mao dẫn xuyên tâm hướng ra ngoài từ tâm giọt ra chu vi biên để bù đắp lượng chất lỏng bốc hơi nhanh ở rìa | Hiệu ứng vết cà phê (coffee-ring effect) lắng đọng hạt keo |
| Góc tiếp xúc không đổi (CCA) | Góc tiếp xúc giữ nguyên, đường chu vi tiếp xúc co dần về tâm giọt | Đường ba pha lỏng-khí-rắn chuyển động trượt tự do trên bề mặt chất nền lý tưởng nhẵn bóng | Sự co ngót đều đặn của giọt chất lỏng không bị ghim mép |
Nghiên cứu của Cazabat và Guéna (2010) chứng minh rằng ở vùng rìa của giọt chất lỏng, građien thế hóa học và thông lượng bốc hơi đạt cực đại do cấu trúc hình học tiếp xúc mỏng, thúc đẩy dòng đối lưu mao dẫn bù trừ từ tâm ra chu vi ngoài.
Phân biệt bốc hơi với sôi và thăng hoa
Để tránh nhầm lẫn giữa các quá trình chuyển thể pha trong nhiệt động lực học, cần phân định rõ ràng các ranh giới vật lý giữa bốc hơi, sôi và thăng hoa:
| Tiêu chí so sánh | Bốc hơi (Evaporation) | Sôi (Boiling) | Thăng hoa (Sublimation) |
|---|---|---|---|
| Pha ban đầu | Pha lỏng | Pha lỏng | Pha rắn |
| Vị trí chuyển pha | Chỉ diễn ra tại bề mặt thoáng tự do | Diễn ra đồng thời trong toàn bộ thể tích và bề mặt lỏng | Diễn ra tại bề mặt tiếp xúc của chất rắn |
| Điều kiện chuyển pha | Mọi nhiệt độ khi chất còn ở thể lỏng dưới điểm sôi | Tại nhiệt độ sôi xác định khi áp suất hơi bão hòa bằng áp suất ngoài | Áp suất môi trường thấp hơn áp suất điểm ba pha và nhiệt độ dưới nhiệt độ điểm ba pha của chất |
| Dấu hiệu trực quan | Không xuất hiện bọt khí bên trong chất lỏng | Xuất hiện bọt khí chuyển động sôi nổi từ đáy và lòng chất lỏng | Chất rắn biến mất dần thành thể khí không qua thể lỏng |
| Áp suất hơi bão hòa | Nhỏ hơn áp suất khí quyển bao quanh | Bằng hoặc lớn hơn áp suất khí quyển đè lên bề mặt | Đạt trạng thái cân bằng hơi trực tiếp với pha rắn |
Phương pháp quan trắc và đo đạc thực nghiệm
Trong mạng lưới quan trắc khí tượng thủy văn, đo đạc lượng bốc hơi bề mặt là yêu cầu bắt buộc để đánh giá cán cân nước vùng và nhu cầu tưới tiêu nông nghiệp. Theo tài liệu tổng quan chuyên ngành của Verhoef và Campbell (2005) trong Encyclopedia of Hydrological Sciences, các phương pháp quan trắc lượng bốc hơi bề mặt tiêu chuẩn quốc tế chủ yếu vận hành dựa trên các dụng cụ sau:
- Ống đo bốc hơi Piche: Được lắp đặt cố định bên trong lều khí tượng ở độ cao chuẩn. Dụng cụ gồm một ống thủy tinh có vạch chia chính xác đến 0.1 mm, đầu dưới gắn giấy thấm chuyên dụng có diện tích bốc hơi xác định, dùng để đo sức bốc hơi của không khí trong điều kiện không chịu tác động trực tiếp của bức xạ mặt trời.
- Thùng bốc hơi Class A: Được lắp đặt ngoài trời trên giá gỗ cách mặt đất để không khí lưu thông bên dưới đáy thùng. Theo Verhoef và Campbell (2005), thùng bốc hơi Class A hình trụ tròn có đường kính tiêu chuẩn 1206 mm và chiều cao thành thùng 254 mm, làm bằng vật liệu thép không gỉ hoặc kim loại mạ kẽm. Mực nước trong thùng được theo dõi định kỳ bằng thước đo đầu nhọn (hook gauge) hoặc giếng đo tĩnh để tính toán lượng nước bốc hơi ròng theo đơn vị milimét.
Vai trò sinh thái, nghiên cứu tại Việt Nam và điều kiện hạn chế
Bốc hơi giữ vị trí trung tâm trong chu trình thủy văn toàn cầu. Theo nghiên cứu tổng quan của Oki và Kanae (2006) công bố trên tạp chí Science, tổng lượng nước bốc hơi trên toàn hành tinh đạt khoảng 502000 km3/năm. Trong cấu trúc này, các đại dương bốc hơi một lượng khổng lồ khoảng 436500 km3/năm (chiếm xấp xỉ 87% tổng thông lượng), đóng vai trò cỗ máy nhiệt khổng lồ hấp thụ năng lượng mặt trời và chuyển tải ẩm vào lục địa. Bề mặt đất liền đóng góp lượng bốc thoát hơi nước khoảng 65500 km3/năm, nuôi dưỡng hệ thống sông ngòi và thảm thực vật sinh quyển.
Tại Việt Nam, sự phân hóa mạnh mẽ của lượng bốc thoát hơi nước theo mùa gió mùa và địa hình ảnh hưởng sâu sắc đến chế độ hạn hán. Nghiên cứu của Le và cộng sự (2020) công bố trên Journal of Hydrology: Regional Studies đã phân tích chuỗi dữ liệu 30 năm trên toàn lãnh thổ Việt Nam bằng cách sử dụng chỉ số bốc thoát hơi nước chuẩn hóa SPEI kết hợp dữ liệu tái phân tích MERRA-2. Công trình này làm rõ biến động bốc thoát hơi nước tại 7 vùng sinh thái, chứng minh rằng sự gia tăng nhiệt độ làm gia tăng mạnh mẽ lượng bốc thoát hơi tiềm năng, làm trầm trọng thêm mức độ gay gắt của các đợt hạn hán lịch sử ở khu vực Nam Trung Bộ và Tây Nguyên.
Trong kỹ thuật công nghiệp, nguyên lý làm mát bằng bốc hơi (evaporative cooling) được khai thác tối đa trong tháp giải nhiệt của các nhà máy điện và hệ thống điều hòa không khí. Công nghệ sấy phun, cô đặc dung dịch và chưng cất chân không đều vận hành dựa trên việc điều khiển tốc độ bốc hơi bề mặt qua kiểm soát nhiệt độ và áp suất buồng công nghệ.
Tuy vậy, việc áp dụng các mô hình bốc hơi có những ranh giới hạn chế rõ rệt. Các công thức kinh điển như Dalton hay Penman chỉ áp dụng chính xác cho bề mặt nước thoáng rộng lớn và đồng nhất. Khi áp dụng cho vùng địa hình chia cắt phức tạp, hồ chứa nước có dòng chảy nhiệt phân tầng hoặc khu vực đô thị có hiệu ứng đảo nhiệt, các giả định về lớp biên khí quyển phẳng không còn duy trì, đòi hỏi phải hiệu chỉnh bằng phương pháp đo tương quan xoáy (eddy covariance) hoặc mô hình truyền nhiệt đối lưu phi tuyến tính phức tạp.