Tóm tắt Luận án Tiến sĩ Kỹ thuật: Nghiên cứu sóng tràn qua đê biển có tường đỉnh ở Bắc bộ

Chia sẻ: _ _ | Ngày: | Loại File: PDF | Số trang:27

Thêm vào BST

Báo xấu

39
lượt xem 5
download

Download Vui lòng tải xuống để xem tài liệu đầy đủ

Mục tiêu nghiên cứu: Nghiên cứu ảnh hưởng của tường đỉnh thấp đến lưu lượng sóng tràn trung bình và tính chất dòng chảy sóng tràn qua đê biển, từ đó góp phần nâng cao độ tin cậy trong tính toán sóng tràn qua đê biển để bổ sung luận cứ khoa học cho tiêu chuẩn kỹ thuật đê biển hiện nay. Mời các bạn tham khảo!

Chủ đề:

Bình luận(0) Đăng nhập để gửi bình luận!

Lưu

Nội dung Text: Tóm tắt Luận án Tiến sĩ Kỹ thuật: Nghiên cứu sóng tràn qua đê biển có tường đỉnh ở Bắc bộ

BỘ GIÁO DỤC VÀ ĐÀO TẠO BỘ NÔNG NGHIỆP & PTNT TRƯỜNG ĐẠI HỌC THỦY LỢI NGUYỄN VĂN THÌN NGHIÊN CỨU SÓNG TRÀN QUA ĐÊ BIỂN CÓ TƯỜNG ĐỈNH Ở BẮC BỘ Chuyên ngành: Xây dựng công trình thủy Mã số chuyên ngành: 62-58-40-01 TÓM TẮT LUẬN ÁN TIẾN SĨ KỸ THUẬT HÀ NỘI, NĂM 2014
Công trình được hoàn thành tại Trường Đại học Thủy lợi Người hướng dẫn khoa học 1: PGS.TS. Nguyễn Bá Quỳ Người hướng dẫn khoa học 2: GS.TS. Ngô Trí Viềng Phản biện 1: GS.TS. Đinh Văn Ưu, Trường Đại học Khoa học tự nhiên, Đại học Quốc gia Hà Nội Phản biện 2: GS.TS.Trần Đình Hợi, Viện Nước, Môi trường và Biến đổi khí hậu Phản biện 3: GS.TS. Trần Đình Hòa, Viện Khoa học Thủy lợi Việt Nam Luận án sẽ được bảo vệ trước Hội đồng chấm luận án họp tại ............................................................................................................................... ............................................................................................................................... vào lúc ….. giờ …. ngày ….. tháng ….. năm …. Có thể tìm hiểu luận án tại các thư viện: - Thư viện Quốc Gia - Thư viện Trường Đại học Thủy lợi
MỞ ĐẦU 1. Tính cấp thiết của đề tài Việt Nam là quốc gia chịu ảnh hưởng nặng nề của biến đổi khí hậu và nước biển dâng. Dọc theo bờ biển là những trung tâm kinh tế, văn hóa quan trọng của cả nước. Đê biển Bắc bộ có cao trình đỉnh tương đối thấp, bề rộng mặt đê nhỏ, mái đê phía biển, phía đồng dốc và hầu hết đê lại trực diện với biển. Theo kết quả thống kê từ các sự cố vỡ đê trong những năm qua thì sóng tràn gây hư hại mặt đê và mái phía đồng là phổ biến. Một giải pháp hữu hiệu để giảm sóng tràn qua đê là xây tường đỉnh thấp đặt trên đỉnh đê, vì việc tôn cao mặt đê hay làm cơ ở phía thượng lưu là rất tốn kém và khó khả thi, đặc biệt là những tuyến đê có hành lang hẹp. Đến nay, các nghiên cứu ảnh hưởng của tường đỉnh thấp, đặc biệt là nghiên cứu tương tác giữa sóng - tường và dòng chảy sóng tràn chưa đầy đủ. Việc hiểu rõ ảnh hưởng của tường đỉnh thấp và thềm trước tường đối với sóng tràn có ý nghĩa khoa học và thực tiễn, bổ sung luận cứ khoa học cho tiêu chuẩn kỹ thuật đê biển hiện nay. Với ý nghĩa đó tác giả chọn đề tài “Nghiên cứu sóng tràn qua đê biển có tường đỉnh ở Bắc bộ” làm đề tài nghiên cứu của mình. 2. Mục tiêu nghiên cứu Nghiên cứu ảnh hưởng của tường đỉnh thấp đến lưu lượng sóng tràn trung bình và tính chất dòng chảy sóng tràn qua đê biển, từ đó góp phần nâng cao độ tin cậy trong tính toán sóng tràn qua đê biển để bổ sung luận cứ khoa học cho tiêu chuẩn kỹ thuật đê biển hiện nay. 3. Đối tượng và phạm vi nghiên cứu Nghiên cứu sóng tràn qua đê biển có tường đỉnh thấp cho đê biển hiện có ở Bắc bộ - Việt Nam. 1
4. Nội dung nghiên cứu Tổng quan về nghiên cứu sóng tràn qua đê biển có tường đỉnh thấp; mô hình vật lý máng sóng nghiên cứu ảnh hưởng của tường đỉnh thấp đến sóng tràn qua đê biển; tương tác giữa sóng - tường và dòng chảy sóng tràn qua đê biển có tường đỉnh thấp; áp dụng kết quả nghiên cứu tính toán sóng tràn qua đê biển Giao Thủy, tỉnh Nam Định. 5. Cách tiếp cận và phương pháp nghiên cứu 5.1. Cách tiếp cận Để đạt được mục tiêu nghiên cứu, tác giả đã tổng hợp, phân tích các công trình nghiên cứu có liên quan trong nước và trên thế giới về sóng tràn qua đê biển có tường đỉnh thấp. Từ đó lựa chọn hướng tiếp cận vừa mang tính kế thừa vừa mang tính sáng tạo và phù hợp với điều kiện Việt Nam. 5.2. Các phương pháp sử dụng trong luận án Phương pháp nghiên cứu tổng quan; phương pháp nghiên cứu thực nghiệm; phương pháp kết hợp mô hình toán và mô hình vật lý; phương pháp nghiên cứu ứng dụng. 6. Ý nghĩa khoa học và thực tiễn Xây dựng tường đỉnh thấp trên đê nhằm nâng cao cao trình đỉnh đê, giảm thiểu sóng tràn qua đê là giải pháp được áp dụng phổ biến do tính khả thi cao và phù hợp với điều kiện kinh tế hiện nay của nước ta. Tuy nhiên, các nghiên cứu hiện tại về sóng tràn qua đê biển có tường đỉnh thấp chưa đầy đủ, đặc biệt là vấn đề tương tác giữa sóng và tường. Do đó, việc hiểu rõ ảnh hưởng của tường đỉnh thấp đến lưu lượng sóng tràn trung bình cũng như tương tác giữa sóng – tường và tính chất dòng chảy sóng tràn sẽ góp phần nâng cao độ tin cậy trong công tác thiết kế đê biển có tường đỉnh ở nước ta. 2
7. Những đóng góp mới của luận án - Làm sáng tỏ bản chất ảnh hưởng của tường đỉnh đến các đặc trưng sóng tràn và chứng minh được tính ưu việt của thềm trước thông qua việc đi sâu phân tích quá trình tương tác sóng – tường; - Xây dựng được công thức thực nghiệm xác định hệ số ảnh hưởng tổng hợp của tường đỉnh thấp trên đê đến lưu lượng sóng tràn trung bình cho trường hợp sóng đều (2-12); - Xây dựng được đường cong quan hệ tường minh giữa chiều cao sóng bắn với các tham số sóng và hình học tường (Hình 2.13); - Xây dựng được một mặt cắt ngang đê biển tường đỉnh có thềm trước hợp lý, hiệu quả, phù hợp với thực tiễn đê biển Bắc bộ - Việt Nam (Hình 4.8). 8. Cấu trúc của luận án Ngoài phần mở đầu, phần kết luận và kiến nghị; luận án được trình bày trong 4 chương bao gồm: Chương 1: Tổng quan nghiên cứu sóng tràn qua đê biển có tường đỉnh thấp; Chương 2: Mô hình vật lý máng sóng nghiên cứu ảnh hưởng của tường đỉnh thấp đến sóng tràn qua đê biển; Chương 3: Tương tác giữa sóng - tường và dòng chảy sóng tràn qua đê biển có tường đỉnh thấp; Chương 4: Áp dụng kết quả nghiên cứu để tính toán sóng tràn qua đê biển Giao Thủy, tỉnh Nam Định. CHƯƠNG 1 TỔNG QUAN NGHIÊN CỨU SÓNG TRÀN QUA ĐÊ BIỂN CÓ TƯỜNG ĐỈNH THẤP 1.1 Tổng quan về nghiên cứu sóng tràn qua đê biển Hiện nay, do những biến động lớn về môi trường, tác động rõ nét của biến đổi khí hậu toàn cầu, tần suất và cường độ thiên tai ngày càng gia tăng, đặc biệt là 3
bão và nước biển dâng làm gia tăng nguy cơ ngập lụt do vỡ đê, trong đó sóng tràn là một trong những tác nhân cơ bản. Vì vậy, nghiên cứu sóng tràn qua đê biển vẫn là vấn đề thời sự trên Thế giới và ở Việt Nam 1.2 Nguyên nhân, cơ chế phá hoại đê biển và giải pháp giảm thiểu 1.2.1 Nguyên nhân hư hỏng đê biển Có nhiều nguyên nhân gây ra sự cố hư hỏng đê biển trong đó sóng tràn gây phá hoại đỉnh và mái đê phía trong là nguyên nhân cơ bản dẫn đến vỡ đê. 1.2.2 Cơ chế phá hoại đê biển do sóng tràn Cơ chế phá hoại đê biển có nhiều, từ phá hoại cục bộ đến phá hoại tổng thể; các nguyên nhân, yếu tố tác động, hệ quả.. cũng rất đa dạng. Qua phân tích, thống kê, kế thừa có chọn lọc các báo cáo, các sự cố vỡ đê ở Bắc bộ và các kết quả nghiên cứu trong và ngoài nước, thấy rằng: cơ chế phá hoại đê biển do sóng tràn, dòng chảy tràn là cơ chế nổi trội nhất. 1.2.3 Giải pháp giảm thiểu sóng tràn cho đê biển Bắc bộ Hiện nay, có khá nhiều giải pháp để giảm thiểu sóng tràn qua đê, các giải pháp công trình và phi công trình phổ biến như: làm đê ngầm phá sóng ở phía trước đê, dùng các cấu kiện khác nhau đặt trước đê để giảm sóng, làm cơ đê phía biển một bậc hoặc nhiều bậc, trồng rừng ngập mặn phía trước đê, nâng cao cao trình đỉnh đê…Tuy nhiên, do hạn chế về không gian xây dựng và điều kiện kinh tế còn khó khăn nên việc xây tường đỉnh thấp trên đê để nâng cao cao trình đỉnh đê, giảm thiểu sóng tràn qua đê là một trong những giải pháp khả thi về kinh tế và kỹ thuật được sử dụng phổ biến hiện nay ở nước ta. 1.3 Tổng quan về đê biển có tường đỉnh thấp ở Miền Bắc Đê biển có tường đỉnh thấp (W/Hs ≤ 0.5) nằm sát mép ngoài mặt đê được sử dụng khá phổ biến cho hệ thống đê biển ở Miền Bắc và được sử dụng ở những nơi mà không còn quỹ đất để tôn cao, hoặc điều kiện kinh tế không cho phép xây đê cao, như: đê bảo vệ các khu đô thị, khu du lịch. 4
1.4 Tổng quan nghiên cứu sóng tràn qua đê biển có tường đỉnh thấp 1.4.1 Ở trên thế giới (TAW 2002) Ở trên thế giới nghiên cứu phổ biến nhất về sóng tràn qua đê có tường đỉnh thấp là TAW (2002) sau này là EurOtop (2007). Ảnh hưởng của tường đỉnh thấp trên đê đến lưu lượng sóng tràn trung bình theo phương pháp TAW (2002) là không rõ ràng thể hiện qua: ảnh hưởng của tường đỉnh đến lưu lượng sóng tràn trung bình theo phương pháp của TAW (2002) là ẩn; tường đỉnh làm tăng độ dốc mái đê quy đổi có thể là nguyên nhân làm tăng lưu lượng sóng tràn, tuy nhiên sau đó lại được chiết giảm bởi hệ số v. Hệ số ảnh chiết giảm của tường đỉnh v chỉ dựa trên một số ít số liệu thí nghiệm và mới chỉ xét đến ảnh hưởng góc nghiêng của mặt tường phía biển, chưa kể đến tính chất tương tác sóng – tường và dòng chảy sóng tràn cũng như kích thước hình học tường. 1.4.2 Ở Việt Nam Tường đỉnh thấp là dạng công trình đặc thù ở nước ta. Do vậy, có rất ít công trình nghiên cứu về sóng tràn qua dạng công trình này. Nghiên cứu ban đầu của Tuấn và cộng sự (2009) đã đưa ra một cách tiếp cận mới trong việc xem xét ảnh hưởng của tường đỉnh thấp trên đê đến sóng tràn. Tuy nhiên trong nghiên cứu này chưa xem xét ảnh hưởng của thềm trước (S = 0). Tiếp theo là nghiên cứu của Tuấn (2013) đã đánh giá được ảnh hưởng của tường đỉnh và thềm trước tường đến lưu lượng sóng tràn qua đê biển. Tuy nhiên, nghiên cứu của Tuấn (2013) chưa xem xét đến tính chất tương tác sóng - tường và sự thay đổi của dòng chảy sóng tràn qua đê khi có sự hiện diện của tường đỉnh. Ngoài ra, ảnh hưởng chiết giảm của tường đối với sóng đều cũng chưa được đề cập trong nghiên cứu. 1.5 Kết luận chương 1 Sóng tràn là một dạng tải trọng cơ bản nhất đối với đê biển, là tác nhân dẫn đến hơ hỏng và sự cố vỡ đê. Tường đỉnh được xem là giải pháp hiệu quả nâng cao cao trình đỉnh, giảm sóng tràn qua đê. Nghiên cứu sóng tràn qua đê biển có tường đỉnh từ TAW(2002) đến Tuấn (2013) vẫn chưa thật đầy đủ. 5
Tiếp tục phát triển nghiên cứu của Tuấn (2013) luận án đã tiến hành thí nghiệm mô hình vật lý máng sóng với sóng đều nhằm xem xét ảnh hưởng của tường đỉnh đến các đặc trưng sóng tràn và đặc biệt là chi tiết về tương tác sóng – tường và tính chất dòng chảy sóng tràn khi có tường. Sử dụng mô hình toán từ đơn giản như mô hình phi tuyến nước nông (NLSW) đến phức tạp hơn như mô hình (RANS-VOF) để xem xét một cách chi tiết ảnh hưởng của tường dưới góc nhìn tương tác sóng - tường và chế độ dòng chảy sóng tràn qua tường. Những kết quả nghiên cứu của luận án sẽ làm sáng tỏ tính chất sóng tràn qua đê khi có tường đỉnh thấp, góp phần nâng cao chất lượng thiết kế đê biển có tường đỉnh thấp ở nước ta. CHƯƠNG 2 MÔ HÌNH VẬT LÝ MÁNG SÓNG NGHIÊN CỨU ẢNH HƯỞNG CỦA TƯỜNG ĐỈNH THẤP ĐẾN SÓNG TRÀN QUA ĐÊ BIỂN 2.1 Mục đích nghiên cứu Xem xét ảnh hưởng của tường đỉnh thấp đến lưu lượng sóng tràn trung bình và tương tác giữa sóng – tường và tính chất dòng chảy sóng tràn khi có tường. 2.2 Cơ sở lý thuyết về tương tự Để mô hình tương tự với nguyên hình một cách hoàn toàn thì cần phải đầy đủ 3 đặc trưng tương tự: hình học, động học và động lực học. Để có được tương tự cơ bản về các yếu tố sóng, mô hình cần làm chính thái, tỷ lệ mô hình cần tuân theo tiêu chuẩn Froude. Trong các thí nghiệm sóng ngắn với mô hình chính thái tiêu chuẩn Froude tự động được thỏa mãn. 2.3 Mô tả thí nghiệm sóng đều 2.3.1 Máng sóng Máng sóng Hà Lan có tổng chiều dài 45m, chiều dài hiệu quả 42m, chiều cao 1.2m, chiều rộng 1.0m. Máy tạo sóng được trang bị hệ thống hấp thụ sóng phản xạ tự động. 6
Hình 2.1 Toàn cảnh máng sóng sử dụng thí nghiệm Máy tạo sóng có thể tạo sóng đều, hoặc ngẫu nhiên theo một số dạng phổ phổ biến ví dụ như JONSWAP. Chiều cao sóng ngẫu nhiên tối đa có thể tạo ra trong máng là 0.3m và chu kỳ 3.0s (Hình 2.1) 2.3.2 Mô hình đê và các tham số thí nghiệm Mô hình đê và các tham số thí nghiệm có tỷ lệ mô hình 1/10, mái đê nhẵn không thấm nước. Đê trong máng sóng có chiều cao 70cm, mái đê phía biển có độ dốc 1/3. Tường đỉnh thấp trên đê có các chiều cao lần lượt là W = 4cm, 6cm và 9cm được làm bởi các khối rời nhau sao cho khi kết hợp lại với nhau thì chiều cao tường đỉnh thấp (W) và chiều rộng thềm trước (S) được thỏa mãn theo yêu cầu của kịch bản thí nghiệm. Bề rộng thềm trước trong thí nghiệm lần lượt là S = 0cm, 10cm và 20cm. Chiều dài bãi trước đê 24.5m và có độ dốc i =1/100 (Hình 2.5). Hình 2.5 Mô hình thí nghiệm sóng đều 2.3.3 Chương trình thí nghiệm Độ sâu nước máng d = 0.60m được chọn cho thí nghiệm. Một hệ thống 3 đầu đo sóng được đặt ở trước chân đê và một đầu đo được đặt cách chân đê 24.5m. 7
Bảng 2.1 Tổng hợp chương trình thí nghiệm sóng đều Số thí Các thông số sóng Rc (m) W (cm) S (cm) nghiệm H (m) T (s) 40 0.16- 0.24 1.5 – 2.5 0.10 0; 4; 6; 9 0; 10; 20 Một camera có độ phân giải cao được đặt vuông góc với tường kính của máng sóng để thu các ảnh với tốc độ 50 ảnh/s để xem xét tương tác giữa sóng - tường và dòng chảy sóng tràn. Thời gian thí nghiệm là 10 con sóng đối với sóng đều trước khi sóng đều bị ảnh hưởng bởi các sóng phản xạ (Bảng 2.1). 2.3.4 Trình tự thí nghiệm và các tham số đo đạc Thời gian chuẩn bị từ tháng 6 năm 2012 đến tháng 8 năm 2012, thời gian tiến hành thí nghiệm chính thức từ tháng 8 năm 2012 đến tháng 9 năm 2012. Các tham số đo đạc bao gồm: chiều cao sóng H, chu kỳ sóng T, lưu lượng sóng tràn trung bình q, chiều cao sóng bắn Hb, chiều dày lớp nước trên đỉnh tường Ht, độ cao lưu không Rc. 2.4 Phân tích kết quả thí nghiệm 2.4.1 Ảnh hưởng của tường đỉnh đến lưu lượng sóng tràn trung bình Hệ số ảnh hưởng tổng hợp của tường đỉnh là tích của các hệ số ảnh hưởng thành phần do chiều cao tường và chiều rộng thềm trước đem lại. 1 1 1 𝑊 1 𝑆 1 𝛾𝑣 = 𝛾 . 𝛾 = (1 + 𝑐1. 𝑅 . 𝜉) . (1 + 𝑐2. 𝐻 . 𝜉) (2-11) 𝑤 𝑠 𝑐 Xác định hệ số ảnh hưởng tổng hợp 𝛾v (thực đo) và 𝛾v (tính toán), lập đường hồi quy giữa 𝛾v (thực đo) và 𝛾v (tính toán), thử dần các giá trị c1, c2 sao cho đường hồi quy phù hợp nhất với đường phân giác y = x, từ đó tác giả xác định được hệ số c1 =1.26, c2 = 1.44 (Hình 2.12) 1 𝑊 1 𝑆 1 = (1 + 1,26. . ) . (1 + 1,44. . ) (2-12) 𝛾𝑣 𝑅𝑐 𝜉 𝐻 𝜉 2.4.2 Ảnh hưởng của tường đến chiều cao sóng bắn Chiều cao sóng bắn được xác định qua phân tích ảnh bằng Matlab. Tác giả đi xây 𝐻𝑏 𝑆.𝐻 dựng tương quan giữa giữa với , phương trình đường cong có dạng y = 𝐻 𝑔.𝑊.𝑇2 1.544e-30.9x với độ hồi quy R2 = 0.624 từ biểu đồ này ta có thể sơ bộ xác định được chiều cao sóng bắn từ yếu tố sóng, tường (Hình 2.13). 8
1.00 R² = 0.603 0.90 0.80  v thực đo 0.70 0.60 0.50 Sóng vỡ 0.40 Sóng không vỡ 0.30 0.30 0.40 0.50 0.60 0.70 0.80 0.90 1.00 v tính toán Hình 2.10 Hệ số ảnh hưởng tổng hợp của tường đỉnh (đo đạc - tính toán) 2.50 2.00 1.50 y = 1.544e-30.9x R² = 0.624 Hb/H 1.00 0.50 0.00 0.00 0.01 0.01 0.02 0.02 0.03 0.03 0.04 0.04 S.H/g.W.T2 𝐻𝑏 𝑆.𝐻 Hình 2.13 Biểu đồ quan hệ 𝐻 với 𝑔.𝑊.𝑇 2 9
2.5 Kết luận chương 2 Luận án đã xây dựng được công thức thực nghiệm xác định hệ số ảnh hưởng tổng hợp của tường đỉnh đến lưu lượng sóng tràn trung bình trong trường hợp sóng đều (2-12) và đường cong quan hệ một cách rõ ràng giữa chiều cao sóng bắn với các yếu sóng, tường (Hình 2-13). CHƯƠNG 3 TƯƠNG TÁC SÓNG – TƯỜNG VÀ DÒNG CHẢY SÓNG TRÀN QUA ĐÊ BIỂN CÓ TƯỜNG ĐỈNH THẤP 3.1 Đặt vấn đề Ở các mức độ chi tiết khác nhau, luận án đã sử dụng từ mô hình đơn giản như mô hình toán dựa trên phương trình phi tuyến nước nông NLSW (Non-Linear Shallow Water) đến mô hình phức tạp như mô hình toán dựa trên phương trình RANS–VOF (Reynolds Averaged Navier Stokes – Volume Of Fluid). Mô hình NLSW nhìn chung là hiệu quả tính toán cao (ví dụ mô phỏng 1000 con sóng chỉ 5-10 phút). Dùng mô hình này cho kết quả tính toán khá tin cậy lưu lượng sóng tràn trung bình qua đê mái thoải và không có tường đỉnh. Đối với sóng tràn qua các kết cấu có hình dạng phức tạp (ví dụ đê có tường đỉnh thì còn nhiều hạn chế). Luận án đã áp dụng mô hình NLSW của Tuấn và Oumeraci (2010) để mô phỏng, tính toán sóng tràn với sóng ngẫu nhiên. Kết quả tính toán được kiểm chứng với kết quả thí nghiệm từ mô hình vật lý. Mô hình RANS–VOF (máng sóng số) có khả năng áp dụng mô phỏng tương tác sóng - tường và dòng chảy với các dạng kết cấu bất kỳ (tường thẳng đứng, tường rỗng…), từ việc tạo biên sóng nguồn, đến tạo sóng tương tự như máng sóng vật lý. Tuy nhiên, hiệu quả tính toán rất thấp, thường mất nhiều giờ trên máy tính thông thường để mô phỏng một số giây của dòng chảy trong thời gian thực, khó khăn khi áp dụng với sóng ngẫu nhiên do thời đoạn tính toán dài. Vì những lý do nêu trên mô hình NLSW của Tuấn và Oumeraci (2010) được sử dụng để kiểm chứng kết quả dự báo lưu lượng sóng tràn trung bình. Kết quả tính toán mô hình toán so sánh với các dữ liệu thu được từ thí nghiệm mô hình vật lý 10
cho sóng ngẫu nhiên. Sự tương tác giữa sóng - tường và dòng chảy, hay nói cách khác các đặc trưng của dòng chảy sóng tràn được xem xét một cách chi tiết cho trường hợp sóng đều bằng cả mô hình vật lý và mô hình máng sóng số. 3.2 Mô hình NLSW (Tuấn và Oumeraci, 2010) 3.2.1 Hệ phương trình cơ bản Mô hình dựa trên dạng bảo toàn của các phương trình NLSW được giải theo sơ đồ của Roe kết hợp với phương pháp xác định thông lượng bậc cao TVD (total variation diminishing), như sau: U  F ( x,U )   S ( x, U ) (3-1) t x Trong đó: các véc-tơ bảo toàn U , F ( x,U ) và các vec-tơ thành phần S ( x,U ) được xác định: h  U ( x)    (3-2)  uh   uh  F ( x, U )   2  (3-3)  u h  gh / 2  2 0  S ( x, U )    (3-4)  gh( Sbx  S f  Sr )  Trong đó: g là gia tốc trọng trường, h là độ sâu dòng chảy, u là vận tốc dòng chảy theo phương ngang, Sbx là mái dốc đáy, Sf là ma sát đáy, Sr là độ dốc do tiêu năng cuộn sóng mặt. 3.2.2 Sóng tràn đối với sóng ngẫu nhiên Tường thẳng đứng không mô tả được trong mô hình NLSW vì trường dòng chảy khi đó sẽ vi phạm các giới hạn nước nông. Do vậy việc điều chỉnh, biến đổi hình học tương đương là cần thiết. Ở đây tác giả sử dụng hai phương án biến đổi đó 11
là: phương án tường tương tương đương (Hình 3.1) và phương án chiều cao lưu không tương đương (Hình 3.2). Sử dụng mực nước trung bình kết hợp với chuỗi sóng theo thời gian đã được đo đạc từ các đầu đo sóng đặt gần với chân đê nhất (biên sóng đầu vào càng gần càng đảm bảo độ chính xác cho mô hình NLSW). Hình 3.1 Mô tả tường thẳng đứng qua mái nghiêng (TAW-2002)-PA1 Hình 3.2 Mô tả tường bằng chiều cao lưu không tương đương – PA2 Tại phía hạ lưu, giá trị mực nước không đổi thấp hơn hẳn đỉnh đê (để tránh bất kỳ ảnh hưởng nào đến sóng tràn) được sử dụng biên đầu ra. Thời gian mô phỏng giống như thí nghiệm mô hình vật lý (1000.Tp ~ 10 phút PC). Nhìn chung, các kết quả của cả hai phương án phù hợp khá tốt với các số liệu thí nghiệm mô hình với độ hồi quy R2 lần lượt là 0.88 và 0.87 cho phương án thứ nhất và phương án thứ hai. Sai số trung bình là 39.8% với độ lệch chuẩn là  56.2%. Tuy nhiên, vẫn tồn tại sai khác khá lớn cho một số trường hợp cụ thể với lưu lượng sóng tràn trung bình nhỏ khi tường đỉnh cao và khi tường không có thềm trước. Điều này 12
là do sự tương tác phức tạp giữa sóng và tường không thể được giải quyết một cách đầy đủ trong mô hình NLSW thông qua phương pháp biến đổi tương đương. Hình 3.3 Kết quả tính toán sóng tràn bằng mô hình NLSW (TAW 2002)-PA1 Hình 3.4 Kết quả tính toán sóng tràn bằng mô hình NLSW (Rc*)-PA2 13
3.3 Mô hình RANS-VOF (COBRAS-UC, máng sóng số) 3.3.1 Giới thiệu máng sóng số Máng sóng số có khả năng mô phỏng tương tác sóng - công trình với các tính năng tương tự như máng sóng vật lý. Dạng công trình có thể mô phỏng là bất kỳ với kết quả mô phỏng các đặc trưng dòng chảy khá chi tiết và đầy đủ. Máng sóng số đã được kiểm nghiệm với nhiều số liệu thí nghiệm mô hình vật lý tỷ lệ lớn, nhỏ và cho kết quả tin cậy. 3.3.2 Hệ phương trình cơ bản Mô hình dựa trên hệ phương trình trung bình Reynolds Navier – Stockes 2 chiều  ui 0 (3-12) xi  ui  ui 1 p 1    ui   uj   gi     uiu j  (3-13) t x j  xi  x j  x j  Và được khép kín bởi hệ phương trình vận chuyển rối k k    k   ui  uj   t      uiu j  (3-14) t x j x j   k x j   x j      t      ui  2 (3-15)  uj       C1 uiu j  C2 t x j x j     x j  k x j k Trong đó: ui là vận tốc trung bình theo phương i (i, j =1, 2 cho dòng chảy hai chiều), p là áp suất dòng chảy,  khối lượng riêng của nước, gi là gia tốc trọng trường theo phương i,  uiuj là ứng suất Reynolds được mô phỏng theo độ nhớt xoáy phi tuyến, Các hệ số kinh nghiệm của mô hình rối k = 1.0,  =1.3, C1 = 1.44, C1 = 1.92;  = / vàt = Cdk2/ (Cd = 0.99) tương ứng là các hệ số nhớt động và nhớt xoáy. COBRAS-UC tính toán dòng chảy trên một lưới chữ nhật không đều, mặt thoáng bất kỳ của dòng chảy được tính toán theo phương pháp thể tích chất lỏng (VOF). 14
3.3.3 Sóng tràn đối với sóng ngẫu nhiên Chuỗi sóng ít nhất 1000 con sóng (1000.Tp = 2200s) cần 75h để mô phỏng với máy tính có cấu hình 3.1GHz-4GB RAM. Do hạn chế về năng lực máy tính nên chỉ có 14 trường hợp thí nghiệm tiêu biểu được xem xét. Đây là sự kết hợp của một điều kiện sóng (Hm0 = 0.10m, Tp = 2.2s) và chiều sâu mực nước d = 0.55m với tất cả các dạng đê với có tường đỉnh cao (W= 6cm và 9cm, có hoặc không có thềm trước). Hình 3.7 Lưu lượng sóng tràn trung bình (Sóng ngẫu nhiên, COBRAS-UC) Lưu lượng sóng tràn trung bình từ thí nghiệm mô hình vật lý được so sánh với kết quả tính toán bởi mô hình COBRAS-UC và mô hình NLSW với cả hai phương án biến đổi tường tương đương. Kết quả từ mô hình NLSW và dữ liệu thí nghiệm có sai số trung bình cho phương án 1 là 60.1%, độ lệch chuẩn  63.2% và phương án 2 là 129.4%, độ lệch chuẩn  100.6%. Kết quả từ mô hình COBRAS-UC và dữ liệu từ thí nghiệm khá khớp với nhau với độ sai số trung bình chỉ là 39.7%, độ lệch chuẩn  24.5%. Rõ ràng, mô hình COBRAS-UC thể hiện ưu điểm hơn so với mô hình NLSW. Tuy vậy, kết quả từ mô hình COBRAS- 15
UC vẫn có thể cho sai số lên tới 63% cho các trường hợp lưu lượng sóng tràn nhỏ (Hình 3.7). 3.3.4 Sóng tràn đối với sóng đều 3.3.4.1 Lưu lượng sóng tràn trung bình Với 40 thí nghiệm cho 10 mô hình vật lý tương ứng, các chuỗi sóng đều được mô phỏng để đánh giá lưu lượng sóng tràn trung bình cho mỗi thí nghiệm. Thời gian tính toán cho mỗi thí nghiệm khoảng 6 giờ được thực hiện trên PC chuẩn (1giờ CPU/10 giây thời gian dòng chảy). Lưu lượng sóng tràn trung bình tính toán theo mô hình COBRAS-UC được so sánh với dữ liệu đo đạc từ thí nghiệm mô hình vật lý. Mô hình COBRAS-UC dự báo khá tin cậy về lưu lượng sóng tràn trung bình qua các dạng tường đỉnh thấp khác nhau (R2 = 0.95). Sai số trung bình chỉ là 23.4% với độ lệch chuẩn là  30.2 %. Sai số đáng kể nhất xảy ra trong một số ít trường hợp khi chiều cao tường lớn (W=9cm, hình tam giác). So sánh với các trường hợp sóng ngẫu nhiên cho thấy sự phù hợp của kết quả tính toán cho sóng đều là tốt hơn (Hình 3.8). Hình 3.8 Lưu lượng sóng tràn trung bình (Sóng đều, COBRAS-UC) 16
3.3.4.2 Tương tác giữa sóng - tường và dòng chảy Sóng tràn bao gồm hai thành phần: thành phần tràn chảy thành dòng (green overtoping) và thành phần sóng bắn (splash overtoping). Để hiểu về ảnh hưởng của tường đỉnh đến sóng tràn, ảnh phân tách từ máy quay được xem xét và phân tích tỉ mỉ cùng với kết quả đạt được từ mô hình toán. Quá trình tương tác sóng – tường được mô tả trong 4 giai đoạn như sau: Giai đoạn 1: Quá trình sóng tràn bắt đầu bằng việc lưỡi sóng va với tường và kết quả là tạo ra sóng bắn trong không khí (Hình 3.13); Giai đoạn 2: Sóng bắn bắt đầu biến mất khi đạt tới một độ cao nhất định trên tường và đổ xuống đỉnh tường, mặt đê. Tại thời điểm cuối của giai đoạn này, một dòng chảy tràn được hình thành (Hình 3.14); Giai đoạn 3: Nếu sóng tiếp tục tràn mạnh trong giai đoạn thứ 3 này, dòng chảy tràn sẽ phát triển tới chiều sâu ngập lớn nhất bên trên đỉnh tường (Hình 3.15); Giai đoạn 4: Sóng chảy tràn theo dòng tiếp tục phát triển cho tới khi sóng rút về phía biển trong giai đoạn cuối cùng (Hình 3.16). Tất cả các thời đoạn tính toán cũng được xem xét kỹ lưỡng để xác định chiều cao sóng bắn lớn nhất bằng mô hình toán và so sánh với kết quả đạt được trong mô hình vật lý. Xét về mặt mô hình toán, mô hình COBRAS-UC có khả năng mô phỏng sóng bắn rời nhưng khi miêu tả các giọt nước riêng biệt hay các tia nhỏ như trong mô hình vật lý thì có độ tin cậy không cao. Phương pháp mô phỏng bề mặt thoáng (VOF) trong mô hình toán không kể đến sức căng bề mặt cũng như việc chia lưới quá thô so với kích thước của từng hạt nước riêng rẽ là các nguyên nhân dẫn tới sự khác biệt giữa kết quả mô hình toán với mô hình vật lý, các hạn chế này chưa thể được khắc phục trong công nghệ mô phỏng dòng chảy hiện nay. Quá trình 4 giai đoạn của sự tương tác giữa sóng - tường từ mô vật lý cho trường hợp thí nghiệm REW6S20_4, được kiểm chứng một cách chi tiết bằng mô hình COBRAS-UC tương ứng. Hai giai đoạn đầu (Hình 3.13, Hình 3.14) kết quả không chính xác bằng hai giai đoạn sau (Hình 3.15, Hình 3.16). 17
Hình 3.13 Sóng bắn khi sóng va vào Hình 3.14 Sóng đổ lên đỉnh tường tường t= 27.1s ( MH vật lý) t=27.3s ( MH vật lý) Hình 3.15 Sóng chảy thành dòng Hình 3.16 Sóng rút t=27.8s t = 27.5s ( MH vật lý) ( MH vật lý) Ở mức độ ít chi tiết hơn, Hình 3.17, Hình 3.18 lần lượt thể hiện các kết quả dự báo bằng mô hình toán về chiều cao sóng bắn lớn nhất (xảy ra trong giai đoạn 1, Hình 3.13) và chiều sâu dòng chảy tràn lớn nhất trên đỉnh tường (xảy ra trong giai đoạn 3, Hình 3.15) được so sánh với các dữ liệu từ thí nghiệm mô hình vật lý. Khi xem xét kết quả từ mô hình vật lý và kết quả dự báo từ mô hình toán (Hình 3.17) với cùng một điều kiện sóng và chiều cao tường thì chiều cao sóng bắn lớn nhất giảm khi chiều rộng thềm trước tường tăng lên. Cũng tương tự, với cùng một điều kiện sóng và chiều rộng thềm thì khi chiều cao tường lớn hơn sẽ cho kết quả sóng bắn cao hơn. 18