Trao đổi nhiệt trong một số quá trình lan truyền của sóng trong môi trường hai pha

ĐẠI HỌC THÁI NGUYÊN TRƢỜNG ĐẠI HỌC KHOA HỌC PHẠM THỊ HẢI TRAO ĐỔI NHIỆT TRONG MỘT SỐ QUÁ TRÌNH LAN TRUYỀN CỦA SĨNG TRONG MƠI TRƢỜNG HAI PHA LUẬN VĂN THẠC SĨ TOÁN HỌC Thái Nguyên - 2014 Số hóa Trung tâm Học liệu – Đại học Thái Nguyên Tai ngay!!! Ban co the xoa dong chu nay!!! http://www.lrc-tnu.edu.vn/ ĐẠI HỌC THÁI NGUYÊN TRƢỜNG ĐẠI HỌC KHOA HỌC PHẠM PHẠM THỊ THỊ HẢI HẢI TRAO TRAO ĐỔI ĐỔI NHIỆT NHIỆT TRONG TRONG MỘT MỘT SỐ SỐ QUÁ QUÁ TRÌNH TRÌNH LAN LAN TRUYỀN TRUYỀN CỦA CỦA SĨNG SĨNG TRONG TRONG MÔI MÔI TRƢỜNG TRƢỜNG HAI HAI PHA PHA Chuyên Chuyên ngành: ngành: Toán Toán ứng ứng dụng dụng Mã Mã số: số: 60 60 46 46 01 01 12 12 LUẬN VĂN THẠC LUẬN VĂN THẠC SĨ SĨ TOÁN TOÁN HỌC HỌC Ngƣời Ngƣời hƣớng hƣớng dẫn dẫn khoa khoa học: học: PGS.TS PGS.TS Nguyễn Nguyễn Văn Văn Tuấn Tuấn Thái Thái Nguyên Nguyên 2014 2014 Số hóa Trung tâm Học liệu – Đại học Thái Nguyên http://www.lrc-tnu.edu.vn/ i LỜI CAM ĐOAN Tôi xin cam đoan kết có luận văn thân thực hướng dẫn thầy giáo PGS.TS Nguyễn Văn Tuấn Ngoài tài liệu tham khảo trích dẫn, số liệu kết có luận văn trung thực chưa có cơng bố cơng trình khác Thái Ngun, tháng năm 2014 Ngƣời thực Phạm Thị Hải Số hóa Trung tâm Học liệu – Đại học Thái Nguyên http://www.lrc-tnu.edu.vn/ ii MỤC LỤC Trang bìa phụ Lời cam đoan……………………………………………………………… i Mục lục… ii Danh mục ký hiệu, chữ viết tắt iv Danh mục hình vẽ v MỞ ĐẦU CHƢƠNG TỔNG QUAN CHƢƠNG CƠ SỞ LÝ THUYẾT 10 2.1 Hệ phương trình thuỷ - Nhiệt động lực học cho môi trường hai pha lỏng 10 2.1.1 Xây dựng mơ hình 10 2.1.2 Hệ phương trình sở 10 2.1.2.1 Hệ phương trình học hỗn hợp hai pha 11 2.1.2.2 Điều kiện đầu điều kiện biên 19 2.2 Phương pháp giải số chương trình tính tốn 20 2.2.1 Phương pháp giải số 20 2.2.2 Chương trình tính tốn 22 CHƢƠNG MỘT SỐ KẾT QUẢ NGHIÊN CỨU VỀ SỰ ẢNH HƢỞNG CỦA SĨNG LÊN Q TRÌNH TRAO ĐỔI NHIỆT GIỮA CÁC PHA TRONG HỖN HỢP CHẤT LỎNG CHỨA BỌT HƠI 23 3.1 Quá trình trao đổi nhiệt pha có sóng ngắn lan truyền hỗn hợp nước ni tơ lỏng 23 3.1.1 Quá trình trao đổi nhiệt pha hỗn hợp nước chứa bọt 23 3.1.1.1 Sự phụ thuộc vào cường độ ban đầu sóng xung kích 24 3.1.1.2 Sự phụ thuộc vào bán kính bọt 26 3.1.1.3 Sự phụ thuộc vào thể tích pha 28 3.1.2 Quá trình trao đổi nhiệt pha hỗn hợp ni tơ lỏng chứa bọt 30 3.1.2.1 Sự phụ thuộc vào cường độ ban đầu sóng xung kích 31 Số hóa Trung tâm Học liệu – Đại học Thái Nguyên http://www.lrc-tnu.edu.vn/ iii 3.1.2.2 Sự phụ thuộc vào bán kính bọt 33 3.1.2.3 Sự phụ thuộc vào thể tích pha 34 3.1.3 Nhận xét 38 3.2.1 Quá trình trao đổi nhiệt pha hỗn hợp nước chứa bọt 38 3.2.1.1 Sự phụ thuộc vào cường độ ban đầu sóng xung kích 39 3.2.1.2 Sự phụ thuộc vào bán kính bọt 41 3.2.1.3 Sự phụ thuộc vào thể tích pha 43 3.2.2 Quá trình trao đổi nhiệt pha hỗn hợp ni tơ chứa bọt 45 3.2.2.1 Sự phụ thuộc vào cường độ ban đầu sóng xung kích 45 3.2.2.2 Sự phụ thuộc vào bán kính bọt 47 3.2.2.3 Sự phụ thuộc vào thể tích pha 49 KẾT QUẢ VÀ BÀN LUẬN 52 Kết luận 52 Kiến nghị nghiên cứu 52 TÀI LIỆU THAM KHẢO 54 Số hóa Trung tâm Học liệu – Đại học Thái Nguyên http://www.lrc-tnu.edu.vn/ iv DANH MỤC CÁC KÝ HIỆU SỬ DỤNG Chỉ số - Trạng thái cân sau sóng tới sau sóng phản xạ hỗn hợp (i= 1, 2) Chỉ số - Trạng thái ban đầu hỗn hợp, pha lỏng, pha khí, khí hồ tan khơng tan (i = 0,12,v,g) B - Hằng số khí c, cp2, cv2 - Nhiệt dung riêng, nhiệt dungriêng áp suất vận tốc không đổi D(1) - Vận tốc sóng j - Cường độ chuyển pha l - Nhiệt độ hoá nước n - Số lượng bọt p - áp suất hỗn hợp pe - Cường độ sóng xung kích ban đầu q - Dịng nhiệt R - Bán kính bọt T - Nhiệt độ hỗn hợp t - Thời gian v, v(i) - Vận tốc hỗn hợp w1 - Vận tốc hướng kính chất lỏng chuyển động xung quanh bọt w1 - Vận tốc màng bọt Xi - Phần khối lượng pha i  - Phần thể tích pha hỗn hợp  - Hệ số sức căng bề mặt  - Toạ độ Lagrange 1, 2 - Hệ số điều chỉnh hỗn hợp 1 - Hệ số nhớt động lực chất lỏng Số hóa Trung tâm Học liệu – Đại học Thái Nguyên http://www.lrc-tnu.edu.vn/ v DANH MỤC CÁC HÌNH VẼ Trang Hình 3.1: Quá trình truyền nhiệt xung quanh bọt từ pha sang pha lỏng áp suất lan truyền hỗn hợp Pmax  1.5 24 Hình 3.2: Quá trình truyền nhiệt xung quanh bọt từ pha sang pha lỏng áp suất lan truyền hỗn hợp P max  24 Hình 3.3: Quá trình truyền nhiệt xung quanh bọt từ pha sang pha lỏng áp suất lan truyền hỗn hợp Pmax  2.5 25 Hình 3.4: Sự giảm nhiệt độ dọc theo ống xung kích xung áp suất lan truyền hỗn hợp nước .25 Hình 3.5: Quá trình truyền nhiệt xung quanh bọt từ pha sang pha sang pha lỏng xung áp suất lan truyền hỗn hợp R0  0.001m 26 Hình 3.6: Quá trình truyền nhiệt xung quanh bọt từ pha sang pha sang pha lỏng xung áp suất lan truyền hỗn hợp R0  0.0012m 27 Hình 3.7:Quá trình truyền nhiệt xung quanh bọt từ pha sang pha sang pha lỏng xung áp suất lan truyền hỗn hợp R0  0.0015m 27 Hình 3.8: Quá trình truyền nhiệt xung quanh bọt từ pha sang pha lỏng xung áp suất lan truyền hỗn hợp  20  2% .28 Hình 3.9: trình truyền nhiệt xung quanh bọt từ pha sang pha lỏng xung áp suất lan truyền hỗn hợp  20  4% 29 Hình 3.10: Quá trình truyền nhiệt xung quanh bọt từ pha sang pha lỏng xung áp suất lan truyền hỗn hợp  20  6% 29 Hình 3.11: Sự giảm nhiệt độ dọc theo chiều dài ống xung kích 30 Hình 3.12: Quá trình truyền nhiệt từ pha sang pha lỏng xung áp suất lan truyền hỗn hợp Pmax  1.5 .31 Hình 3.13: Quá trình truyền nhiệt từ pha sang pha lỏng xung áp suất lan truyền hỗn hợp Pmax  .31 Hình 3.14: Quá trình truyền nhiệt từ pha sang pha lỏng xung áp suất lan truyền hỗn hợp Pmax  2.5 32 Hình 3.15: Sự giảm nhiệt độ dọc theo ống xung kích xung áp suất lan truyền hỗn hợp nước .33 Số hóa Trung tâm Học liệu – Đại học Thái Nguyên http://www.lrc-tnu.edu.vn/ vi Hình 3.16: Quá trình truyền nhiệt từ pha sang pha lỏng xung áp suất lan truyền hỗn hợp ni tơ R0  0.001m 33 Hình3.17: Quá trình truyền nhiệt từ pha sang pha lỏng xung áp suất lan truyền hỗn hợp ni tơ R0  0.0012m 34 Hình 3.18: Quá trình truyền nhiệt từ pha sang pha lỏng xung áp suất lan truyền hỗn hợp ni tơ R0  0.0015m 34 Hình 3.19: Quá trình truyền nhiệt xung quanh bọt từ pha sang pha lỏng xung áp suất lan truyền hỗn hợp  20  2% .35 Hình 3.20: Quá trình truyền nhiệt xung quanh bọt từ pha sang pha lỏng xung áp suất lan truyền hỗn hợp  20  4% .35 Hình 3.21: Quá trình truyền nhiệt xung quanh bọt từ pha sang pha lỏng xung áp suất lan truyền hỗn hợp  20  6% 36 Hình 3.22: Sự giảm nhiệt độ dọc theo ống xung kích xung áp suất lan truyền hỗn hợp nước nitơ 37 Hình 3.23: Sự giảm nhiệt độ xung quanh bọt xung áp suất lan truyền hỗn hợp nước nito, tương ứng với đường cong (1), (2) .37 Hình 3.24: Quá trình truyền nhiệt xung quanh bọt từ pha sang pha lỏng cường độ sóng dài P max  1.4 39 Hình 3.25: Quá trình truyền nhiệt xung quanh bọt từ pha sang pha lỏng cường độ sóng dài Pmax  1.5 40 Hình 3.26: Quá trình truyền nhiệt xung quanh bọt từ pha sang pha lỏng cường độ sóng dài Pmax  1.6 .40 Hình 3.27: Sự giảm nhiệt độ dọc theo ống xung kích hỗn hợp nước 41 Hình 3.28: Quá trình truyền nhiệt xung quanh bọt từ pha sang pha lỏng bán kính bọt R0  0.001m 42 Hình 3.29: Quá trình truyền nhiệt xung quanh bọt từ pha sang pha lỏng bán kính bọt R0  0.0012m 42 Hình 3.30: Quá trình truyền nhiệt xung quanh bọt từ pha sang pha lỏng bán kính bọt R0  0.0015m 42 Hình 3.31: Quá trình truyền nhiệt xung quanh bọt từ pha sang pha lỏng thể tích pha  20  2% 43 Số hóa Trung tâm Học liệu – Đại học Thái Nguyên http://www.lrc-tnu.edu.vn/ vii Hình 3.32: Quá trình truyền nhiệt xung quanh bọt từ pha sang pha lỏng thể tích pha  20  4% 44 Hình 3.33: Quá trình truyền nhiệt xung quanh bọt từ pha sang pha lỏng thể tích pha  20  5% 44 Hình 3.34: Quá trình truyền nhiệt xung quanh bọt từ pha sang pha lỏng cường độ sóng dài Pmax  1.4 .45 Hình 3.35: Quá trình truyền nhiệt xung quanh bọt từ pha sang pha lỏng cường độ sóng dài P max  1.5 46 Hình 3.36: Quá trình truyền nhiệt xung quanh bọt từ pha sang pha lỏng cường độ sóng dài P max  1.6 46 Hình 3.37: Sự giảm nhiệt độ dọc theo ống xung kích hỗn hợp ni tơ 47 Hình 3.38: Quá trình truyền nhiệt xung quanh bọt từ pha sang pha lỏng bán kính bọt R0  0.001m 47 Hình 3.39: Quá trình truyền nhiệt xung quanh bọt từ pha sang pha lỏng bán kính bọt R0  0.0012m 48 Hình 3.40: Quá trình truyền nhiệt xung quanh bọt từ pha sang pha lỏng bán kính bọt R0  0.0015m 48 Hình 3.41: Quá trình truyền nhiệt xung quanh bọt từ pha sang pha lỏng thể tích pha  20  2% 49 Hình 3.42: Quá trình truyền nhiệt xung quanh bọt từ pha sang pha lỏng thể tích pha  20  4% 50 Hình 3.43: Quá trình truyền nhiệt xung quanh bọt từ pha sang pha lỏng thể tích pha  20  5% 50 Hình 3.44: Sự giảm nhiệt độ dọc theo ống xung kích xung áp suất lan truyền hỗn hợp nước ni tơ 51 Số hóa Trung tâm Học liệu – Đại học Thái Nguyên http://www.lrc-tnu.edu.vn/ MỞ ĐẦU Hỗn hợp chất lỏng chứa bọt - môi trường hai pha phức tạp Môi trường thường gặp nhiều trình tự nhiên, cơng nghệ hóa học, vật lý ứng dụng nhiều ngành công nghiệp lượng như: khai thác, vận chuyển chế biến dầu khí Tuy nhiên, từ năm 1950 trở lại đây, việc nghiên cứu dòng hai pha lỏng - khí tiến hành cách có hệ thống lý thuyết thực nghiệm Trong hỗn hợp dạng trình trao đổi nhiệt - khối lượng tượng quan trọng tách rời trường hợp tồn sóng xung kích lan truyền hỗn hợp xảy tượng tương tác pha hỗn hợp Mơi trường hỗn hợp hai pha có tính chất đặc biệt chỗ hỗn hợp có kết hợp tính chất phi tuyến vật lý mạnh, tán sắc hao tán lượng nên biểu đồ mơ tả sóng có nhiều dạng Chính vậy, thay đổi điều kiện thuỷ động lực dẫn đến thay đổi cấu trúc sóng, tính chất vật lý nhiệt, q trình tương tác pha Tính chất đặc trưng hỗn hợp chất lỏng chứa bọt khí trình động lực học xuất biến dạng cục hỗn hợp thay đổi thể tích mơi trường thay đổi thể tích bọt Khả mức độ co nén bọt phụ thuộc mạnh vào trao đổi nhiệt khối lượng pha lỏng pha khí Sự xuất đồng thời lượng biến dạng dẫn tới sóng có cấu trúc khác Ngồi truyền sóng áp suất mơi trường dẫn đến khả hóa ngưng tụ pha khí, từ dẫn đến thay đổi chủ yếu cấu trúc môi trường Do hỗn hợp chất lỏng chứa bọt khí có tính chất đặc biệt hỗn hợp xuất nhiều lĩnh vực công nghiệp lượng, cơng nghệ hóa học q trình tự nhiên cho nên, hiểu biết tượng xuất sóng áp suất lan truyền qua chất lỏng chứa bọt xảy q trình tương tác sóng cần thiết Số hóa Trung tâm Học liệu – Đại học Thái Nguyên http://www.lrc-tnu.edu.vn/ 41 Pmax T T0 1.4 1.025551 1.5 1.030692 1.6 1.037882 Có thể thấy sóng lan truyền hỗn hợp nhiệt độ hỗn hợp giảm dọc theo chiều dài ống.Ví dụ xét vị trí cách đầu ống xung kích 0.02, 0.2, 0.4, 0.6m màng bọt ta có kết sau: 1.04 (a) (b) 1.035 T/T0 1.03 1.025 (c) 1.02 1.015 1.01 1.005 0.1 0.2 0.3 0.4 0.5 0.6 l (m) Hình 3.27: Sự giảm nhiệt độ dọc theo ống xung kích hỗn hợp nước Pmax Ở đường cong (a), (b), (c) tương ứng với Pmax  1.6 , Pmax  1.5  Đồng thời thấy trường nhiệt độ nhận lớp chất lỏng xung quanh bọt giảm dần Vậy với điều kiện bán kính bọt thể tích pha thay đổi cường độ ban đầu sóng xung kích cường độ sóng xung kích tăng bọt co nén mạnh kéo theo trao đổi nhiệt tăng, giảm nhiệt độ dọc theo ống xung kích mạnh 3.2.1.2 Sự phụ thuộc vào bán kính bọt Các kết nhận hình từ 3.28 đến 3.30 đường cong biểu diễn trình trao đổi nhiệt pha phụ thuộc vào bán kính bọt hỗn hợp nước chứa bọt Ở ta lấy cường độ sóng dài P max  1.4 , thể tích pha  20  4% bán nhận giá trị là: R0  0.001m ; 0.0012m 0.0015m Số hóa Trung tâm Học liệu – Đại học Thái Nguyên http://www.lrc-tnu.edu.vn/ 42 Trên hình 3.28 tương ứng với R0  0.001m 1.03 1.025 T/T0 1.02 1.015 4 1.01 1.005 (a) (c ) (b) 0.002 0.004 0.006 0.008 0.01 0.012 0.014 0.016 t (s) Hình 3.28: Quá trình truyền nhiệt xung quanh bọt từ pha sang pha lỏng bán kính bọt R0  0.001m Trên hình 3.29 tương ứng với R0  0.0012m 1.03 1.025 T/T0 1.02 1.015 (a) (b) 4 1.01 1.005 (c) 0.002 0.004 0.006 0.008 0.01 0.012 0.014 0.016 t (s) Hình 3.29: Quá trình truyền nhiệt xung quanh bọt từ pha sang pha lỏng bán kính bọt R0  0.0012m Trên hình 3.30 tương ứng với R0  0.0015m 1.03 1.025 T/T0 1.02 1.015 1.01 1.005 (a) (b) (c ) 0.002 0.004 0.006 0.008 0.01 0.012 0.014 t (s) Hình 3.30: Quá trình truyền nhiệt xung quanh bọt từ pha sang pha lỏng bán kính bọt R0  0.0015m Số hóa Trung tâm Học liệu – Đại học Thái Nguyên http://www.lrc-tnu.edu.vn/ 43 Từ kết thu ta có nhận xét sau: Với điều kiện thể tích pha có cường độ sóng dài tác động vào hỗn hợp thay đổi bán kính bọt có ảnh hưởng đến q trình trao đổi nhiệt pha, bán kính bọt tăng nhiệt độ thu màng bọt chất lỏng tăng lên Cụ thể xét vị trí cách đầu ống xung kích khoảng   0.02m ta thu nhiệt độ T sau: T0 T T0 R0 0.001m 0.0012m 0.0015m 1.025374 1.025438 1.025567 Qua kết ta thấy, với điều kiện thể tích pha cường độ sóng xung kích thay đổi bán kính bọt ảnh hưởng bán kính bọt lên q trình trao đổi nhiệt pha không nhiều 3.2.1.3 Sự phụ thuộc vào thể tích pha So sánh kết nhận trình trao đổi nhiệt xung quanh bọt từ pha sang pha lỏng điều kiện đầu điều kiện biên không đổi với P max  1.4 R0  0.001m thay đổi thể tích pha  20  2% ,  20  4%  20  5% ta thu kết sau đây: Hình 3.31 ứng với  20  2% 1.03 1.025 T/T0 1.02 1.015 1 4 1.01 1.005 (a) (b) (d) (c ) 0.002 0.004 0.006 0.008 0.01 0.012 t (s) Hình 3.31: Quá trình truyền nhiệt xung quanh bọt từ pha sang pha lỏng thể tích pha  20  2% Số hóa Trung tâm Học liệu – Đại học Thái Nguyên http://www.lrc-tnu.edu.vn/ 44 Hình 3.32 ứng với  20  4% 1.03 1.025 T/T0 1.02 1.015 4 1.01 1.005 (a) (c ) (b) 0.002 0.004 0.006 0.008 0.01 0.012 0.014 0.016 t (s) Hình 3.32: Quá trình truyền nhiệt xung quanh bọt từ pha sang pha lỏng thể tích pha  20  4% Hình 3.33 ứng với  20  5% T/T0 1.03 1.025 1.02 1.015 1.01 1.005 (b) (a) 0.002 0.004 0.006 (c ) 0.008 0.01 0.012 0.014 t (s) Hình 3.33: Quá trình truyền nhiệt xung quanh bọt từ pha sang pha lỏng thể tích pha  20  5% Qua ta thấy thể tích pha tăng nhiệt độ trao đổi nhiệt pha gần khơng thay đổi Cụ thể xét vị trí cách đầu ống xung kích khoảng   0.02m cách tâm bọt khoảng r *  ta có kết sau:  20 2% 4% 5% T T0 1.025551 1.025574 1.025581 Vậy qua nhận xét ta thấy trình trao đổi nhiệt từ pha sang pha lỏng hỗn hợp nước chứa bọt phụ thuộc nhiều vào cường độ Số hóa Trung tâm Học liệu – Đại học Thái Nguyên http://www.lrc-tnu.edu.vn/ 45 ban đầu sóng xung kích tác động lên hỗn hợp phụ thuộc vào thay đổi bán kính bọt thể tích pha 3.2.2 Quá trình trao đổi nhiệt pha hỗn hợp ni tơ chứa bọt Môi trường ni tơ lỏng chứa bọt sử dụng nhiệt độ T  77,350 K điều kiện áp suất P0  0.1MPa , hỗn hợp chứa ống khảo sát nằm ngang Trong phần ta xét ảnh hưởng cường độ sóng dài, bán kính bọt thể tích pha đến q trình trao đổi nhiệt pha hỗn hợp ni tơ chứa bọt 3.2.2.1 Sự phụ thuộc vào cƣờng độ ban đầu sóng xung kích Trong hình từ 3.34 đến 3.36 đường cong biểu diễn trình trao đổi nhiệt pha phụ thuộc vào cường độ ban đầu sóng xung kích Trên hình kết nhận lấy cường độ sóng dài Pmax  1.4 , Pmax  1.5 Pmax  1.6 có bán kính bọt R0  0.001m có thể tích pha  20  2% Trên hình đường cong từ đến thuộc nhóm đường cong (a), (b), (c), (d) kết nhiệt độ nhận cách tâm bọt r *  1;5;10;15 (trong r *  r , r khoảng cách từ tâm bọt tới vị trí xác định nhiệt) R0 Kết hình 3.34 tương ứng với Pmax  1.4 1.05 T/T0 1.04 1.03 1.02 1.01 (a ) 0.001 (b) 0.002 0.003 (c ) 0.004 0.005 (d) 0.006 0.007 0.008 0.009 0.01 t (s) Hình 3.34: Quá trình truyền nhiệt xung quanh bọt từ pha sang pha lỏng cường độ sóng dài Pmax  1.4 Số hóa Trung tâm Học liệu – Đại học Thái Nguyên http://www.lrc-tnu.edu.vn/ 46 Kết hình 3.35 tương ứng với Pmax  1.5 1.06 1.05 T/T0 1.04 1.03 1.02 1.01 (a) 0.001 0.002 0.003 (d) (c) (b) 0.004 0.005 0.006 0.007 0.008 0.009 0.01 t (s) Hình 3.35: Quá trình truyền nhiệt xung quanh bọt từ pha sang pha lỏng cường độ sóng dài P max  1.5 Kết hình 3.36 tương ứng với Pmax  1.6 1.07 1.06 T/T0 1.05 1.04 1.03 1.02 1.01 (a) 0.001 (b) 0.002 0.003 (c) 0.004 (d) 0.005 0.006 0.007 0.008 0.009 0.01 t (s) Hình 3.36: Quá trình truyền nhiệt xung quanh bọt từ pha sang pha lỏng cường độ sóng dài P max  1.6 Qua kết biểu diễn cho thấy cường độ ban đầu P max nhận giá trị tăng dần từ 1.4, 1.5, đến 1.6 nhiệt độ thu màng bọt tăng lên Cụ thể xét vị trí cách đầu ống xung kích khoảng 0.02m ta nhận kết sau: P max T T0 1.4 1.038386 1.5 1.047951 1.6 1.057635 Số hóa Trung tâm Học liệu – Đại học Thái Nguyên http://www.lrc-tnu.edu.vn/ 47 Đồng thời có giảm nhiệt độ dọc theo chiều dài ống xung kích trường nhiệt độ xung quanh bọt giảm ví dụ xét vị trí cách đầu ống xung kích 0.02, 0.2, 0.4 0.6m màng bọt ta có kết sau 1.07 (a) 1.06 (b) (c) T/T0 1.05 1.04 1.03 1.02 1.01 0.1 0.2 0.3 0.4 0.5 0.6 l (m) Hình 3.37: Sự giảm nhiệt độ dọc theo ống xung kích hỗn hợp ni tơ Ở đường cong (a), (b), (c) tương ứng với cường độ sóng tác động 1.6; 1.5; 1.4 Vậy với điều kiện bán kính bọt thể tích pha khác cường độ sóng dài tác động vào hỗn hợp cường độ sóng xung kích tăng kéo theo bọt co nén mạnh dẫn đến trao đổi nhiệt pha tăng giảm nhiệt độ dọc theo chiều dài ống xung kích mạnh 3.2.2.2 Sự phụ thuộc vào bán kính bọt Các kết nhận hình từ 3.37 đến 3.40 đường cong biểu diễn trình trao đổi nhiệt pha phụ thuộc vào bán kính bọt hỗn hợp ni tơ chứa bọt Ở ta lấy cường độ sóng dài P max  1.4 , thể tích pha  20  2% bán kính bọt nhận giá trị là: R0  0.001m ; 0.0012m 0.0015m Trên hình 3.38 tương ứng với R0  0.001m 1.05 1.04 T/T0 1.03 1.02 1.01 (a ) 0.001 (b) 0.002 0.003 (c ) 0.004 0.005 (d) 0.006 0.007 0.008 0.009 0.01 t (s) Hình 3.38: Quá trình truyền nhiệt xung quanh bọt từ pha sang pha lỏng bán kính bọt R0  0.001m Số hóa Trung tâm Học liệu – Đại học Thái Nguyên http://www.lrc-tnu.edu.vn/ 48 Trên hình 3.39 tương ứng với R0  0.0012m 1.04 1.03 T/T0 1.02 1.01 (a) (b) 0.002 0.004 (c) 0.006 (d) 0.008 0.01 0.012 0.014 t (s) Hình 3.39: Quá trình truyền nhiệt xung quanh bọt từ pha sang pha lỏng bán kính bọt R0  0.0012m Trên hình 3.49 tương ứng với R0  0.0015m 1.04 T/T0 1.03 1.02 1.01 (a ) (b) (c) (d) 0.002 0.004 0.006 0.008 0.01 0.012 0.014 t (s) ) Hình 3.40: Quá trình truyền nhiệt xung quanh bọt từ pha sang pha lỏng bán kính bọt R0  0.0015m Từ kết thu ta có nhận xét sau: Với điều kiện thể tích pha có cường độ sóng dài tác động vào hỗn hợp thay đổi bán kính bọt có ảnh hưởng đến q trình trao đổi nhiệt pha, bán kính bọt tăng nhiệt độ thu màng bọt gần không thay đổi Cụ thể xét vị trí cách đầu ống xung kích khoảng   0.02m ta thu nhiệt độ T sau: T0 Số hóa Trung tâm Học liệu – Đại học Thái Nguyên http://www.lrc-tnu.edu.vn/ 49 T T0 R0 0.001m 1.038386 0.0012m 1.036828 0.0015m 1.035113 Qua kết ta thấy, với điều kiện thể tích pha cường độ sóng xung kích thay đổi bán kính bọt thay đổi bán kính bọt có ảnh hưởng đến q trình trao đổi nhiệt pha hỗn hợp ni tơ, nhiên thay đổi nhiệt không nhiều 3.2.2.3 Sự phụ thuộc vào thể tích pha So sánh kết nhận trình trao đổi nhiệt xung quanh bọt từ pha sang pha lỏng hỗn hợp ni tơ lỏng chứa bọt điều kiện đầu điều kiện biên không đổi với P max  1.4 R0  0.001m thay đổi thể tích pha hơi,  20  2% ,  20  4%  20  5% ta thu kết sau đây: Hình 3.41 ứng với  20  2% 1.05 T/T0 1.04 1.03 1.02 1.01 (a ) 0.001 (b) 0.002 0.003 (c ) 0.004 0.005 (d) 0.006 0.007 0.008 0.009 0.01 t (s) Hình 3.41: Quá trình truyền nhiệt xung quanh bọt từ pha sang pha lỏng thể tích pha  20  2% Số hóa Trung tâm Học liệu – Đại học Thái Nguyên http://www.lrc-tnu.edu.vn/ 50 Hình 3.42 ứng với  20  4% 1.04 T/T0 1.03 1.02 1.01 (a) (b) 0.002 (c) 0.004 (d) 0.006 0.008 0.01 0.012 0.014 t (s) Hình 3.42: Quá trình truyền nhiệt xung quanh bọt từ pha sang pha lỏng thể tích pha  20  4% Hình 3.43 ứng với  20  5% 1.04 1.03 T/T0 1.02 1.01 (a) 0.002 (b) (d) (c) 0.004 0.006 0.008 0.01 0.012 0.014 0.016 t (s) Hình 3.43: Quá trình truyền nhiệt xung quanh bọt từ pha sang pha lỏng thể tích pha  20  5% Từ kết ta thấy, thay đổi thể tích pha hỗn hợp Ni tơ trình trao đổi nhiệt pha khơng thay đổi nhiều Ví dụ xét vị trí cách đầu ống xung kích khoảng   0.02m cách tâm bọt khoảng r *  màng bọt ta có kết sau:  20 2% 4% 5% T T0 1.038386 1.037757 1.037562 Số hóa Trung tâm Học liệu – Đại học Thái Nguyên http://www.lrc-tnu.edu.vn/ 51 Trên số kết nghiên cứu trình trao đổi nhiệt pha có sóng xung kích lan truyền hỗn hợp ni tơ chứa bọt Khi xét đến ảnh hưởng tham số đặc trưng điều kiện đầu, điều kiện biên cấu trúc hỗn hợp tham số vật lý nhiệt đến q trình trao đổi nhiệt pha q trình truyền sóng ta thấy trình truyền nhiệt hỗn hợp ni tơ phụ thuộc nhiều vào cường độ ban đầu sóng dài tác động lên hỗn hợp, phụ thuộc vào thể tích bán kính bọt Sau nghiên cứu trình trao đổi nhiệt pha hỗn hợp nước ni tơ ta có so sánh sau: Trên hình 3.44 giảm nhiệt độ dọc theo chiều dài ống xung kích hỗn hợp nước ni tơ, ta lấy điều kiện đầu điều kiện biên sau: Cường độ ban đầu P max  1.4 , bán kính bọt R0  0.001m thể tích pha  20  2% hình vẽ đường cong (a) ứng với kết ni tơ đường cong (b) ứng với kết nước 1.05 (a) T/T0 1.04 (b) 1.03 1.02 1.01 0.1 0.2 0.3 0.4 0.5 0.6 l (m) Hình 3.44: Sự giảm nhiệt độ dọc theo ống xung kích xung áp suất lan truyền hỗn hợp nước ni tơ Qua đó, ta thấy giảm nhiệt độ nước mạnh ni tơ Từ kết nhận thấy, điều kiện đầu, điều kiện biên, sử dụng áp suất ban đầu tác động vào hỗn hợp lỏng truyền nhiệt từ pha sang pha lỏng hỗn hợp khác khác Từ kết thấy q trình trao đổi nhiệt pha hỗn hợp nước mạnh so với hỗn hợp ni tơ Số hóa Trung tâm Học liệu – Đại học Thái Nguyên http://www.lrc-tnu.edu.vn/ 52 KẾT QUẢ VÀ BÀN LUẬN Kết luận Trên sở kết nghiên cứu trình trao đổi nhiệt pha hỗn hợp chất lỏng chứa bọt có sóng xung kích lan truyền, đưa số kết luận sau: - Luôn xảy tượng trao đổi nhiệt pha hỗn hợp - Trường nhiệt độ lớp lỏng xung quanh bọt giảm khoảng cách từ lớp chất lỏng tới tâm bọt tăng - Quá trình trao đổi nhiệt pha hỗn hợp phụ thuộc tính chất vật lý nhiệt hỗn hợp Tuy nhiên, trình trao đổi nhiệt pha giảm dần theo chiều dài ống xung kích - Q trình trao đổi nhiệt khối lượng pha hỗn hợp phụ thuộc mạnh vào cường độ ban đầu sóng xung kích vào nồng độ thể tích pha cịn phụ thuộc yếu vào bán kính bọt Từ kết phân tích cho thấy q trình trao đổi nhiệt pha hỗn hợp nước mạnh so với hỗn hợp nitơ lỏng Kiến nghị nghiên cứu Trong nội dung luận văn, dựa sở hệ phương trình thuỷ nhiệt động lực học với việc sử dụng thuật toán chương trình tính phù hợp nghiên cứu q trình trao đổi nhiệt pha số hỗn hợp chất lỏng chứa bọt như: hỗn hợp nước chứa bọt hơivà nitơ lỏng chứa bọt Các kết cho thấy ảnh hưởng cường độ ban đầu sóng, nồng độ thể tích pha tính chất vật lý nhiệt hỗn hợp lên trình trao đổi nhiệt pha cho thấy tượng tách rời trình động lực học sóng Tuy nhiên, tranh tồn cảnh q trình sóng mơi trường chất lỏng chứa bọt nhiềt vấn đề phức tạp mà nhà khoa học giới nước quan tâm ví dụ như: Phân tích chế độ làm việc độ trạm lượng, phân tích Số hóa Trung tâm Học liệu – Đại học Thái Nguyên http://www.lrc-tnu.edu.vn/ 53 điều kiện hư hỏng đảm bảo an toàn khai thác nhà máy điện nguyên tử, sử dụng buồng bọt khảo sát hạt bản, phân tích tượng xâm thực máy tuôcbin, sử dụng để khai thác dầu mỏ, vận chuyển chế biến dầu khí… Đó vấn đề nghiên cứu theo hướng luận văn Số hóa Trung tâm Học liệu – Đại học Thái Nguyên http://www.lrc-tnu.edu.vn/ 54 TÀI LIỆU THAM KHẢO Tiếng Việt Dương Ngọc Hải (1989), “ảnh hưởng trao đổi nhiệt - chất lên động lực học sóng mơi trường hai pha”, Tạp chí Cơ học, 11(3), pp 28-32 Dương Ngọc Hải (2001), “Năng lượng nguyên tử - an toàn hạt nhân số vấn đề nhiệt thuỷ học”, Tuyển tập hội nghị khoa học toàn quốc Cơ học kỹ thuật, HàNội, pp 68 - 74 Nguyễn Văn Tuấn, Dương Ngọc Hải (2006) “ Quá trình lan truyền tương tác sóng ngắn số hỗn hợp lỏng-hơi” , Tuyển tập cơng trình Hội nghị Cơ học Thủy khí Tồn quốc, Vũng tàu , 26-28/7/2006 pp.515-526 Dương Ngọc Hải, Nguyễn Văn Tuấn, (2007), “ Sự truyền nhiệt pha trình lan truyền sóng xung kích số hỗn hợp lỏng-hơi” , Tuyển tập cơng trình khoa học Hội nghị Cơ học Thủy khí Tồn quốc, Huế , 2628/7/2007, pp 515-526 Nguyễn Thị Thu Hằng (2010), “Nghiên cứu trao đổi nhiệt khối lượng pha q trình truyền sóng số hỗn hợp chất lỏng chứa bọt.” Luận văn thạc sĩ Cơ học kỹ thuật , trường đại học kỹ thuật công nghiệp 11/2010 Trần Quốc Nghị (2011), “Nghiên cứu q trình lan truyền sóng áp suất hỗn hợp chất lỏng hai pha” Luận văn thạc sĩ Toán học ứng dụng, trường đại học khoa học 11/2011 Tiếng Anh Duong Ngoc Hai and Nguyen Van Tuan (2002), “Waves reflected by solid wall in the mixture of liquid with vapour bubbles”, J Mechanics Vol 24, No 3, pp 167-180 Duong Ngọc Hai and Nguyen Van Tuan, (2008), “Interphase Heat-Mass Transfer in Bubbly Liquid and its Influence on Wave Propagation Processes”, Proceedings of the 2nd International Forum on Heat Transfer, September 17-19, 2008, Tokyo, Japan, paper ID N128 Elperin T and Fominykh A (1995), “Heat and Mass Transfer During Gas Số hóa Trung tâm Học liệu – Đại học Thái Nguyên http://www.lrc-tnu.edu.vn/ 55 Hydrate Formation in Gas-Liquid Slug Flow”, Pro of the 2nd Int Conf on Multiphase Flow’95-Kyoto, Kyoto, Japan, pp IF2_15-IF2_22 10 Fujii T and Akagawa k (1995), “Water hammers phenomena in one-component two-phase bubbly flow”, Proc of the Int Conf on Multiphase Flow’ 95, Kyoto, Japan, pp - 22 11 Gasenko V.G (1977), “The Structure of Stasionary Shock Waves in Gas-Liquid Mixture with Thermal Relaxation”, in: Thermophysical Investigation, Institute of Thermophysics SD Academy of Sciences of the USSR, Novosibirsk, pp.42-46 12 Gubaidulin A.A., Beregova O.Sh and Bekishev S.A (1998), “Features of NonStationary Shock Waves Propagation in Non-Newtonian Liquid with Gas Bubbles”, Pro of Int Conf on Multiphase Flow, ICMF’98, Lyon, France, from CD - ICMF, Sesion 3.2, Unit 249 13 Kubie J (2000), “Velocity of lony Bubbles in Horizontally Oscillating Vertical Piper”, Int J Multiphase Flow Vol 26, pp 339-349 14 Kwidzinki R., Karda D and Pribaturin N.A (1998), “Experimental investigation of structure of stationary shock wave and its interaction with transient impulse of pressure in two-phase flow”, Proc of Int Conf on Multiphase Flow ICMF’ 98, Lyon, France, from CD - ICMF, Sesion 3.2, Unit 353 15 Nakoryakov V.E., Pokusaev B.G., Shreiber I.R., Kuznetsov V.V., Malykh N.V (1975), “Experimental Investigation of Shock Waves in a Fluid with Gas Bubbles” in: Wave Processes in Two-Phase Systems, Institute of Thermophysics SD Academy of Sciences of the USSR, Novosibirsk, pp 54-97 16 Gubaidulin A.A (1995), “Shock Wave Phenomena in Bubbles Liquids”, Pro of the 2nd Int Conf on Multiphase Flow’95-Kyoto, pp PH2_3-PH2_8, 17 Nakoryako V.E., Pokusaev B.G and shreiber I.R (1983), Wave Propagation Processes in Gas and Vapour-Liquid Media, Institute of Thermophysics, Novosibirsk, pp 109-133 18 Nigmatulin R.I (1990), Dynamics of Multiphase Media Hemisphere, publ Corp, Washington Số hóa Trung tâm Học liệu – Đại học Thái Nguyên http://www.lrc-tnu.edu.vn/