Header Page of 126 ĐẠI HỌC THÁI NGUYÊN TRƯỜNG ĐẠI HỌC KHOA HỌC NGUYỄN HUY VIỆT NGHIÊN CỨU HIỆN TƯỢNG SÓNG VA TRONG CHẤT LỎNG CHỨA BỌT HƠI LUẬN VĂN THẠC SỸ TOÁN HỌC THÁI NGUYÊN - 2013 Footer Page of 126 Header Page of 126 ĐẠI HỌC THÁI NGUYÊN TRƯỜNG ĐẠI HỌC KHOA HỌC NGUYỄN HUY VIỆT NGHIÊN CỨU HIỆN TƯỢNG SÓNG VA TRONG CHẤT LỎNG CHỨA BỌT HƠI Chuyên ngành: Toán ứng dụng Mã số: 60 46 01 12 LUẬN VĂN THẠC SỸ TOÁN HỌC Người hướng dẫn khoa học: PGS.TS NGUYỄN VĂN TUẤN THÁI NGUYÊN - 2013 Footer Page of 126 Header Page of 126 MỤC LỤC LỜI CẢM ƠN Danh sách kí hiệu chữ viết tắt MỞ ĐẦU CHƢƠNG I TỔNG QUAN CHƢƠNG II - CƠ SỞ LÝ THUYẾT 12 2.1 HIỆN TƢỢNG NƢỚC VA 12 2.1.1 Nƣớc va đóng khoá tức thời …………………………………………….14 2.1.2 Nƣớc va đóng khoá từ từ……………………………………………… 17 2.2 HIỆN TƢỢNG SÓNG VA TRONG HỖN HỢP CHẤT LỎNG CHỨA BỌT HƠI 19 2.2.1 Hệ phƣơng trình thủy nhiệt động lực học hỗn hợp chất lỏng chứa bọt hơi……………………………………………………………………………… 20 2.2.2 Phƣơng pháp giải số………………………………………………………….25 2.2.3 Chƣơng trình tính toán:…………………………………………………… 27 Chƣơng III 30 MỘT SỐ KẾT QUẢ NGHIÊN CỨU HIỆN TƢỢNG SÓNG VA TRONG HỖN HỢP CHẤT LỎNG CHỨA BỌT HƠI 30 KẾT LUẬN 41 KIẾN NGHỊ VỀ NHỮNG NGHIÊN CỨU TIẾP THEO 42 TÀI LIỆU THAM KHẢO 43 Số hóa trung tâm học liệu Footer Page of 126 http://www.lrc.tnu.edu.vn/ Header Page of 126 LỜI CẢM ƠN Để hoàn thành luận văn Thạc sĩ em xin đƣợc gửi lời cảm ơn chân thành tới Ban giám hiệu, phòng Đào tạo, Khoa Toán - Tin, thầy cô giáo trƣờng Đại học Khoa học - Đại học Thái Nguyên thầy giáo công tác Viện Toán học Việt Nam, thầy cô giáo nhiệt tình truyền đạt cho em kiến thức quý báu suốt trình học tập hoàn thành luận văn Thạc sĩ Đặc biệt, em xin đƣợc gửi lời biết ơn chân thành tới Phó Giáo sƣ Tiến sĩ Nguyễn Văn Tuấn Ngƣời trực tiếp bảo, hƣớng dẫn động viên em suốt trình nghiên cứu hoàn thành luận văn Xin đƣợc gửi lời cảm ơn đến Ban giám hiệu, bạn đồng nghiệp trƣờng Trung Tâm Giáo dục thƣờng xuyên Mƣờng La - Tỉnh Sơn La quan tâm, tạo điều kiện thuận lợi để hoàn thành khóa học Xin cảm ơn anh chị em học viên lớp cao học K5C đoàn kết, đùm bọc giúp đỡ toàn khóa học Cuối xin đƣợc gửi lời biết ơn sâu sắc đến ngƣời thân gia đình tôi, ngƣời động viên, khuyến khích giúp đỡ suốt trình học tập Thành đạt đƣợc quà mà muốn dành tặng cho gia đình thân yêu Thái Nguyên, tháng năm 2013 Học viên Nguyễn Huy Việt Số hóa trung tâm học liệu Footer Page of 126 http://www.lrc.tnu.edu.vn/ Header Page of 126 Danh sách kí hiệu chữ viết tắt B - Hằng số khí Nhiệt dung riêng, nhiệt dung riêng áp suất vận tốc không đổi c, c p , cv j - Cƣờng độ chuyển pha l - Nhiệt hoá chất lỏng n - Số lƣợng bọt p,p0 - Áp suất hỗn hợp, áp suất ban đầu pi - Áp suất pha i pe - Cƣờng độ sóng xung kích ban đầu q - Dòng nhiệt R,a - Bán kính bọt R0, a0 - Bán kính ban đầu bọt T,T0 - Nhiệt độ hỗn hợp, nhiệt độ ban đầu Ti - Nhiệt độ pha i t - Thời gian v, vi - Vận tốc hỗn hợp, vận tốc pha thứ i w1 xung quanh bọt w1a - Vận tốc hƣớng kính chất lỏng chuyển động Vận tốc màng bọt , , P,W Đại lƣợng không thứ nguyên mật độ, nhiệt độ, áp suất vận tốc màng bọt i - Phần thể tích pha i hỗn hợp  i0 - Phần thể tích ban đầu pha i Số hóa trung tâm học liệu Footer Page of 126 http://www.lrc.tnu.edu.vn/ Header Page of 126  - Hệ số sức căng bề mặt  - Tọa độ Lagrange 1 ,  - Hệ số hiệu chỉnh hỗn hợp 1 - Hệ số nhớt động lực chất lỏng i ,  - Mật độ pha i, mật độ hỗn hợp  i0 ,  - Mật độ thực pha i, mật độ ban đầu hỗn hợp  Tỷ số nhiệt dung riêng áp suất không đổi thể tích không đổi b1 , b2 , b3 , b4 - Cƣờng độ xung áp suất Số hóa trung tâm học liệu Footer Page of 126 http://www.lrc.tnu.edu.vn/ Header Page of 126 MỞ ĐẦU Chất lỏng không đồng hay hỗn hợp chất lỏng nhiều pha môi trƣờng 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ý… đƣợc ứ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í, thủy điện,… Tuy nhiên, từ 1950 trở lại đây, việc nghiên cứu dòng hai pha khí lỏng đƣợc tiến hành cách có hệ thống lý thuyết thực nghiệm Trong hỗn hợp trình trao đổi nhiệt chất tƣợng quan trọng tách rời, tƣợng sóng va chất lỏng chứa bọt So với môi trƣờng chất lỏng đồng pha môi trƣờng hỗn hợp hai pha khác xa tính chất vật lý, thể 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 trình 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 thủy động lực dẫn đến thay đổi cấu trúc sóng, tính chất vật lý nhiệt, 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 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 Ngoài ra, truyền sóng áp suất môi trƣờng nhƣ 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 có tính chất đặc biệt nhƣ trên, hỗn hợp xuất nhiều lĩnh vực công nghiệp Số hóa trung tâm học liệu Footer Page of 126 http://www.lrc.tnu.edu.vn/ Header Page of 126 lƣợng, công nghệ hóa học trình tự nhiê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 trình tƣơng tác sóng cần thiết Căn vào tình hình phát triển nghiên cứu tƣợng sóng va chất lỏng chứa bọt nƣớc giới, mục tiêu đề tài bao gồm vấn đề sau: Nghiên cứu trình tăng áp suất van, van đóng chặn dòng chảy hai pha chuyển động ống với vận tốc v0 Để thực đƣợc mục tiêu nghiên cứu trên, thực đề tài có tên: Nghiên cứu tƣợng sóng va chất lỏng chứa bọt Nội dung bao gồm phần sau: Chƣơng I: TỔNG QUAN Trình bày tổng quan phát triển xu hƣớng phát triển lĩnh vực sóng va hỗn hợp lỏng - hơi, tƣợng xảy sóng lan truyền hỗn hợp Chƣơng II: CƠ SỞ LÝ THUYẾT Trong phần trình bày tƣợng nƣớc va chất lỏng pha, sở trình bày mô hình vật lý cho tƣợng dòng chảy hai pha xây dựng hệ phƣơng trình thủy nhiệt động lực học để mô tả tƣợng Sử dụng chƣơng trình tính điều kiện biên phù hợp để giải hệ phƣơng trình Chƣơng III: Một số kết nghiên cứu tƣợng sóng va hỗn hợp chất lỏng chứa bọt Số hóa trung tâm học liệu Footer Page of 126 http://www.lrc.tnu.edu.vn/ Header Page of 126 Trong phần trình bày số kết nghiên cứu trình tăng áp suất van đóng, tƣợng sóng va hỗn hợp chất lỏng chứa bọt Nghiên cứu ảnh hƣởng vận tốc hỗn hợp, thời gian đóng van cấu trúc hỗn hợp lên trình tăng áp suất Bản luận văn đƣợc thực trƣờng Đại học Khoa học - Đại học Thái Nguyên hoàn thành từ ý tƣởng nghiên cứu, hƣớng dẫn tận tình giúp đỡ mặt khoa học Phó Giáo sƣ - Tiến sĩ Nguyễn Văn Tuấn Tôi xin đƣợc bày tỏ lòng biết ơn chân thành tới thầy Số hóa trung tâm học liệu Footer Page of 126 http://www.lrc.tnu.edu.vn/ Header Page 10 of 126 CHƢƠNG I TỔNG QUAN Nhƣ nêu phần mở đầu, hỗn hợp hai pha chất lỏng chứa bọt khí có tính chất đặc biệt, kết hợp tính chất phi tuyến vật lý mạnh, phân tán hao tán lƣợng, nên biểu đồ mô tả sóng có nhiều dạng Chính thay đổi điều kiện thủy động lực, dẫn đến thay đổi cấu trúc sóng trình tuơng tác pha Cho nên hiểu biết tƣợng bất thƣờng xuất sóng áp suất lan truyền qua chất lỏng chứa bọt cần thiết để giải loạt toán thực tế nhƣ để phân tích chế độ làm việc độ thiết bị lƣợng, phân tích tình hƣ hỏng đảm bảo an toàn khai thác nhà máy điện nguyên tử, phân tích tƣợng xâm thực máy tuốc bin, hệ thống truyền dẫn thuỷ lực, để ứng dụng công nghiệp khai thác, vận chuyển chế biến dầu khí, công nghệ hóa học, Dẫn đến cấp thiết việc nghiên cứu thủy động lực học trình sóng môi trƣờng chất lỏng có bọt Môi trƣờng hỗn hợp chất lỏng với bọt khí hòa tan ngƣng tụ (hay khí không hòa tan không ngƣng tụ) thú vị chỗ chúng đƣợc cấu thành từ ba yếu tố chính: tính phi tuyến, tán sắc trình hao tán lƣợng [1, 21] Bức tranh sóng có nhiều dạng dễ dàng thay đổi cách thay đổi điều kiện thủy động lực, cấu trúc tính chất vật lý nhiệt hỗn hợp với trình tƣơng tác pha Điều đặc biệt hỗn hợp chất lỏng chứa bọt trình thủy động lực xuất lƣợng biến dạng cục hỗn hợp thay đổi thể tích môi trƣờng, thay đổi chủ yếu thay đổi thể tích bọt hỗn hợp tính chất dễ co lại hay giãn nở khí (hoặc hơi) bọt [3], [6] Số hóa trung tâm học liệu Footer Page 10 of 126 http://www.lrc.tnu.edu.vn/ Header Page 33 of 126 3.5 P / Po 2.5 2 1.5 0.005 0.01 0.015 t (s) Hình 3.1:Quá trình lan truyền sóng va đường ống có van đóng trường hợp thể tích pha  20  0.1% ; vận tốc hỗn hợp v 1m / s ,thời gian đóng van hoàn toàn 7m/s.Các đường cong từ 1- profil áp suất biểu diễn tăng áp suất giai đoạn tăng áp, tương ứng với vị trí xác định áp suất hỗn hợp ống xung kích, tính từ van đóng ống là:  = 0; 1.75; 3.51; 4.72; 7.01; 10.15; 13.64 16.00m Từ kết trình bày hình 3.1 Khi thời gian đóng van hoàn toàn ms giá trị lớn áp suất nhận đƣợc van đóng (  = m) 1.72 MPa Sau sóng áp suất lan truyền ngƣợc lại bị giảm dần theo hƣớng trục, mà nguyên nhân ảnh hƣởng quan trọng tƣợng chuyển khối chuyển nhiệt pha Khoảng thời gian bắt đầu xuất hiện tƣợng tăng áp van tới sóng lan truyền tới cuối đƣờng ống (  = 16.00 m) 17.88 ms, nhƣ vận tốc truyền sóng trung bình nhận đƣợc vào khoảng 904.36 m/s Khi sóng tăng áp lan truyền tới cuối đƣờng ống, chênh áp đột ngột ống ống, nên hỗn hợp bắt đầu bị chảy ngƣợc áp suất giảm xuống Sự giảm áp suất chảy ngƣợc bắt đầu lan truyền từ cuối ống van đóng (đây giai đoạn giảm áp sóng va) Số hóa trung tâm học liệu Footer Page 33 of 126 http://www.lrc.tnu.edu.vn/ Header Page 34 of 126 3,5 P / Po 2,5 1,5 0,005 0,01 0,015 0,02 t (s) Hình 3.2:Quá trình lan truyền sóng va đường ống có van đóng trường hợp thể tích pha  20  0.1% ; vận tốc hỗn hợp v 1m / s ,thời gian đóng van hoàn toàn 9m/s.Các đường cong từ 1- profil áp suất biểu diễn tăng áp suất giai đoạn tăng áp, tương ứng với vị trí xác định áp suất hỗn hợp ống xung kích, tính từ van đóng ống là: x = 0; 1.75; 3.51; 4.72; 7.01; 10.15; 13.64 16.00m Từ kết trình bày hình 3.2 Khi thời gian đóng van hoàn toàn ms giá trị lớn áp suất nhận đƣợc van đóng 1.684 MPa Thời gian sóng lan truyền tới cuối đƣờng ống 19.34 ms, nhƣ vận tốc truyền sóng trung bình nhận đƣợc vào khoảng 836.09 m/s 3.5 P / Po 2.5 2 1.5 0.005 0.01 0.015 0.02 t (s) Hình 3.3: Quá trình lan truyền sóng va đường ống có van đóng trường hợp thể tích pha  20  0.1% ; vận tốc hỗn hợp v 1m / s ,thời gian đóng van hoàn toàn 11m/s.Các đường cong từ 1- profil áp suất biểu diễn tăng áp suất giai đoạn tăng áp, tương ứng với vị trí xác định áp suất hỗn hợp Số hóa trung tâm học liệu Footer Page 34 of 126 http://www.lrc.tnu.edu.vn/ Header Page 35 of 126 ống xung kích, tính từ van đóng ống là: x = 0; 1.75; 3.51; 4.72; 7.01; 10.15; 13.64 16.00m Từ kết trình bày hình 3.3 Khi thời gian đóng van hoàn toàn 11 ms giá trị lớn áp suất nhận đƣợc van đóng 1.681 MPa Thời gian sóng lan truyền tới cuối đƣờng ống 20.62 ms, vận tốc truyền sóng trung bình nhận đƣợc vào khoảng 784.19 m/s Nhƣ vậy, tăng dần thời gian đóng van, giá trị áp suất cực đại van đóng thay đổi nhỏ nhƣng nhận đƣợc chậm Tuy nhiên lại có ảnh hƣởng lớn tới vận tốc truyền sóng trung bình hỗn hợp Trong hình 3.4 biểu đồ biểu diễn kết tính toán tăng áp suất hỗn hợp bị đóng van đột ngột Trong trƣờng hợp điều kiện đầu tham số vật chất hỗn hợp đƣợc cho nhƣ sau: thể tích pha  20  0.1% ; vận tốc hỗn hợp v  0.5 m / s ; thời gian đóng van hoàn toàn  *  ms Các đƣờng cong liên tục từ - profil áp suất biểu diễn tăng áp suất giai đoạn tăng áp, tƣơng ứng với vị trí đƣợc xác định áp suất hỗn hợp ống xung kích, tính từ van đóng ống là:  = 0; 1.75; 3.51; 4.72; 7.01; 10.15; 13.64 16.00 m P / Po 1,8 1,6 1,4 1,2 0,005 0,01 0,015 0,02 t (s) Hình 3.4:Quá trình lan truyền sóng va đường ống có van đóng trường hợp thể tích pha  20  0.1% ; vận tốc hỗn hợp v  0.5 m / s ,thời gian đóng van hoàn toàn  *  ms Các đường cong từ 1- profil áp suất biểu diễn tăng áp suất giai đoạn tăng áp, tương ứng với vị trí xác định áp suất Số hóa trung tâm học liệu Footer Page 35 of 126 http://www.lrc.tnu.edu.vn/ Header Page 36 of 126 hỗn hợp ống xung kích, tính từ van đóng ống là:  = 0; 1.75; 3.51; 4.72; 7.01; 10.15; 13.64 16.00 m Từ kết trình bày hình 3.4 Khi thời gian đóng van hoàn toàn ms giá trị lớn áp suất nhận đƣợc van đóng 1.083 MPa Thời gian sóng lan truyền tới cuối đƣờng ống 22.21 ms, nhƣ vận tốc truyền sóng trung bình nhận đƣợc vào khoảng 728.05 m/s Do vậy, giảm vận tốc hỗn hợp áp suất cực đại nhận đƣợc van đóng giảm Trong hình 3.5 biểu đồ biểu diễn kết tính toán tăng áp suất hỗn hợp đóng van đột ngột.Trong trƣờng hợp điều kiện đầu tham số vật chất hỗn hợp đƣợc cho nhƣ sau: thể tích pha  20  0.2 % ; vận tốc hỗn hợp v 1 m / s ; thời gian đóng van hoàn toàn   ms Các đƣờng cong liên tục từ - profil áp suất biểu diễn tăng áp suất giai đoạn tăng áp, tƣơng ứng với vị trí đƣợc xác định áp suất hỗn hợp ống xung kích, tính từ van đóng ống là:  = 0; 0,5; 1,5; 2,5; 3,8; 5,3; 6,3; 7,0; 10,7 mét 2,2 P/Po 1,8 1,6 1,4 1,2 0,005 0,01 0,015 0,02 t(s) Hình 3.5:Quá trình lan truyền sóng va đường ống có van đóng trường hợp thể tích pha  20  0.2 % ; vận tốc hỗn hợp v 1 m / s ,thời gian đóng van hoàn toàn   ms Các đường cong từ 1- profil áp suất biểu diễn tăng áp suất giai đoạn tăng áp, tương ứng với vị trí xác định áp suất hỗn hợp ống xung kích, tính từ van đóng ống là: x = 0; 0,5; 1,5; 2,5; 3,8; 5,3; 6,3; 7,0; 10,7 mét Số hóa trung tâm học liệu Footer Page 36 of 126 http://www.lrc.tnu.edu.vn/ Header Page 37 of 126 Từ kết trình bày hình 3.5 Khi thời gian đóng van hoàn toàn ms giá trị lớn áp suất nhận đƣợc van đóng 1.309 MPa Trong hình 3.6 biểu đồ biểu diễn kết tính toán tăng áp suất hỗn hợp đóng van đột ngột Trong trƣờng hợp điều kiện đầu tham số vật chất hỗn hợp đƣợc cho nhƣ sau: thể tích pha  20  0.2 % ; vận tốc hỗn hợp v 1 m / s ; thời gian đóng van hoàn toàn   9ms Các đƣờng cong liên tục từ - profil áp suất biểu diễn tăng áp suất giai đoạn tăng áp, tƣơng ứng với vị trí đƣợc xác định áp suất hỗn hợp ống xung kích, tính từ van đóng ống là:  = 0; 0,5; 1,5; 2,5; 3,8; 5,3; 6,3; 7,0; 10,7 mét 2,4 2,2 P/Po 1,8 1,6 1,4 1,2 0,005 0,01 0,015 0,02 t(s) Hình 3.6:Quá trình lan truyền sóng va đường ống có van đóng trường hợp thể tích pha  20  0.2 % ; vận tốc hỗn hợp v 1 m / s ,thời gian đóng van hoàn toàn   9ms Từ kết trình bày hình 3.6 Khi thời gian đóng van hoàn toàn ms giá trị lớn áp suất nhận đƣợc van đóng 1.302 MPa Trong hình 3.7 biểu đồ biểu diễn kết tính toán tăng áp suất hỗn hợp đóng van đột ngột.Trong trƣờng hợp điều kiện đầu tham số vật chất hỗn hợp đƣợc cho nhƣ sau: thể tích pha  20  0.2 % ; vận tốc hỗn hợp v 1 m / s ; thời gian đóng van hoàn Số hóa trung tâm học liệu Footer Page 37 of 126 http://www.lrc.tnu.edu.vn/ Header Page 38 of 126 toàn   11ms Các đƣờng cong liên tục từ - profil áp suất biểu diễn tăng áp suất giai đoạn tăng áp, tƣơng ứng với vị trí đƣợc xác định áp suất hỗn hợp ống xung kích, tính từ van đóng ống là:  = 0; 0,5; 1,5; 2,5; 3,8; 5,3; 6,3; 7,0; 10,7 mét 2.2 P/Po 1.8 1.6 1.4 1.2 0.005 0.01 0.015 0.02 t(s) Hình 3.7: Quá trình lan truyền sóng va đường ống có van đóng trường hợp thể tích pha  20  0.2 % ; vận tốc hỗn hợp v 1 m / s ,thời gian đóng van hoàn toàn   11ms Các đường cong từ 1- profil áp suất biểu diễn tăng áp suất giai đoạn tăng áp, tương ứng với vị trí xác định áp suất hỗn hợp ống xung kích, tính từ van đóng ống là:  = 0; 0,5; 1,5; 2,5; 3,8; 5,3; 6,3; 7,0; 10,7 mét Từ kết trình bày hình 3.7 Khi thời gian đóng van hoàn toàn 11 ms giá trị lớn áp suất nhận đƣợc van đóng 1.267 MPa Nhƣ vậy, từ kết biểu diễn hình ta thấy tăng dần thời gian đóng van, giá trị áp suất cực đại van đóng thay đổi nhỏ nhƣng nhận đƣợc chậm Trong hình 3.8 biểu đồ biểu diễn kết tính toán tăng áp suất hỗn hợp đóng van đột ngột.Trong trƣờng hợp điều kiện đầu tham số vật chất hỗn hợp đƣợc cho nhƣ sau: thể tích pha  20  0.2 % ; vận tốc hỗn hợp v  0,5m / s ; thời gian đóng van hoàn toàn   7ms Các đƣờng cong liên tục từ - profil áp suất biểu Số hóa trung tâm học liệu Footer Page 38 of 126 http://www.lrc.tnu.edu.vn/ Header Page 39 of 126 diễn tăng áp suất giai đoạn tăng áp, tƣơng ứng với vị trí đƣợc xác định áp suất hỗn hợp ống xung kích, tính từ van đóng ống là:  = 0; 0,5; 1,5; 2,5; 3,8; 5,3; 6,3; 7,0; 10,7 mét 1,6 1,5 P/Po 1,4 1,3 1,2 1,1 0,005 0,01 0,015 t(s) Hình 3.8: Quá trình lan truyền sóng va đường ống có van đóng trường hợp thể tích pha  20  0.2 % ; vận tốc hỗn hợp v  0,5m / s ,thời gian đóng van hoàn toàn   ms Các đường cong từ 1- profil áp suất biểu diễn tăng áp suất giai đoạn tăng áp, tương ứng với vị trí xác định áp suất hỗn hợp ống xung kích, tính từ van đóng ống là:  = 0; 0,5; 1,5; 2,5; 3,8; 5,3; 6,3; 7,0; 10,7 mét Từ kết trình bày hình 3.8 Khi thời gian đóng van hoàn toàn ms giá trị lớn áp suất nhận đƣợc van đóng 1.046 MPa Trong hình 3.9 biểu đồ biểu diễn kết tính toán tăng áp suất hỗn hợp đóng van đột ngột.Trong trƣờng hợp điều kiện đầu tham số vật chất hỗn hợp đƣợc cho nhƣ sau: thể tích pha  20  0.2 % ; vận tốc hỗn hợp v  0,5m / s ; thời gian đóng van hoàn toàn   9ms Các đƣờng cong liên tục từ - profil áp suất biểu diễn tăng áp suất giai đoạn tăng áp, tƣơng ứng với vị trí đƣợc xác định áp suất hỗn hợp ống xung kích, tính từ van đóng ống là:  = 0; 0,5; 1,5; 2,5; 3,8; 5,3; 6,3; 7,0; 10,7 mét Số hóa trung tâm học liệu Footer Page 39 of 126 http://www.lrc.tnu.edu.vn/ Header Page 40 of 126 1,6 1,5 P/Po 1,4 1,3 1,2 1,1 0,005 0,01 0,015 t(s) Hình 3.9: Quá trình lan truyền sóng va đường ống có van đóng trường hợp thể tích pha  20  0.2 % ; vận tốc hỗn hợp v  0,5m / s ,thời gian đóng van hoàn toàn   9ms Các đường cong từ 1- profil áp suất biểu diễn tăng áp suất giai đoạn tăng áp, tương ứng với vị trí xác định áp suất hỗn hợp ống xung kích, tính từ van đóng ống là:  = 0; 0,5; 1,5; 2,5; 3,8; 5,3; 6,3; 7,0; 10,7 mét Từ kết trình bày hình 3.9 Khi thời gian đóng van hoàn toàn ms giá trị lớn áp suất nhận đƣợc van đóng 1.034 MPa Trong hình 3.10 biểu đồ biểu diễn kết tính toán tăng áp suất hỗn hợp đóng van đột ngột.Trong trƣờng hợp điều kiện đầu tham số vật chất hỗn hợp đƣợc cho nhƣ sau: thể tích pha  20  0.2 % ; vận tốc hỗn hợp v  0,5m / s ; thời gian đóng van hoàn toàn   11ms Các đƣờng cong liên tục từ - profil áp suất biểu diễn tăng áp suất giai đoạn tăng áp, tƣơng ứng với vị trí đƣợc xác định áp suất hỗn hợp ống xung kích, tính từ van đóng ống là:  = 0; 0,5; 1,5; 2,5; 3,8; 5,3; 6,3; 7,0; 10,7 mét Số hóa trung tâm học liệu Footer Page 40 of 126 http://www.lrc.tnu.edu.vn/ Header Page 41 of 126 1,6 1,5 P/Po 1,4 1,3 1,2 1,1 0,005 0,01 0,015 t(s) Hình 3.10: Quá trình lan truyền sóng va đường ống có van đóng trường hợp thể tích pha  20  0.2 % ; vận tốc hỗn hợp v  0,5m / s ,thời gian đóng van hoàn toàn   11ms Các đường cong từ 1- profil áp suất biểu diễn tăng áp suất giai đoạn tăng áp, tương ứng với vị trí xác định áp suất hỗn hợp ống xung kích, tính từ van đóng ống là:  = 0; 0,5; 1,5; 2,5; 3,8; 5,3; 6,3; 7,0; 10,7 mét Từ kết trình bày hình 3.10 Khi thời gian đóng van hoàn toàn 11 ms giá trị lớn áp suất nhận đƣợc van đóng 1.022 MPa Các kết nghiên cứu trình tăng áp suất van đóng dòng chảy hai pha phụ thuộc vào vận tốc hỗn hợp, thời gian đóng van nồng độ thể tích pha tổng kết bảng sau: Vận tốc hỗn Thể tích pha hợp v (m/s)  20  = ms  = ms  = 11 ms 0,1% 1,72 1,684 1,681 0,2% 1,309 1,302 1,267 0,5 0,2% 1,046 1,034 1,022 Áp suất cực đại trên van Pmax (MPa) Bảng 3.1 Số hóa trung tâm học liệu Footer Page 41 of 126 http://www.lrc.tnu.edu.vn/ Header Page 42 of 126 Từ kết đƣợc trình bày bảng 3.1 thấy: - Khi tăng thể tích pha hỗn hợp áp suất cực đại nhận đƣợc van đóng giảm - Khi vận tốc hỗn hợp tăng áp suất cựuc đại nhận đƣợc van đóng tăng - Thời gian đóng van hoàn toàn không ảnh hƣởng nhiều tới áp suất cực đại nhận đƣợc van.(Tuy nhiên thời gian đóng van phải nhỏ nửa pha sóng va ban đầu) Số hóa trung tâm học liệu Footer Page 42 of 126 http://www.lrc.tnu.edu.vn/ Header Page 43 of 126 KẾT LUẬN Luận văn trình bày kết nghiên cứu tăng áp suất giai đoạn tăng áp hỗn hợp nƣớc chứa bọt hơi, hỗn hợp chuyển động ống bị đóng van đột ngột Sự ngừng chảy tăng áp suất van đóng sau đƣợc lan truyền phía ngƣợc lại Các kết tính toán đƣợc kiểm tra cách so sánh chúng với số kết thực nghiệm Fujii T & Akagawa cho thấy phù hợp Trên sở kết nhận đƣợc đƣa số kết luận nhƣ sau: Giá trị áp suất cực đại nhận đƣợc van đóng, phụ thuộc chủ yếu vào thời gian đóng van hoàn toàn, thời gian đóng van nhỏ thời gian nửa pha sóng va ban đầu Tuy nhiên, giá trị cực đại nhận đƣợc chậm tăng thời gian đóng van hoàn toàn Khi tăng thể tích pha hỗn hợp, giá trị cực đại áp suất nhận đƣợc van vận tốc truyền sóng trung bình hỗn hợp giảm Khi vận tốc hỗn hợp tăng, giá trị cực đại áp suất nhận đƣợc van tăng Khi đóng van hoàn toàn, cƣờng độ sóng áp suất hỗn hợp giảm dần, lan truyền ngƣợc lại theo hƣớng trục Nguyên nhân chủ yếu xuất trình tƣơng tác sóng dần tới, tƣợng chuyển khối lƣợng truyền nhiệt pha hỗn hợp Số hóa trung tâm học liệu Footer Page 43 of 126 http://www.lrc.tnu.edu.vn/ Header Page 44 of 126 KIẾN NGHỊ VỀ NHỮNG NGHIÊN CỨU TIẾP THEO Trong nội dung luận văn, dựa sở hệ phƣơng trình thủy 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 trình sóng va môi trƣờng hỗn hợp nƣớc chứa bọt hơi, đƣợc tính chất phi tuyến tính phân tán hao tán lƣợng môi trƣờng đặc biệt Tuy nhiên, tranh toàn cảnh trình sóng va môi trƣờng lỏng - nhiều 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 điều kiện hƣ hỏng đảm bảo an toàn khai thác nhà máy điện nguyên tử, nghiên cứu tƣợng xâm thực ảnh hƣởng sóng va cánh tuốc bin nhà máy thủy điện hay tƣợng sóng va hệ thống truyền dẫn thủy lực,… Đó vấn đề cần thiết nghiên cứu theo hƣớng đề tài Số hóa trung tâm học liệu Footer Page 44 of 126 http://www.lrc.tnu.edu.vn/ Header Page 45 of 126 TÀI LIỆU THAM KHẢO Lan Đau L.D & Lifsitx E.M (2001), Thuỷ động lực học, NXB KH & KT, Hà Nội Dƣơng Ngọc Hải (1987), “Sóng xung kích hỗn hợp dị thể lỏng với chất khí hơi”, Tạp chí Cơ học, 9(4), pp 3- 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 Dƣơng Ngọc Hải, Nguyễn Văn Tuấn (2002), “Hiện tƣợng sóng va hỗn hợp lỏng - hơi”, Tuyển tập công trình hội nghị Cơ học toàn quốc lần thứ 7, Hà Nội, pp 177 - 184 Dƣơng Ngọc Hải, Nguyễn Văn Tuấn (2002), “Sóng phản xạ tƣờng cứng tƣơng tác sóng số hỗn hợp lỏng - hơi”, Tuyển tập Công trình Hội nghị Cơ học toàn quốc lần thứ 7, Hà Nội, pp 165 - 176 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í Toàn quốc, Vũng tàu , 26-28/7/2006 pp.515526 Fujii T et al Water hammers phenomena in one-component two-phase bubbly flow Proc of the Int Conf on Multiphase Flow’ 91- Tsukuba, Japan, 1991 Trần Quốc Nghị (2011), “Nghiên cứu 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 Nguyễn Đức Trung (2012), “Nghiên cứu trình lan truyền, tương tác sóng số 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/2012 10 Duong Ngoc Hai and Nguyen Van Tuan (2002), “Waves reflected by Số hóa trung tâm học liệu Footer Page 45 of 126 http://www.lrc.tnu.edu.vn/ Header Page 46 of 126 solid wall in the mixture of liquid with vapour bubbles”, J Mechanics Vol 24, No 3, pp 167-180 11 Duong Ngoc Hai, Nguyen Van Tuan (1999) , “Shock Waves in some Mixture of Liquid and Gas of two Components”, Trainning-scientific workshop French-Vietnamese, Ha Noi, pp 79 – 86 12 Nguyen Van Tuan, Duong Ngoc Hai, (2004) “ Waves reflected by Solid Wall and Interaction in Vapour Bubbly Liqiuds” The 6th International Conference on Nuclear Thermal Hydraulics, Operations and Safety (NUTHOS-6), Nara, Japan, October 4-8, 2004, Paper ID N6P016 13 Gubaidulin A.A., Beregova O.Sh and Bekishev S.A (1998), “Features of Non-Stationary 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 14 Kubie J (2000), “Velocity of lony Bubbles in Horizontally Oscillating Vertical Piper”, Int J Multiphase Flow Vol 26, pp 339-349 15 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 16 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 1719, 2008, Tokyo, Japan, paper ID N128 17 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 Số hóa trung tâm học liệu Footer Page 46 of 126 http://www.lrc.tnu.edu.vn/ Header Page 47 of 126 18 Fujikawa S and Maerefat M (1991), “A Molecular Theory of the Phase Change from Vapour to Liquid”, Pro of the Int Conf on Multiphase Flow’ 91, Tsukuba, Japan, pp 175-179 19 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, 20 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 21 Nigmatulin (1990), “Dynamics of Multiphase Media”, Hemisphere, publ.Corp, Washington 22 Nguyên Cảnh Cầm, Thủy lực tập 2, nhà xuất giáo dục Số hóa trung tâm học liệu Footer Page 47 of 126 http://www.lrc.tnu.edu.vn/ ... va, trạm thủy điện lớn có đƣờng ống dẫn dài, xây dựng tháp điều áp 2.2 HIỆN TƢỢNG SÓNG VA TRONG HỖN HỢP CHẤT LỎNG CHỨA BỌT HƠI Chất lỏng hai pha hỗn hợp chất lỏng chứa bọt (đƣợc gọi hỗn hợp lỏng. .. THUYẾT Trong phần tác giả trình bày tƣợng nƣớc va (hay gọi tƣợng sóng va) chất lỏng pha (chất lỏng nƣớc) Sau tác giả trình bày mô hình vật lý tƣợng sóng va cho hỗn hợp chất lỏng chứa bọt hơi, hệ... 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 trình tƣơng tác sóng cần thiết Căn vào tình hình phát triển nghiên cứu tƣợng sóng va chất lỏng chứa bọt nƣớc giới, mục tiêu