Chúng ta biết rằng bài tốn thủy triều là một bài tốn mơ tả sóng dài nên cần phải kiểm nghiệm mơ hình bài tốn đối với trường hợp sóng đứng. Đây là một trong những bài toán kiểm nghiệm đơn giản nhưng lại rất có hiệu ích trong việc kiểm tra tính bảo tồn năng lượng của sơ đồ số. Bài tốn đặt ra là khi có một sóng dài truyền vào trong kênh thẳng có độ sâu khơng đổi và kín một đầu. Giả thiết sóng phản xạ hồn tồn và khơng có tác dụng của lực Coriolis thì do kết hợp của sóng tới và sóng phản xạ, trường sóng đứng sẽ xuất hiện trong kênh với các điểm nút và điểm bụng phân bố theo quy luật. Nếu suy giảm số của mơ hình số là nhỏ thì tại các điểm bụng ta phải thu được các dao động có biên độ ít nhất là gấp hai lần biên độ của sóng tới. Tại các điểm nút có các dao động nhỏ. Nếu tính phi tuyến được mơ phỏng tốt thì các dao động tại điểm nút khơng bị triệt tiêu mà sẽ xuất hiện sóng thứ cấp có tần số gấp hai lần tần số của sóng tới.
Hình 3.3: Kênh dùng để thử nghiệm
Với bài tốn kiểm tra thứ hai này, kênh được chọn có độ dài l =1000m, độ rộng
m
b=50 , độ sâu h=10m, sóng tới có chu kỳ T =30s, biên độ A=20cm. Lưới tính trong kênh được chia thành các lưới hình chữ nhật 2.5m×2.5m. Biên tại cuối kênh được chọn là biên không thấm (phản xạ tồn phần). Thời gian tích phân được chọn là 250s đủ để sóng truyền hết hai lần chiều dài kênh (xem hình 3.3). Các kết quả mơ phỏng dao động mực nước tại các điểm bụng và điểm nút (cách nhau xen kẽ bằng một phần tư độ dài sóng L≈T gh) được trình bày trên các hình 3.4 và 3.5:
Hình 3.4: Dao động mực nước tại điểm bụng và dao động sóng tới
Hình vẽ 3.4 trình bày mực nước ngay tại điểm bụng ngay phía trước tường đứng và dao động của sóng tới tại biên vào đầu kênh. Trên hình vẽ ta thấy, sau một khoảng thời gian khoảng 100s sóng truyền tới tường và hiện tượng phản xạ xuất hiện. Dao động tại điểm sát tường có biên độ gấp hai lần biên độ sóng tới và ổn định sau một khoảng thời gian rất ngắn. Sau khoảng 200s, sóng phản xạ từ tường kết hợp với sóng tới tại biên tạo thành một sóng kết hợp, lúc này dao động của sóng tại biên khơng cịn là dao động của riêng sóng tới nữa, ta có thể thấy rõ sự biến đổi của đường đứt nét trên hình sau của dao động tại biên.
Hình 3.5: Dao động mực nước tại điểm nút và dao động tại biên
Hình vẽ 3.5 trình bày dao động của mực nước tại điểm cách tường một khoảng 5/4 độ dài sóng tới L, tương ứng với điểm gần điểm nút. Ta thấy sau khi sóng tới đi
qua một khoảng thời gian ngắn cỡ 130s thì sóng phản xạ từ tường quay lại gặp sóng tới tại điểm này và tạo thành sóng kết hợp có dao động nhỏ hơn rất nhiều dao động của sóng tới. Đặc biệt ta có thể quan sát thấy dao động có tần số cao tại điểm này do sự kết hợp phi tuyến giữa sóng tới và sóng phản xạ chính vì vậy các dao động mực nước vấn tồn tại ở điểm nút.
3.4 Kiểm nghiệm với bài tốn thủy triều trong Biển Đơng
Mơ hình số đã kiểm nghiệm ở trên được áp dụng tính tốn thử cho truyền triều trong Biển Đơng. Vùng biển được chia thành các lưới hình vng với ∆x=1/5 độ
5 / 1
=
∆y độ, các trục được chọn xấp xỉ hướng theo kinh tuyến và vĩ tuyến trên hệ tọa
độ Đề Các. Biên cứng được xấp xỉ với điều kiện thông lượng bằng không. Các biên lỏng được cho các dao động mực nước theo các hằng số điều hịa của 4 sóng. Có ba biên lỏng được sử dụng cho mơ hình là eo Đài Loan, eo Basi và eo Malaca. Các hằng số điều hòa trên các eo biển được lấy theo bảng thủy triều Admiralty Tidal Tables.
Hình 3.6: Dao động mực nước tại khu vực Hòn Dấu (bắt đầu từ 22h ngày 4/4/2004)
Kết quả mơ phỏng dao động mực nước tại Hịn Dấu được so sánh với số liệu thực đo tại trạm Hịn Dấu trong tháng 4 năm 2004. Hình 3.6 trình bày dao động mực nước triều tại khu vực Hòn Dấu trong tháng 4 năm 2004. Trên hình vẽ ta thấy dao động mực nước tính tốn khá phù hợp với số liệu thực đo. Dao động mực nước thể hiện rõ tính nhật triều thuần nhất và xu thế thay đổi của dao động mực nước tính tốn rất phù hợp với thực tế. Trên hình vẽ so sánh ta cũng thấy có sự sai khác nhỏ về biên độ triều giữa mô phỏng số và số liệu thực đo, ngun nhân có thể giải thích do vị trí lấy kết quả mơ phỏng sai lệch khá lớn so với vị trí có số liệu vì lưới là 1/5 độ và hơn nữa do lưới mơ phỏng cịn khá thô nên chưa mô tả hết được những ảnh hưởng phi tuyến của vùng nước nơng và sự thay đổi địa hình phức tạp lên dao động mực nước triều.
Hình 3.7: Dao động mực nước tại khu vực Hịn Ngư (bắt đầu từ 00h ngày 7/1/2001)
Kết quả mô phỏng dao động mực nước tại Hòn Ngư được so sánh với số liệu thực đo tại trạm Hịn Ngư trong tháng 1 năm 2001. Hình 3.7 trình bày dao động mực nước triều tại khu vực Hòn Ngư trong tháng 1 năm 2001. Trên hình vẽ ta thấy dao động mực nước tính tốn khi triều lên khá phù hợp với số liệu thực đo. Trên hình vẽ so sánh ta
cũng thấy có sự sai khác về biên độ triều giữa mô phỏng số và số liệu thực đo, ngun nhân có thể giải thích do chưa hiệu chỉnh đúng được tham số ma sát và ảnh hưởng do địa hình chưa đúng. Tuy nhiên triều lên mơ phỏng khá đúng, pha triều cũng khá chính xác.
Hình 3.8: Dao động mực nước tại khu vực Vũng Tàu (bắt đầu từ 22h ngày 7/1/2000)
Kết quả mô phỏng dao động mực nước tại Vũng Tàu được so sánh với số liệu thực đo tại trạm Vũng Tàu trong tháng 1 năm 2000. Hình 3.8 trình bày dao động mực nước triều tại khu vực Vũng Tàu trong tháng 1 năm 2000. Trên hình vẽ ta thấy dao động mực nước tính tốn khá phù hợp với số liệu thực đo. Dao động mực nước thể hiện rõ tính bán nhật triều không đều và xu thế thay đổi của dao động mực nước tính tốn rất phù hợp với thực tế. Trên hình vẽ so sánh ta cũng thấy có sự sai khác nhỏ về biên độ triều giữa mô phỏng số và số liệu thực đo, nguyên nhân có thể giải thích do những yếu tố đã nêu ở phần trên.
Hình 3.9: Dao động mực nước tại khu vực Qui Nhơn (bắt đầu từ 22h ngày 23/5/2002)
Kết quả mô phỏng dao động mực nước tại Qui Nhơn được so sánh với số liệu thực đo tại trạm Qui Nhơn trong tháng 5 năm 2002. Hình 3.9 trình bày dao động mực nước triều tại khu vực Qui Nhơn trong tháng 5 năm 2002. Trên hình vẽ ta thấy dao động mực nước tính tốn khá phù hợp với số liệu thực đo. Dao động mực nước thể hiện rõ tính nhật triều khơng đều và xu thế thay đổi của dao động mực nước tính tốn rất phù hợp với thực tế. Trên hình vẽ so sánh ta cũng thấy có sự sai khác nhỏ về biên độ triều giữa mô phỏng số và số liệu thực đo.
3.5 Kết luận tóm lược
Kết quả của việc kiểm chứng cho ba bài tốn nói trên, bài tốn vỡ đập và bài toán truyền triều trong kênh và bài tốn truyền triều Biển Đơng cho thấy: Mơ hình số có thể mơ phỏng tốt trường hợp lan truyền sóng dài và mơ phỏng được các tương tác giữa các sóng với sự suy giảm số khơng đáng kể. Do đó ta có cơ sở tin tưởng vào việc sử dụng phương pháp CIP cho bài tốn sóng dài thực tế bao gồm cả bài toán truyền triều và bài toán nước dâng do bão.
Chương 4 - TÍNH TỐN NƯỚC DÂNG DO BÃO
4.1 Khái quát vùng nghiên cứu4.1.1 Vị trí địa lí 4.1.1 Vị trí địa lí
Vịnh Bắc Bộ nằm ở phía tây bắc Biển Đơng, là một vịnh nửa kín nước nơng. Ba mặt vịnh được bao bọc bởi lục địa Việt Nam ở phía tây, Trung Quốc ở phía bắc và đảo Hải Nam ở phía đơng. Vịnh ăn thơng ra Biển Đơng qua hai cửa, cửa chính nằm giữa bờ phía nam đảo Hải Nam và bờ biển Trung Bộ, Việt Nam và cửa thứ hai qua eo biển Quỳnh Châu nằm giữa bán đảo Lôi Châu, Trung Quốc và mũi Oanh Ca, đảo Hải Nam. Diện tích vịnh Bắc Bộ vào khoảng 126.250 km2; Chiều dài khoảng trên 500 km; chiều rộng nơi rộng nhất khoảng 310 km, nơi hẹp nhất khoảng 220 km. Hình 4.1 trình bày sơ lược độ sâu vùng nghiên cứu.
4.1.2 Lưới tính
Vùng nghiên cứu được chia thành các lưới tính tốn hình chữ nhật, có kích thước 30
/ 1
= =dy
dx độ. Trục oy theo chiều dọc theo kinh tuyến được chia thành 181 điểm tính (từ vĩ độ 16 đến vĩ độ 22), trục ox theo chiều ngang song song với vĩ tuyến được chia thành 151 điểm tính (từ kinh độ 105 đến kinh độ 110). Tọa độ gốc có kinh độ
0 0 =104.97
x , vĩ độ 0
0 =15.97
y .
4.1.3 Đặc điểm nước dâng do bão vùng Vịnh Bắc Bộ
Dọc theo ven biển Việt Nam thì nước dâng do bão có xu hướng giảm dần từ bắc vào nam. Mực nước dâng cao nhất có thể xẩy ra là lớn nhất chỉ xuất hiện ở vùng Vịnh Bắc Bộ. Hơn nữa đây là vùng mà có số liệu đo nước dâng do bão cụ thể. Trong lịch sử đã đo được mực nước dâng lớn nhất là 3.5m.
4.2 Số liệu tính tốn mơ hình
Trong những năm gần đây, một số cơ quan bộ ngành đã tiến hành đo đạc khảo sát nước dâng do bão khi bão xẩy ra và sau khi đổ bộ. Dựa vào những kết quả khảo sát số liệu thực tế của hai cơn bão số 6 và 7 năm 2005. Chúng tôi lựa chọn hai cơn bão này để tính nước dâng do bão. Mục đích là nhằm có số liệu để so sánh. Đồng thời cũng sử dụng để so sánh với kết quả của mơ hình Delft 3D của Viện thủy lực Hà Lan.
Cơn bão số 6
Bão số 6 có tên quốc tế là Vicente hình thành tại phía nam quần đảo Hồng Sa từ một áp thấp nhiệt đới. Bão số 6 có quỹ đạo di chuyển rất phức tạp, lúc đầu tưởng chừng bão đổ bộ vào các tỉnh Nam Bộ nhưng sau đó bão đi theo đường vịng cung và chuyển hướng lên phía bắc, lúc vào gần đến bờ tưởng chừng bão đổ bộ vào giữa tỉnh Quảng Bình và Hà Tĩnh nhưng sau đó lại đi dọc bờ trước khi đổ bộ vào thành phố Vinh - Nghệ An. Bão số 6 là cơn bão mạnh có tốc độ di chuyển tương đối chậm, độ giảm áp ở tâm tại thời điểm bão gần bờ có lúc đến 985 mlb (bảng 4.1), tốc độ gió giật 25 m/s.
Bảng 4.1: Các tham số cơn bão Vicente 2005 (cơn bão số 6). Ngày/giờ Kinh độ (độ) Vĩ độ (độ) Rmax (km) Vmax (m/s) Pmin (mlb) P∞ (mlb) 15/09/05 00:18 113.6 12.9 35 30 1004 1013 16/09/05 00:00 112.7 12.2 35 30 1004 1013 16/09/05 00:06 112.9 11.1 35 30 1000 1013 16/09/05 00:12 113.2 10.7 35 30 992 1013 16/09/05 00:18 114.1 12.1 35 30 990 1013 17/09/05 00:00 114.4 13.6 35 35 990 1013 17/09/05 00:06 113.8 15.2 35 35 990 1013 17/09/05 00:12 111.7 16.4 35 35 990 1013 17/09/05 00:18 110.2 16.8 35 35 990 1013 18/09/05 00:00 108.3 16.8 35 40 985 1013 18/09/05 00:06 106.9 17.4 35 40 992 1013 18/09/05 00:12 105.1 18.7 35 30 994 1013 18/09/05 00:18 102.4 19.3 35 30 998 1013 19/09/05 00:00 100 20 35 30 1004 1013 Cơn bão số 7
Cơn bão số 7 có tên quốc tế là Damrey được coi là mạnh nhất trong vòng 9 năm qua. Bão số 7 được hình thành tại vùng biển phía bắc của Philippin từ một áp thấp nhiệt đới. Bão số 7 có quỹ đạo di chuyển tương đối giống bão số 2, tức là tây và tây - tây bắc, một hướng di chuyển rất thông thường của các cơn bão đổ bộ vào Việt Nam. Bão số 7 được những trung tâm dự báo bão trên thế giới có những dự báo xa rất chính xác cả về quỹ đạo di chuyển và cường độ bão. Ngay sau khi hình thành bão số 7 đã di chuyển nhanh và cường độ tăng mạnh, có thời điểm khi đổ vào Vịnh Bắc Bộ, độ giảm áp ở tâm tụt xuống tới 955 mlb, tốc độ gió giật 42 m/s, lúc vào gần bờ Việt Nam bão vẫn duy trì cường độ rất mạnh trước khi đổ bộ với áp ở tâm là 975 mlb, tốc độ gió giật tới 30 m/s (xem bảng 4.3). Mặc dù vậy những thiệt hại về đê biển do bão số 7 gây ra là rất lớn, đặc biệt là tại tỉnh Nam Định.
Bảng 4.2: Các tham số cơn bão Damrey 2005 (cơn bão số 7). Ngày/giờ Kinh độ(độ) Vĩ độ(độ) Rmax(km) Vmax(m/s) (mlb)Pmin P∞
(mlb) 24/09/05 19:00 114 19.6 35 30.42 975 1013 24/09/05 22:00 113.7 19.5 35 1013 25/09/05 01:00 113.4 19 35 1013 25/09/05 04:00 113.2 18.9 35 1013 25/09/05 07:00 113 18.9 35 25.35 980 1013 25/09/05 10:00 112.6 18.9 35 38.02 960 1013 25/09/05 13:00 112.3 18.9 35 1013 25/09/05 16:00 112.1 18.9 35 48.16 960 1013 25/09/05 19:00 111.8 19 35 40.56 955 1013 25/09/05 22:00 111.3 19.1 35 40.56 955 1013 26/09/05 01:00 110.7 19.1 35 40.56 955 1013 26/09/05 04:00 110.3 19 35 38.02 955 1013 26/09/05 07:00 109.9 18.8 35 38.02 955 1013 26/09/05 10:00 109.5 18.9 35 38.02 955 1013 26/09/05 13:00 109 18.9 35 38.02 955 1013 26/09/05 16:00 108.6 19 35 38.02 955 1013 26/09/05 19:00 108.1 19.2 35 30.42 975 1013 26/09/05 22:00 107.8 19.3 35 1013 27/09/05 01:00 107.5 19.4 35 30.42 975 1013 27/09/05 04:00 107.1 19.5 35 30.42 975 1013 27/09/05 07:00 106.3 19.6 35 30.42 975 1013 27/09/05 10:00 105.6 19.7 35 27.88 980 1013 27/09/05 13:00 105 19.7 35 25.35 985 1013 27/09/05 16:00 35 1013 27/09/05 19:00 103 19.5 35 17.74 996 1013 4.3 Các kết quả tính tốn
Q trình tính tốn trường gió và áp suất được bắt đầu bằng việc chuẩn bị số liệu đầu vào, bao gồm thời gian, vị trí tâm bão, độ giảm áp tại tâm, tốc độ gió cực đại, bán kính gió cực đại. Khoảng cách giữa hai lần số liệu đầu vào là 6 giờ.
Trường gió và áp suất trong hai cơn bão được mô phỏng bằng mơ hình bão của Isozaki (trình bày trong chương 2). Các kết quả mơ phỏng gió và áp cho hai cơn bão được trình bày trên các hình 4.2 đến hình 4.11.
Hình 4.2: Trường véc tơ vận tốc gió và áp lúc 18h ngày 17/9/2005 (cơn bão số 6 năm 2005 - Vicente)
Hình 4.3: Trường véc tơ vận tốc gió và áp lúc 00h ngày 18/9/2005 (cơn bão số 6 năm 2005 - Vicente)
Hình 4.4: Trường véc tơ vận tốc gió và áp lúc 06h ngày 18/9/2005 (cơn bão số 6 năm 2005 - Vicente)
Hình 4.5: Trường véc tơ vận tốc gió và áp lúc 12h ngày 18/9/2005 (cơn bão số 6 năm 2005 - Vicente)
Hình 4.6: Trường véc tơ vận tốc gió và áp lúc 00h ngày 26/9/2005 (cơn bão số 7 năm 2005 - Damrey)
Hình 4.7: Trường véc tơ vận tốc gió và áp lúc 06h ngày 26/9/2005 (cơn bão số 7 năm 2005 - Damrey)
Hình 4.8: Trường véc tơ vận tốc gió và áp lúc 12h ngày 26/9/2005 (cơn bão số 7 năm 2005 - Damrey)
Hình 4.9: Trường véc tơ vận tốc gió và áp lúc 18h ngày 26/9/2005 (cơn bão số 7 năm 2005 - Damrey)
Hình 4.10: Trường véc tơ vận tốc gió và áp lúc 00h ngày 27/9/2005 (cơn bão số 7 năm 2005 - Damrey)
Hình 4.11: Trường véc tơ vận tốc gió và áp lúc 06h ngày 27/9/2005 (cơn bão số 7 năm 2005 - Damrey)
Từ hình 4.2 đến hình 4.5 là kết quả tính phân bố trường gió và áp suất trong cơn