Phát triển mô hình sóng-dòng qua cặp mô hình FEM và SWAN tại vịnh Nha Trang

Mô hình sóng SWAN (phiên bản 41.31) với lưới phi cấu trúc đã được áp dụng thử nghiệm vào tính sóng tại vùng biển Nha Trang. Đặc biệt, có thể sử dụng cùng mạng lưới tam giác phi cấu trúc trong tính toán dòng chảy bằng phương pháp phần tử hữu hạn (FEM). Bài viết này trình bày về việc phát triển mô hình sóng-dòng qua cặp mô hình FEM và SWAN tại vịnh Nha Trang.

52 T V Chung, N M Tiến, Võ Văn Quang / Tạp chí Khoa học Cơng nghệ Đại học Duy Tân 06(43) (2020) 52-59 06(43) (2020) 52-59 Phát triển mơ hình sóng-dịng qua cặp mơ hình FEM SWAN vịnh Nha Trang Development of a wave-current model through coupling of FEM and SWAN models in Nha Trang bay Trần Văn Chung*, Ngô Mạnh Tiến, Võ Văn Quang Tran Van Chung*, Ngo Manh Tien, Vo Van Quang Viện Hải dương học, Viện Hàn lâm Khoa học Công nghệ Việt Nam (VAST) Institute of Oceanography, VAST (Ngày nhận bài: 23/11/2020, ngày phản biện xong: 05/12/2020, ngày chấp nhận đăng: 20/12/2020) Tóm tắt Mơ hình sóng SWAN (phiên 41.31) với lưới phi cấu trúc áp dụng thử nghiệm vào tính sóng vùng biển Nha Trang Đặc biệt, sử dụng mạng lưới tam giác phi cấu trúc tính tốn dịng chảy phương pháp phần tử hữu hạn (FEM) Một kết hợp hệ thống sóng-dịng phát triển báo kết hợp hai mơ hình, mơ hình thủy động lực theo phương pháp phần tử hữu hạn (FEM) mơ hình sóng gần bờ (SWAN) Điều hữu ích giải tốn tương tác sóng - dịng sóng - sóng Những kết bước đầu cho thấy, cách tiếp cận ứng dụng mơ hình tính sóng hợp lý, ổn định cho bước thời gian cho mạng lưới làm mịn, mang đặc trưng địa phương vùng nghiên cứu Từ khóa: Mơ hình kết hợp; SWAN; Lưới phi cấu trúc; Sai phân hữu hạn; Phần tử hữu hạn Abstract The wave model SWAN (version 41.31) with unstructured grid has been applied for Nha Trang waters In particular, the unstructured meshes (triangular mesh) can be used in order to calculate the flow by finite element method (FEM) A coupled wave-current system has been developed in the present paper based on two open source community models, the Hydrodynamic model by finite element method (FEM) and the Simulating WAves Nearshore (SWAN) This is very helpful in solving the problems of interactions between wave-current and wave-wave It indicates that the SWAN model approach is reasonable, stable for any time step while permitting local mesh refinements in interested areas Keywords: Coupled model; SWAN; Unstructured grid; Finite difference; Finite Element Mở đầu Sóng dịng chảy thường tương tác rõ rệt vùng nước nông điều kiện trường độ sâu phức tạp đóng vai trị vận chuyển trầm tích, tiến hóa hình thái pha trộn chất nhiễm (Rodriguez nnk., 1995 [1]; Li Johns, 1998 [2]; Bever MacWilliams, 2013 [3]) Các tương tác phi tuyến phức tạp (Olabarrieta nnk., 2011 [4]; Roland nnk., 2012 [5]; Benetazzo nnk., 2013 [6]) Một mặt, gradient ứng suất xạ kết hợp với sóng trọng lực, tổng lực đẩy * Corresponding Author: Tran Van Chung; Institute of Oceanography, VAST Email: tvanchung@gmail.com T V Chung, N M Tiến, Võ Văn Quang / Tạp chí Khoa học Công nghệ Đại học Duy Tân 06(43) (2020) 52-59 53 ngang tác động sóng lên bãi biển khu vực gần bờ làm thay đổi dòng chảy dọc bờ (Longuet-Higgins, 1970 [7]; Garcez-Faria nnk., 2000 [8]) ảnh hưởng đến mực nước trung bình cách tăng mực nước gần bờ giảm mực nước gần điểm sóng đổ (Longuet-Higgins Stewart, 1964 [9]; Guza Thornton, 1981 [10]) Trong đó, sóng vùng sóng đổ tăng cường pha trộn theo phương ngang lực cản đáy Mặt khác, thay đổi mực nước dịng chảy ảnh hưởng đến chuyển động phát triển sóng (Dutour-Sikiric nnk., 2013 [11]; Allard nnk., 214 [12]) Bài báo thực ứng dụng mơ hình SWAN (phiên 41.31, cập nhật tới thời điểm này) với lưới phi cấu trúc vào vùng nghiên cứu có đường bờ biến đổi phức tạp, biên mở rộng, xử lý biên phức tạp không thuận lợi cho lan truyền sóng Phiên sử dụng khác đơi chút so với lưới phi cấu trúc với kỹ thuật lặp lặp lại bốn hướng Gauss-Seidel tương tự từ phiên cấu trúc SWAN, đòi hỏi thích nghi lõi tính tốn Điểm bật thuật tốn lưới phi cấu trúc khơng theo phương pháp thể tích hữu hạn phương pháp phần tử hữu hạn theo phương pháp sai phân hữu hạn truyền thống Với lộ trình thực đây, mơ hình giữ trình vật lý số học cấu trúc mã số mơ hình lưới cấu trúc SWAN, chạy mạng lưới phi cấu trúc Phiên cập nhật từ trang web: http://swanmodel.sourceforge.net/ Các phiên SWAN ứng dụng thành công kết hợp với mơ hình thủy động lực theo phương pháp phần tử hữu hạn (FEM) Các cơng trình cơng bố điển hình, áp dụng mơ hình SWAN với lưới phi cấu trúc áp dụng thử nghiệm vào tính tốn đặc trưng sóng vùng biển Đầm Bấy (Nha Trang) có xét đến tương tác sóng - dịng theo pha triều lên pha triều xuống (Trần Văn Chung nnk., 2015 [13]) trường sóng ảnh hưởng chế độ gió mùa Ninh Thuận - Bình Thuận (Trần Văn Chung nnk., 2016 [14]) Đây phiên ứng dụng tương tác sóng - dịng ảnh hưởng chế độ gió mùa khu vực khơng thuận lợi cho lan truyền sóng Tài liệu phương pháp 2.1 Mơ tả mơ hình 2.1.1 Mơ hình SWAN Các mơ hình sóng sử dụng nghiên cứu mơ hình sóng hệ thứ ba SWAN 41.31 (Mơ sóng gần bờ (SWAN Simulating Waves Nearshore); mơ hình phổ sóng SWAN tính tốn phát triển phổ mật độ tác động N sử dụng phương trình cân tác động (Booij nnk., 1999 [15]): (1) Với Stot = Sin + Swc + Sn14 + Sbot + Sbrk + Snl3 (2) Vế bên trái phương trình (1), số hạng biểu diễn thay đổi mật độ tác động sóng theo thời gian, lan truyền sóng theo địa lý khơng gian (với véc tơ vận tốc nhóm sóng – véc tơ dòng chảy), khúc xạ độ sâu dòng chảy gây (với cθ vận tốc lan truyền theo hướng không gian θ) chuyển dịch tần số radian σ thay đổi lấy trung bình dịng chảy độ sâu (với vận tốc lan truyền cσ) Vế bên phải biểu diễn cho trình thành tạo, tiêu tán phân phối lại lượng sóng Trong nước sâu, ba số hạng nguồn phát sử dụng Đây chuyển giao lượng từ gió đến sóng, Sin, tiêu tán lượng sóng sóng bạc đầu, Swc, chuyển đổi phi tuyến lượng sóng tương tác bốn (bốn sóng), 54 T V Chung, N M Tiến, Võ Văn Quang / Tạp chí Khoa học Công nghệ Đại học Duy Tân 06(43) (2020) 52-59 Snl4 Trong vùng nước nông, tiêu tán ma sát đáy, Sbot, độ sâu gây đổ vỡ, Sbrk, ba tương tác phi tuyến (ba sóng), Snl3, cho ước lượng thêm vào Chi tiết mở rộng công thức q trình tìm thấy trong: Ris (1997) [16], Booij nnk (1999) [15] Holthuijsen (2007) [17] Đối với toán đặt hợp lý, điều kiện biên phải cung cấp Các thành phần sóng đến biên phía biển quy định phổ hai chiều Tại biên khép kín, ví dụ: biên đường bờ biển biên bên, lượng sóng tiêu tán hấp thụ hồn toàn giữ lại vùng địa lý tương ứng Các biên không gian tần số định tương ứng σmin σmax Các biên hấp thụ đầy đủ, phần chẩn đốn σ-4 thêm vào tần số cắt cụt cao, sử dụng để tính tốn thành phần phi tuyến tương tác sóng-sóng tính tốn tồn thơng số sóng Từ định hướng khơng gian vùng vịng trịn khép kín, khơng cần thiết có điều kiện biên 2.1.2 FEM model Các cơng trình nghiên cứu số trị cho mơ hình thủy động lực học theo phương pháp phần tử hữu hạn (FEM), xây dựng triển khai cơng trình nghiên cứu Bùi Hồng Long Trần Văn Chung (2009 [18]) Đánh giá tính hợp lý mơ hình qua so sánh với số liệu thực tế trạm mực nước nghiên cứu chế độ dòng chảy cho Vịnh Bắc Bộ (Trần Văn Chung Bùi Hồng Long (2015) [19]) Trong đó, đáng ý bước đầu thực so sánh kết nghiên cứu cấu trúc dịng chảy hai mơ hình FEM mơ hình Ecosmo kiểm chứng thực tế đo đạc vùng nghiên cứu Bình Cang - Nha Trang (Trần Văn Chung Bùi Hồng Long (2014)) [20] sở liệu dòng chảy đầu vào quan trọng cho việc chạy liên kết tương tác sóng- dịng mơ hình cặp FEM-SWAN 2.2 Cơ sở liệu Cơ sở liệu mơ hình gồm: (1) Trường độ sâu: Trường độ sâu cập nhật khuôn khổ đề tài tỉnh Khánh Hịa: “Xác định khu vực có khả cải tạo, phát triển bãi tắm nhân tạo đề xuất phương án bảo vệ bãi tắm tự nhiên vịnh Nha Trang” (2015 2016) với tỉ lệ 1/25.000 cho tồn vùng 1/10.000 cho ven bờ (2) Thơng tin trường chảy dịng chảy tầng mặt trích xuất từ mơ hình tính thủy động lực ba chiều theo phương pháp phần tử hữu hạn cho vùng biển Bình Cang - Nha Trang, với mạng lưới tính trùng với mạng lưới tính sóng theo cơng trình nghiên cứu Trần Văn Chung Bùi Hồng Long công bố năm 2014 [20] mô sở cập nhật liệu nghiên cứu (3) Thông tin đặc trưng sóng (độ cao, chu kỳ hướng) biên khơi cập nhật từ: http://oos.soest.hawaii.edu/erddap/griddap/NWW3_Global_Best.html với độ phân giải 0,5 độ theo (số liệu tính thống kê từ ngày 07/11/ 2010 đến tháng 31/12/2019) (4) Nguồn số liệu vận tốc gió cập nhật từ sở liệu tái phân tích mơ hình khí hậu tồn cầu CFSR (Climate Forecast System Reanalysis) thuộc trung tâm dự báo môi trường NCEP (National Centers for Environmental Prediction) (NCEP CFSR) với tần suất 01 số liệu Với phiên (NCEP CFSv2 - NCEP Climate Forecast System Version (CFSv2)): số liệu lấy từ năm 01/01/2011 đến 31/12/2019 theo lưới 0,20 theo phương ngang Ngồi ra, để phân tích số liệu gió địa phương trạm Nha Trang, báo sử dụng nguồn số liệu gió trạm đo khí tượng với chuỗi số liệu đo từ 1987 đến 2007, tần suất đo số liệu lần ngày 1, 7, 13, 19 Dưới vài kết phân tích từ liệu gió có được: T V Chung, N M Tiến, Võ Văn Quang / Tạp chí Khoa học Công nghệ Đại học Duy Tân 06(43) (2020) 52-59 55 N NNW Tốc độ gió (m/s) NNE NW NE WNW ENE W 2 - >4 - >6 - >8 - 10 >10 - 12 >12 - 14 >14 - 16 >16 - 18 >18 - 20 >20 E 0% 5% 10% 15% 20% WSW 25% ESE SW SE SSW SSE S Hình Hoa gió khu vực Nha Trang Theo kết phân tích gió cho 21 năm khu vực Nha Trang, thấy rõ đặc trưng gió mang tính địa phương với biến đổi tốc độ hướng lớn so với chế độ gió mùa chung, cụ thể gió mùa Tây Nam với hướng gió Nha Trang Đơng Nam (SE) chiếm ưu Theo hình 6, thấy rằng, khu vực Nha Trang thường xun có hướng gió (N, NNE, NE, SE) Căn vào kết phân tích chi tiết cho cụ thể năm, trường gió mùa Đơng Bắc tác động đến khu vực nghiên Hình Mạng lưới nghiên cứu tác động sóng, dịng cho vịnh Nha Trang cứu, vịnh Nha Trang chịu tác động hướng gió Đơng Bắc, Bắc Đơng Bắc, Đơng Đơng Bắc Bắc Khi trường gió mùa Tây Nam xuất hiện, vịnh Nha Trang chịu tác động hướng gió Đơng Nam, Đơng Đơng Nam Nam Đông Nam 2.3 Khu vực nghiên cứu Khu vực nghiên cứu chọn từ kinh độ 109,141oE đến 109,321oE; vĩ độ từ 12,125oN đến 12,462oN (Hình 2) Phương pháp giải sử dụng phương pháp phần tử hữu hạn, mạng lưới tính mạng lưới tam giác (Hình 3) Với mạng lưới tính tam giác thiết lập với góc cực tiểu 25o; tổng diện tích mặt thống cho tính tốn 354,28 km2 Số điểm tính mạng lưới tam giác 7.421, với tổng số tam giác 13.908 Diện tích tam giác lưới tính có giá trị nhỏ 4.714m2, trung bình 25.473m2, lớn 39.999m2 Hình Mạng lưới làm khớp đồ Google Earth Hình Sơ đồ trạm đo kiểm định hiệu chỉnh mơ hình, “LT” trạm đo liên tục 01 ngày đêm Dòng chảy sóng đo đồng thời máy AWAC (Acoustic WAve and Current profiler) số liệu đề tài AST06.03/18-19, đo vào tháng 9/2018) 56 T V Chung, N M Tiến, Võ Văn Quang / Tạp chí Khoa học Cơng nghệ Đại học Duy Tân 06(43) (2020) 52-59 Kết nghiên cứu Kết tính tốn đặc trưng sóng theo dạng tương tác sóng tổng hợp với tác động đồng thời trường gió, dịng chảy tầng mặt độ cao bề mặt biển Mô đầy đủ tương tác làm sáng tỏ chế độ sóng thực tế vùng biển Nha Trang, theo chế độ gió mùa với tác động tổng thể lực tạo sóng quan trọng Từ có nhìn tổng quan tác động sóng ảnh hưởng tương tác qua lại sóng - dịng q trình thủy động lực vịnh 3.1 Trường sóng điển hình mùa gió Đơng Bắc Hình Phân bố dịng chảy tầng mặt vào mùa gió Đơng Bắc Về hình dạng phân bố có phân vùng độ cao sóng rõ rệt, phía bắc vịnh Nha Trang chịu tác động rõ trường gió Đơng Bắc điểm che chắn phía nam vịnh làm suy yếu trường sóng vào phía nam vịnh Sự suy giảm sóng tác động che chắn địa hình rõ ràng, chi tiết lan truyền độ cao sóng thể Hình 6, chu kỳ sóng Hình hướng sóng lan truyền Hình Điểm đặc biệt cần lưu ý, sóng phía ngồi bắc vịnh tương đối lớn mức độ ảnh hưởng vào vùng ven bờ bị tiêu tán suy giảm Hình kết phân tích dịng chảy trung bình tầng mặt khu vực Nha Trang vào mùa gió Đơng Bắc mơ hình FEM Đây thơng tin điều kiện ban đầu cho dịng chảy mơ trường gió mùa Đơng Bắc, nghiên cứu tương tác sóng - dịng mơ hình SWAN Ở cơng trình này, chúng tơi tập trung vào phân tích tính hợp lý q trình lan truyền sóng ảnh hưởng dịng chảy Bộ chương trình thể kết đặc trưng sóng (độ cao, chu kỳ hướng sóng) cải tiến việt hóa từ chương trình thể kết tính SWAN cho lưới phi cấu trúc chương trình mã nguồn mở Matlab “plotunswan.m” Hình Độ cao sóng có nghĩa (m) tương tác sóng dịng ảnh hưởng mùa gió Đơng Bắc lượng sóng mạnh nên mức độ ảnh hưởng khơng cịn lớn phía ngồi khơi Với tác động yếu tố dòng chảy mực nước hướng lan truyền sóng có đổi hướng so với hướng gió túy, kết mô cho thấy giá trị độ cao sóng lớn mạng lưới tính đạt 2,18m với chu kỳ sóng 6,3s (độ cao sóng trung bình 1,0m, chu kỳ trung bình 3,7s) có độ sâu khoảng 26,1m, với hướng lan truyền sóng 224,7o, vị trí đạt độ cao sóng lớn xung quanh vị trí (109,311996oE; 12,273005oN) (hình 6-8) T V Chung, N M Tiến, Võ Văn Quang / Tạp chí Khoa học Công nghệ Đại học Duy Tân 06(43) (2020) 52-59 57 Hình Chu kỳ trung bình (s) tương tác tổng hợp ảnh hưởng mùa gió Đơng Bắc Hình Hướng lan truyền sóng (0) tương tác tổng hợp ảnh hưởng mùa gió Đơng Bắc 3.2 Trường sóng điển hình mùa gió Tây Nam Hình Phân bố dịng chảy tầng mặt vào mùa gió Tây Nam Theo kết mơ phỏng, mùa gió Tây Nam với sóng tác động đến vịnh Nha Trang không đáng kể Kết mô thể chi tiết độ cao lan truyền sóng có nghĩa Hình 10, chu kỳ sóng Hình 11 hướng lan truyền Hình 12 Nói chung, vịnh Nha Trang vịnh kín nên yếu tố tác động sóng ảnh hưởng đến biến đổi đường bờ khơng lớn Hình 10 Độ cao sóng có nghĩa (m) tương tác sóng - dịng ảnh hưởng mùa gió Tây Nam Thơng tin sóng mùa gió Tây Nam đáng ý sau: độ cao sóng có nghĩa lớn mạng lưới tính đạt 1,4m, chu kỳ sóng 5,1s (trung bình độ cao sóng 0,4m, chu kỳ sóng trung bình 2,1s) với hướng lan truyền sóng đạt độ cao sóng lớn 211,5o khu vực có độ sâu khoảng 3,8m xung quanh vị trí (109,298725oE; 12,226684oN) (Hình 10 - 12) 58 T V Chung, N M Tiến, Võ Văn Quang / Tạp chí Khoa học Công nghệ Đại học Duy Tân 06(43) (2020) 52-59 Hình 11 Chu kỳ trung bình (s) tương tác tổng hợp ảnh hưởng mùa gió Tây Nam 3.3 So sánh kết tính tốn số liệu thực đo ngồi trường Để kiểm tra tính đắn mơ hình FEM, cơng trình Trần Văn Chung Bùi Hồng Long (2014) [28] thực so sánh với số liệu thực đo dòng chảy vùng nghiên cứu Trong cơng trình này, để kiểm chứng hiệu chỉnh mơ hình thơng qua mối tương tác sóng - dịng, kết khảo sát Độ cao sóng có nghĩa Trạm Đo đạc (m) Tính tốn (m) Sai số tương đối (%) Lt 0,47 0,67 0,42 0,2 0,49 0,57 0,50 0,20 4,7 14,6 18,1 2,6 Hình 12 Hướng lan truyền sóng (0) tương tác tổng hợp ảnh hưởng mùa gió Tây Nam đồng thời thực máy đo sóng dòng 20 trạm mặt rộng (từ ngày 7-9/9/2018) 01 trạm liên tục (11-12/09/2018) vịnh Nha Trang, thể Hình So sánh kết nghiên cứu mơ hình với trạm đo kiểm tra tính khả thi mơ hình ứng dụng vịnh Nha Trang (Bảng 1) Bảng 1: So sánh số kết tính tốn đo đạc trường trạm đo sóng điển hình Chu kỳ trung bình sóng Đo đạc Tính tốn Sai số tương đối (s) (s) (%) 4,81 4,09 15,1 1,73 2,20 27,2 5,06 4,08 19,4 3,62 3,24 10,3 Nhận xét thảo luận Tác động sóng đến vịnh Nha Trang thường khơng lớn Để đánh giá tác động sóng tổng hợp, nghiên cứu ảnh hưởng tác động trường gió, dịng chảy tầng mặt lên tác động sóng nhằm đánh giá tác động sóng lên khu vực ven bờ vịnh Nha Trang Theo kết nghiên cứu, khu vực chịu tác động mùa gió Đơng Bắc, độ cao sóng có nghĩa trung bình đạt 1,0m, chu kỳ trung bình 3,7s Trong mùa gió Tây Nam ảnh hưởng sóng khơng đáng kể, độ cao sóng có Hướng lan truyền sóng Đo đạc (o) Tính tốn (o) Độ lệch (o) 275,42 124,69 0,63 84,13 272,10 108,84 341,80 82,22 3,3 15,9 18,8 1,9 nghĩa trung bình độ cao sóng 0,4m, chu kỳ sóng trung bình 2,1s So sánh kết tính tốn với thực tế, sai số tương đối trung bình khoảng 10% (thấp 2,6% lớn 18,1%) cho độ cao sóng có nghĩa trung bình khoảng 18% (thấp 10,3% lớn 27,2%) cho chu kỳ sóng trung bình Độ lệch hướng sóng tính tốn đo đạc trung bình khoảng 10o (thấp 1,9o lớn 18,8o) cho hướng lan truyền sóng Các ứng dụng mơ hình sóng SWAN thể khả mơ trường sóng vùng T V Chung, N M Tiến, Võ Văn Quang / Tạp chí Khoa học Công nghệ Đại học Duy Tân 06(43) (2020) 52-59 59 biển thềm lục địa, đầm cửa sông ven biển mà điển hình chứng minh cơng trình với vùng biển với nhiều đảo lớn nhỏ, biên mở phức tạp vịnh Nha Trang cách hiệu ổn định Đồng thời ứng dụng đủ linh hoạt phép liên kết chặt chẽ hệ thống mạng lưới với mơ hình dịng chảy FEM (mơ hình dịng chảy ba chiều phi tuyến theo phương pháp phần tử hữu hạn) Điều hữu ích giải tốn tương tác sóng dịng sóng - sóng Những kết nghiên cứu cho thấy, cách tiếp cận ứng dụng mơ hình tính sóng hợp lý, ổn định cho bước thời gian cho mạng lưới làm mịn, mang đặc trưng địa phương vùng nghiên cứu Lời cảm ơn: Bài báo sử dụng nguồn tài liệu từ đề tài mã số VAST06.03/18-19 “Nghiên cứu chế lưu giữ lại trình phát tán nguồn giống cá (trứng cá - cá bột) vịnh Nha Trang” đề tài với tỉnh Khánh Hòa: “Xác định khu vực có khả cải tạo, phát triển bãi tắm nhân tạo đề xuất phương án bảo vệ bãi tắm tự nhiên vịnh Nha Trang” (2015 - 2016) Các tác giả xin gởi lời cảm ơn chân thành đến Ban chủ nhiệm đề tài đồng nghiệp nhóm nghiên cứu góp ý hỗ trợ chúng tơi hồn thành báo Tài liệu tham khảo [1] Rodriguez, A., Sanchez-Arcilla, A., Redondo, J.M., Bahia, E., Sierra, J.P 1995 Pollutant dispersion in the nearshore region: modelling and measurements Water Sci Technol., 32, pp 169-178 [2] Li, Z., Johns B., 1998 A three-dimensional numerical model of surface waves in the surf zone and longshore current generation over a plane beach Estuar Coast Shelf Sci., 47, pp 395-413 [3] Bever, A.J., MacWilliams, M.L., 2013 Simulating sediment transport processes in San Pablo Bay using coupled hydrodynamic, wave, and sediment transport models Mar Geol., 345, pp 235-253 [4] Olabarrieta, M., Warner, J.C., Kumar, N., 2011 Wave-current interaction in willapa bay J Geophys Res., 116 [5] Roland, A., Zhang, Y.J., Wang, H.V., Meng, Y., Teng, Y.-C., Maderich, V., Brovchenko, I., DutourSikiric, M., Zanke, U., 2012 A fully coupled 3D wave-current interaction model on unstructured grids J Geophys Res., 117 [6] Benetazzo, A., Carniel, S., Sclavo, M., Bergamasco, A., 2013 Wave–current interaction: effect on the wave field in a semi-enclosed basin Ocean Model., 70, pp 152-165 [7] Longuet-Higgins, M.S., 1970 Longshore currents generated by obliquely incident sea waves J Geophys Res., 75 (33), pp 6778-6789 [8] Garcez-Faria, A.F., Thornton, E.B., Lippmann, T.C., Stanton, T.P., 2000 Undertow over a barred beach J Geophys Res., 105 (C7), pp 16999-17010 [9] Longuet-Higgins, M.S., Stewart, R.W., 1964 Radiation stresses in water waves; a physical discussion, with applications Deep Sea Res Oceanogr Abstr., 11 (4), pp 529-562 [10] Guza, R.T., Thornton, E.B., 1981 Wave set-up on a natural beach J Geophys Res., 86 (C5), pp 4133-4137 [11] Dutour-Sikiric, M., Roland, A., Janekovic, I., Tomazic, I., Kuzmic, M., 2013 Coupling of the Regional Ocean modeling system (ROMS) and wind wave model Ocean Model., 72, pp 59-73 [12] Allard, R., Rogers, E., Martin, P., Jensen, T., Chu, P., Campbell, T., Dykes, J., Smith, T., Choi, J., Gravois, U., 2014 The US navy coupled oceanwave prediction system Ocean, 27 (3), pp 93-103 [13] Trần Văn Chung, Nguyễn Trương Thanh Hội, Trần Văn Bình, Tơ Duy Thái Ngơ Mạnh Tiến, 2015 Tính tốn thử nghiệm đặc trưng sóng khu vực Đầm Bấy (vịnh Nha Trang) mơ hình SWAN lưới phi cấu trúc Tuyển tập nghiên cứu biển, tập 21, số Chuyên đề vịnh Nha Trang, tỉnh Khánh Hòa, 13-20 [14] Trần Văn Chung, Nguyễn Hữu Hn, Nguyễn Trương Thanh Hội, 2016 Tính tốn đặc trưng sóng khu vực Ninh Thuận – Bình Thuận mơ hình SWAN lưới phi cấu trúc Tạp chí Khoa học Cơng nghệ biển, Hà Nội 16(2), ISSN 1859 3097, 107 - 114 [15] Booij, N., Ris, R.C., Holthuijsen, L.H., 1999 A third-generation wave model for coastal regions Part I: Model description and validation J Geophys Res., 104 (C4) (1999), pp 7649–7666 [16] Ris, R C., Holthuijsen, L H., and Booij, N., 1999 A third‐generation wave model for coastal regions: Verification Journal of Geophysical Research: Oceans, 104(C4): 7667-7681 [17] Holthuijsen, L H., 2010 Waves in oceanic and coastal waters Cambridge University Press 404 p [18] Bui Hong Long, Tran Van Chung, 2009 (online 2010) Calculations of tidal currents in Vanphong bay using the finite element method Advances in Natural Science, Vol 10, No 4, 495-478 [19] Trần Văn Chung, Bùi Hồng Long, 2015 Một số kết tính tốn dịng chảy Vịnh Bắc Bộ mơ hình ba chiều phi tuyến Tạp chí Khoa học Công nghệ biển, Hà Nội, 15(4), ISSN 1859 - 3097, 320 - 333 [20] Trần Văn Chung, Bùi Hồng Long, 2014 Đặc trưng thủy động lực vực nước Bình Cang - Nha Trang qua mơ hình FEM ECOSMO Tạp chí Khoa học Cơng nghệ biển, Hà Nội, 14(4), ISSN 1859-3097, 320 - 331  ... tả mơ hình 2.1.1 Mơ hình SWAN Các mơ hình sóng sử dụng nghiên cứu mơ hình sóng hệ thứ ba SWAN 41.31 (Mơ sóng gần bờ (SWAN Simulating Waves Nearshore); mơ hình phổ sóng SWAN tính tốn phát triển. .. sóng tác động đến vịnh Nha Trang không đáng kể Kết mô thể chi tiết độ cao lan truyền sóng có nghĩa Hình 10, chu kỳ sóng Hình 11 hướng lan truyền Hình 12 Nói chung, vịnh Nha Trang vịnh kín nên yếu... nghiên Hình Mạng lưới nghiên cứu tác động sóng, dòng cho vịnh Nha Trang cứu, vịnh Nha Trang chịu tác động hướng gió Đơng Bắc, Bắc Đông Bắc, Đông Đông Bắc Bắc Khi trường gió mùa Tây Nam xuất hiện, vịnh