
Mô hình tính toán mô phỏng quá trình bất ổn định Kelvin-Helmholtz
lượt xem 1
download

Trong bài báo này chúng tôi dùng mô hình tính toán dòng chảy hai pha để mô phỏng sự hình thành và phát triển của sự bất ổn định Kelvin-Helmholtz. Bài toán được xét trên một kênh hở tuần hoàn theo hướng của trục Ox . Kích thước của kênh là 3m chiều dài và 3m chiều cao, lớp trên là một lớp nước và lớp dưới là một lớp bùn cát.
Bình luận(0) Đăng nhập để gửi bình luận!
Nội dung Text: Mô hình tính toán mô phỏng quá trình bất ổn định Kelvin-Helmholtz
- Tuyển tập Hội nghị Khoa học thường niên năm 2015. ISBN : 978-604-82-1710-5 MÔ HÌNH TÍNH TOÁN MÔ PHỎNG QUÁ TRÌNH BẤT ỔN ĐỊNH KELVIN-HELMHOLTZ Nguyễn Đức Hậu1, Nguyễn Thị Lý1, Sylvain Guillou2 và Kim Dan Nguyen3 Đại học Thủy lợi, email: ndhau.dhtl@tlu.edu.vn 2 2 Université de Caen, email: sylvain.guillou@unicaen.fr 3 Saint-Venant Laboratory for Hydraulics, email: dan.nguyen@saint-venant-lab.fr 1. GIỚI THIỆU CHUNG 0 1000kg / m3 , tỉ khối của hạt chất rắn là Trong bài báo này chúng tôi dùng mô hình 1 1360kg / m và độ nhớt động lực 3 tính toán dòng chảy hai pha để mô phỏng sự 0 10Pas . Ở điều kiện ban đầu áp đặt vận hình thành và phát triển của sự bất ổn định tốc nằm ngang u ( z ) của dòng chảy và nồng Kelvin-Helmholtz. Bài toán được xét trên độ khối ( z ) . Chúng được xác định bởi một kênh hở tuần hoàn theo hướng của trục (xem hình 1). Ox . Kích thước của kênh là 3m chiều dài và 3m chiều cao, lớp trên là một lớp nước và lớp dưới là một lớp bùn cát. Ở điều kiện ban đầu áp đặt vận tốc nằm ngang của dòng chảy. Trường vận tốc có một điểm uốn tại giao diện của nước và bùn cát. Mô hình được tính toán trên lưới chia đều với số nút lưới là 301 301 . 2. PHƯƠNG PHÁP NGHIÊN CỨU (a) Mô hình hai pha Trong mô hình dòng chảy hai pha chúng tôi sử dụng mô hình Euler-Euler để biểu diễn các hệ phương trình dòng chảy của các hạt lỏng và các hạt rắn. Xuất phát từ hệ phương trình Navier-Stokes, chúng tôi dùng các công thức của Drew & Lahey [2] để xây dựng các phương trình bảo toàn khối lượng và phương trình chuyển động. b) Cấu hình bài toán nghiên cứu Bài toán được tính toán trên một kênh hở tuần hoàn theo hướng của trục Ox . Kích thước của kênh được xác định là hình vuông với 3m chiều dài và 3m chiều cao, lớp trên là một lớp nước và lớp dưới là một lớp bùn cát. Tỉ khối của nước được xác định là Hình 1. Trường vận tốc và nồng độ khối 97
- Tuyển tập Hội nghị Khoa học thường niên năm 2015. ISBN : 978-604-82-1710-5 zh u( z) U 0 F , ; zh ( z ) 0 1 F , ; 1 y 1 erf , y 0 F , y 1 1 1 erf y ở đó erf là “error function”, là độ dày của phần mà vận tốc có thay đổi, h là độ dày của lớp cát phía dưới, là hằng số. Trong bài toán này ta chọn h 1m , 0.3m và 0.1 . 3. KẾT QUẢ NGHIÊN CỨU Trong [5] chúng tôi đã nghiên cứu với trường hợp lưới thô 101101 nút lưới ở đó chúng ta đã thấy rằng sự bất ổn định Kelvin-Helmholtz chỉ xuất hiện trong trường hợp số Richardson Ri nhỏ hơn 0.25 điều này phù hợp với các kết quả nghiên cứu trước đó của Milles [3]. Trong [6] chúng tôi đã thực hiện mô phỏng với bốn loại lưới thô, trung bình, lưới mịn và lưới siêu mịn. So sánh các kết quả mô phỏng nhận được đối với các loại lưới trong bài báo này chúng tôi lựa chọn lưới mịn (F) với 301 301 nút lưới để nghiên cứu. Trong mô phỏng này bước thời gian được chọn là 5 104 (s), vận tốc U 0 2.5 m/s tương ứng với hằng số Richardson Ri 0.113 . Đường kính của hạt chất rắn được mô phỏng là 50 m (cát mịn). Hình 2 cho chúng ta thấy đường mức của nồng độ khối mô phỏng tại các thời điểm khác nhau. Từ thời điểm t 17 s ở bề mặt xoáy bắt đầu được hình thành. Chú ý rằng trong các thời điểm trước đó (trước t 17 s) ở bề mặt vẫn còn ổn định. Sự bất ổn định được tiếp tục hình thành và phát triển lớn lên trong khoảng thời gian từ 17s đến 20.5s. Xoáy phát triển lớn nhất tại thời điểm 20.5s. Hình 2. Nồng độ khối 98
- Tuyển tập Hội nghị Khoa học thường niên năm 2015. ISBN : 978-604-82-1710-5 Hình 3 cho chúng ta thấy đường mức của dày, các vùng “noyau” và “périphérique” vẫn độ nhớt mô phỏng tại các thời điểm ứng với chưa phát triển hoàn toàn. Chúng ta có thể hình 2. Bắt đầu từ thời điểm t 18 s độ nhớt quan sát khá rõ 3 vùng đó tại thời điểm 20s: bắt đầu xuất hiện tuy nhiên chúng ta còn khó Vùng “cordon” đã mảnh hơn và hai vùng có thể phân biệt ba vùng “noyau”, “noyau” và “périphérique” phát triển mạnh “périphérique” và “cordon”. Ở thời điểm 19s mẽ. Tại thời điểm t= 20.5s các vùng phát vùng “cordon” đã xuất hiện nhưng vẫn còn triển hoàn toàn và được quát sát rõ nét nhất. Hình 3. Độ nhớt 4. KẾT LUẬN [3] Milles, J. W. (1961), On the stability of heterogeneous shear flows, J. Fluid Mech, Nhờ mô hình tính toán dòng chảy hai pha 10, 496-508. chúng tôi đã mô phỏng sự hình thành và phát [4] Nguyen D.H., Guillou S., Nguyen K.D., triển của quá trình bất ổn định Kelvin- Pham Van Bang D., Chauchat J. (2012), Helmholtz xuất hiện ở bề mặt bùn cát khi tốc Simulation of dredged sediment releases độ dòng đủ lớn (số Richardson Ri nhỏ hơn into homogeneous water using a two-phase 0.25 ) và có điểm uốn tại bề mặt. Trong các model. Advances in Water Resources 48 nghiên cứu tiếp theo chúng ta còn có thể xem (2012) 102–112. xét sự ảnh hưởng của kích thước hạt chất rắn http://dx.doi.org/10.1016/j.advwatres.2012. đến sự hình thành và phát triển sự bất ổn định 03.009. Kelvin-Helmholtz. [5] Nguyễn Đức Hậu, Nguyễn Thị Lý, Sylvain Guillou, Kim Dan Nguyen (2013), Nghiên 5. TÀI LIỆU THAM KHẢO cứu sự bất ổn định Kelvin-Helmholtz ở bề [1] Caulfield, C. P., and W. R. Peltier (2000), mặt bùn cát, Tuyển tập hội nghị thường The anatomy of the mixing transition in niên 2013 của trường Đại học Thủy lợi. homogeneous and stratified free shear layers, [6] Nguyễn Đức Hậu, Nguyễn Thị Lý, Sylvain Journal of Fluid Mechanics, 413, 1-47. Guillou, Kim Dan Nguyen (2014), Nghiên [2] Drew D.A. and Lahey R.T.: Analytical cứu sự ảnh hưởng của lưới trong mô hình Modelling of Multiphase Flow, in Particulate tính toán sự bất ổn định Kelvin-Helmholtz, Two-Phase Flow, ed. M.C. Roco, Tuyển tập hội nghị thường niên 2014. Butterworth-Heinemann, Boston, (1993). Trường Đại học Thủy lợi. 99
- Tuyển tập Hội nghị Khoa học thường niên năm 2015. ISBN : 978-604-82-1710-5 100

CÓ THỂ BẠN MUỐN DOWNLOAD
-
MÔ PHỎNG MONTE CARLO
8 p |
419 |
98
-
Bài giảng Mô hình hóa môi trường: Chương 4 - GV. Trương Thị Thu Hương
14 p |
226 |
49
-
MÔ HÌNH CHẤT LƯỢNG NƯỚC - CHƯƠNG 5
4 p |
176 |
29
-
Bài giảng Quy hoạch tuyến tính: Chương 3 - ThS. Nguyễn Văn Phong (2016 - BT)
17 p |
237 |
16
-
Bài giảng Quy hoạch tuyến tính: Chương 2 - ThS. Nguyễn Văn Phong (2016 - BT)
12 p |
156 |
7
-
Bài giảng Mô hình hồ chứa - Sông ngoài
72 p |
83 |
6
-
Dự báo biến động đáy biển, đường bờ từ hoạt động của dự án “Nạo vét, khơi thông luồng hàng hải khu kinh tế Vân Phong, kết hợp tận thu cát nhiễm mặn để xuất khẩu”
8 p |
56 |
3
-
Tính toán nước dâng do bão các tỉnh ven biển khu vực Vịnh Bắc Bộ
29 p |
55 |
3
-
Nghiên cứu vận động khối nước ô nhiễm kênh vùng triều ứng với trường hợp thay đổi lưu lượng nguồn và vị trí đặt kênh
11 p |
29 |
2
-
Đánh giá khả năng sử dụng dữ liệu mưa vệ tinh để mô phỏng dòng chảy bằng mô hình thủy văn HYPE, áp dụng cho lưu vực sông Sê San
10 p |
9 |
1
-
Tính độ dẫn điện trong môi trường vết nứt bão hòa sử dụng phân bố fractal
3 p |
1 |
1
-
Tính phổ sóng điều hòa bậc cao của phân tử CO tương tác với laser hồng ngoại sóng trung bằng phương pháp TDSE
9 p |
49 |
1
-
Nghiên cứu mô phỏng nước dâng do gió mùa khu vực biển Tây Nam Bộ
14 p |
5 |
1
-
Nghiên cứu xác định thông số nhiệt động học sản phẩm nổ nhiên liệu metan trong động cơ nổ
10 p |
3 |
1
-
Ứng dụng mô hình HEC - RESSIM tính toán điều tiết hệ thống hồ chứa trên lưu vực sông Ba
3 p |
2 |
1
-
Nghiên cứu chế độ thủy động lực vùng cửa Sông Gianh, tỉnh Quảng Bình
9 p |
12 |
1
-
Đánh giá ngập lụt hạ lưu do vỡ đập hồ Đồng Mỏ, tỉnh Vĩnh Phúc
3 p |
3 |
1


Chịu trách nhiệm nội dung:
Nguyễn Công Hà - Giám đốc Công ty TNHH TÀI LIỆU TRỰC TUYẾN VI NA
LIÊN HỆ
Địa chỉ: P402, 54A Nơ Trang Long, Phường 14, Q.Bình Thạnh, TP.HCM
Hotline: 093 303 0098
Email: support@tailieu.vn
