Luận án tiến sĩ HUS nghiên cứu xây dựng mô hình số trị ba chiều cho vùng biển nước nông ven bờluận án TS hải dương học1 07 07

Tác giả luận án xin bày tỏ lòng biết ơn tới PGS TS Đinh Văn Ƣu, TS Bùi Xuân Thông, tận tình hƣớng dẫn tơi thực hồn thành luận án; Cảm ơn chuyên gia ngồi trƣờng đọc thảo đóng góp ý kiến quý báu cho luận án Luận án nhận đƣợc hỗ trợ tích cực đề tài cấp Nhà nƣớc : ‘Luận chứng khoa học mơ hình phát triển kinh tế – sinh thái số đảo, cụm đảo lựa chọn vùng biển ven bờ Việt Nam‘, mã số KC-09-12 thuộc Chƣơng trình Biển GS TS Lê Đức Tố làm chủ nhiệm Những động viên, giúp đỡ giúp chúng tơi vƣợt qua khó khăn để hồn chỉnh luận án NGUYỄN MINH HUẤN MỤC LỤC Trang bìa phụ Lời cam đoan Mục lục Danh mục ký hiệu, chữ viết tắt Danh mục bảng Danh mục hình vẽ, đồ thị MỞ ĐẦU LUAN VAN CHAT LUONG download : add luanvanchat@agmail.com CHƢƠNG 1: TỔNG QUAN TÌNH HÌNH NGHIÊN CỨU DỊNG CHẢY BA CHIỀU VÙNG BIỂN VEN BỜ 1.1 Các phƣơng pháp nghiên cứu dòng chảy ba chiều 1.2 Tình hình nghiên cứu giới nƣớc 20 1.3 Các vấn đề nghiên cứu luận án 23 CHƢƠNG 2: CƠ SỞ TỐN HỌC CỦA MƠ HÌNH THUỶ ĐỘNG LỰC BA CHIỀU 26 2.1 Hệ phƣơng trình thủy động lực ba chiều hệ toạ độ sigma 26 2.2 Các phƣơng pháp tham số hố q trình trao đổi rối 31 2.3 Các điều kiện biên điều kiện đầu 45 2.4 Tác động sóng trọng lực bề mặt 57 CHƢƠNG 3: PHƢƠNG PHÁP SỐ GIẢI HỆ PHƢƠNG TRÌNH THỦY ĐỘNG LỰC 63 3.1 Sơ đồ lƣới tính số ký hiệu 63 3.2 Thuật giải hệ phƣơng trình bảo toàn động lƣợng 67 3.3 Thuật giải hệ phƣơng trình biến vơ hƣớng 97 3.4 Thuật giải hệ phƣơng trình đặc trƣng rối 109 3.5 Các bƣớc tính tốn sơ đồ khối tổng quát 116 CHƢƠNG 4: CÁC KẾT QUẢ THỬ NGHIỆM MƠ HÌNH THỦY ĐỘNG LỰC BA CHIỀU 116 4.1 Kết thử nghiệm so sánh sơ đồ bình lƣu 116 4.2 Kết thử nghiệm so sánh sơ đồ khép rối 126 4.3 Kết thử nghiệm tác động sóng trọng lực 133 bề mặt CHƢƠNG 5: KẾT QUẢ ỨNG DỤNG MƠ HÌNH TÍNH TỐN HỒN LƢU VÙNG BIỂN QUẢNG NAM KHU VỰC HỘI AN - CÙ LAO CHÀM LUAN VAN CHAT LUONG download : add luanvanchat@agmail.com 5.1 Chế độ khí tƣợng thủy văn vùng biển Quảng Nam khu vực 136 Hội An – Cù Lao Chàm 5.2 Chế độ hoàn lƣu vùng biển nghiên cứu qua kết tính 137 tốn mơ hình KÕt luận 146 Danh mục công trình tác giả 169 Tài liệu tham khảo 172 173 LUAN VAN CHAT LUONG download : add luanvanchat@agmail.com DANH MỤC CÁC KÝ HIỆU, CÁC CHỮ VIẾT TẮT Ý nghĩa Ký hiệu a1 – a15 aw Ahx Ahy Ah Av Các hệ số phương trình trạng thái Biên độ sóng trọng lực Thành phần bình lưu theo phương ngang hướng X Thành phần bình lưu theo phương ngang hướng Y Thành phần bình lưu tổng cộng theo phương ngang Thành phần đối lưu theo phương thẳng đứng Ahx Thành phần bình lưu 2D hướng X Ahy Thành phần bình lưu 2D hướng Y Ah Thành phần bình lưu 2D tổng cộng b c Cg cp Độ Vận tốc sóng trọng lực Vận tốc nhóm sóng Nhiệt dung riêng đẳng áp nước biển CDs Hệ số ma sát gió CDb Hệ số ma sát đáy CE CH Dhx Dhy Dh Dxx Dyx Dxy Dyy Dv Hệ số Dalton Hệ số Stanton Thành phần khuếch tán theo phương ngang hướng X Thành phần khuếch tán theo phương ngang hướng Y Thành phần khuếch tán tổng cộng theo phương ngang Thành phần khuếch tán u hướng X Thành phần khuếch tán u hướng Y Thành phần khuếch tán v hướng X Thành phần khuếch tán v hướng Y Thành phần khuếch tán tổng cộng (u, v biến vô hướng) Dxx Thành phần khuếch tán U hướng X Dyx Thành phần khuếch tán U hướng Y Dxy Thành phần khuếch tán V hướng X Dyy Thành phần khuếch tán V hướng Y f Tham số Coriolis i LUAN VAN CHAT LUONG download : add luanvanchat@agmail.com F Fhar Hàm toạ độ sigma Hàm điều hoà g h H I J Gia tốc trọng trường Độ sâu bề mặt đáy so với mực nước trung bình Độ sâu tổng cộng Bức xạ mặt trời mơi trường nước biển Tốn tử chuyển đổi hệ toạ độ Đề Các sigma k kw l L Năng lượng rối Số sóng Độ dài xáo trộn Độ dài tỉ lệ Mt = t3D/t2D Số bước thời gian 2D (thành phần áp) bước tính 3D (thành phần tà áp) Mtot Mx Số bước thời gian 2D tổng cộng Đại lượng hiệu chỉnh hướng X hệ phương trình truyền tải biến vơ hướng My N Đại lượng hiệu chỉnh hướng Y hệ phương trình truyền tải biến vơ hướng Đại lượng hiệu chỉnh hướng thẳng đứng hệ phương trình truyền tải biến vơ hướng Tần số Brunvaisala Nx Ny Nz Số lượng nút lưới hướng X Số lượng nút lưới hướng Y Số lượng nút lưới hướng Z Nt Số bước thời gian 3D p Áp suất p Thành phần thăng giáng mạch động áp suất Pa Qi Áp suất khơng khí Thành phần gradien áp suất tà áp baroclin R Biến đặc trưng Riemann tới khỏi biên Ri S Sm, Sh, Su, Sb T Số Richardson Độ muối Hàm ổn định sơ đồ khép kín rối Nhiệt độ nước Ta TVD Nhiệt độ khơng khí Thuật giải suy giảm biến động tổng cộng Mz ii LUAN VAN CHAT LUONG download : add luanvanchat@agmail.com uF Vận tốc dòng chảy theo phương X tham gia vào q trình bình lưu biến vơ hướng U, V, W u, v, w u*s, u*b Uw Các thành phần vận tốc tức thời Các thành phần vận tốc trung bình theo thời gian Vận tốc ma sát bề mặt đáy Vận tốc quỹ đạo sóng u Thành phần thăng giáng mạch động thành phần vận tốc u UF Giá trị dịng chảy tích phân U theo độ sâu lấy trung bình bước thời gian tính 3D Vận tốc dịng chảy theo phương Y tham gia vào trình bình lưu biến vô hướng vF v Thành phần thăng giáng mạch động thành phần vận tốc v VF x, y , z Giá trị dịng chảy tích phân V theo độ sâu lấy trung bình bước thời gian tính 3D Các trục toạ độ Đề Các x1, x2, x3 Các trục toạ độ Đề Các w ~ w Thành phần thăng giáng mạch động thành phần vận tốc w zr Cao độ từ đáy biển tới trung tâm ô lưới sát đáy x1 Độ lớn bước lưới tính hướng X trung tâm x2 Độ lớn bước lưới tính hướng Y trung tâm ô x3 Độ lớn bước lưới tính hướng thẳng đứng trung tâm ô k Độ lớn bước lưới hướng thẳng đứng hệ toạ Thành phần vận tốc theo phương thẳng đứng hệ toạ độ sigma độ sigma trung tâm t2D Bước thời gian tính thành phần 2D t3D Bước thời gian tính thành phần 3D x Toán tử sai phân hướng X y Toán tử sai phân hướng Y z Toán tử sai phân hướng thẳng đứng b Độ dày lớp biên đáy a Hệ số xác định tính chất sơ đồ sai phân thành phần đối lưu thẳng đứng, a = hiện; a = ẩn v Hệ số xác định tính chất sơ đồ sai phân thành phần khuếch tán thẳng đứng, v = hiện; v = ẩn iii LUAN VAN CHAT LUONG download : add luanvanchat@agmail.com (r) Hàm trọng số sơ đồ sai phân ngược dòng Lax-Wendroff thành phần thơng lượng bình lưu theo phương ngang sơ đồ sai phân ngược dòng trung tâm thành phần thông lượng đối lưu theo phương thẳng đứng  Tần số quay trái đất T Hệ số nhớt rối theo phương thẳng đứng H Hệ số khuếch tán động lượng theo phương ngang T Hệ số khuếch tán rối theo phương thẳng đứng  Mật độ nước 11, 21, 12, 22 Các tenxơ ứng suất thành phần theo phương ngang s1, s2 Tenxơ ứng suất bề mặt thoáng theo phương ngang b1, b2 Tenxơ ứng suất bề mặt đáy theo phương ngang c Ứng suất biến dạng tiếp tuyến đáy dòng chảy gây w,max Ứng suất biến dạng tiếp tuyến đáy cực đại sóng gây b,max Ứng suất biến dạng tiếp tuyến đáy cực đại tổng cộng sóng dịng chảy gây  Toạ độ sigma  Hệ số sóng trọng lực  Tần số sóng trọng lực n Tần số phân triều  Thế vận tốc sóng trọng lực n0 Pha ban đầu phân triều n Góc vị đặc trưng phân triều  Cao độ mặt thoáng so với mực nước trung bình n Biên độ phân triều 0p, 0m Các giá trị ban đầu hệ số nhớt rối  Cường độ tiêu tán lượng rối T Hệ số dãn nở nhiệt S Hệ số dãn nở muối DANH MỤC CÁC HÌNH VẼ Hình 2.1 Hệ toạ độ chuyển đổi sigma theo chiều thẳng đứng iv LUAN VAN CHAT LUONG download : add luanvanchat@agmail.com Hình 2.2 Biến động hàm ổn định động lượng biến vô hướng theo số Richardson Ri Hình 2.3 Hình 3.1 Hình 3.2 Hình 3.3 Kết so sánh mơ hình dự báo ứng suất biến dạng tiếp tuyến đáy trung bình cực đại tác động sóng dịng chảy Sơ đồ lưới tính Sơ đồ xác định thành phần thơng lượng bình lưu Sơ đồ khối tổng qt mơ hình dịng chảy chiều Hình 4.1 Hình 4.2 Trạng thái ban đầu phân bố nồng độ muối dòng chảy Phân bố nồng độ muối bề mặt chu kỳ 0,5T, T, 1,5T 2T Sơ đồ ngược dịng – trường hợp A Hình 4.3 Phân bố nồng độ muối bề mặt chu kỳ 0,5T, T, 1,5T 2T Sơ đồ Lax – Wendroff – trường hợp B Hình 4.4 Phân bố nồng độ muối bề mặt chu kỳ 0,5T, T, 1,5T 2T Sơ đồ suy giảm biến động tổng cộng TVD với hàm giới hạn – trường hợp C Hình 4.5 Hình 4.6 Hình 4.7 Hình 4.8 Hình 4.9 Vị trí ban đầu fron muối Biến trình phát triển dòng chu kỳ triều trường hợp A Biến trình phát triển dịng chu kỳ triều trường hợp B Biến trình phát triển dòng chu kỳ triều trường hợp C Biến trình phát triển dịng chu kỳ triều trường hợp D chảy phân bố độ muối chảy phân bố độ muối chảy phân bố độ muối chảy phân bố độ muối Trường ứng suất đáy với đặc trưng sóng H = 0m; T = A - Địa hình đáy khơng đổi z = 5m; B - Địa hình đáy biến đổi z = 2x2/3 m Trường ứng suất đáy với đặc trưng sóng H = 0,5m; T = A - Địa hình đáy khơng đổi z = 5m; B - Địa hình đáy biến đổi z = 2x2/3 m 0s Hình 4.12 Trường ứng suất đáy với đặc trưng sóng H = 1,5m; T A - Địa hình đáy khơng đổi z = 5m; B - Địa hình đáy biến đổi z = 2x2/3 m 7s Hình 5.1 Hình 5.2 Hình 5.3 Bản đồ khu vực Quảng Nam – Hội An – Cù Lao Chàm Hoa gió trung bình tháng Biến trình mực nước thực đo đảo Cù Lao thời gian từ ngày 26-4-2002 đến 1-5-2002 Trường sóng gió trung bình tháng Kết so sánh giá trị tính tốn trạm đo dịng chảy liên tục TK1 Hình 4.10 Hình 4.11 Hình 5.4 Hình 5.5a = 7s Chàm [m] thực đo v LUAN VAN CHAT LUONG download : add luanvanchat@agmail.com Kết so sánh giá trị trạm đo dịng chảy liên tục TK2 Hình 5.6 Hình 5.7a Kết so sánh giá trị tính tốn thực đo mực nước Hồn lưu tháng VII đại diện mùa khô theo phương ngang tầng mặt (0,2H) Hồn lưu tháng VII đại diện mùa khơ theo phương ngang tầng (0,5H) Hình 5.7b tính tốn thực đo Hình 5.5b Hình 5.7c Hình 5.8 Hình 5.9a Hồn lưu tháng VII đại diện mùa khơ theo phương ngang tầng đáy (1H) Phân bố độ muối tháng VII đại diện mùa khơ Hồn lưu tháng XII đại diện mùa mưa theo phương ngang tầng mặt (0,2H) Hình 5.9b Hồn lưu tháng XII đại diện mùa mưa theo phương ngang tầng (0,5H) Hoàn lưu tháng XII đại diện mùa mưa theo phương ngang tầng đáy (1H) Hình 5.9c Hình 5.10 Hình 5.11 Hình 5.12 Hình 5.13 Phân bố độ muối tháng XII đại diện mùa mưa Phân bố dịng thăng mùa khơ (tháng VII) Phân bố dòng thăng tháng XII đại diện mùa mưa Phân bố dòng chảy theo mặt cắt: A tháng VII đại diện mùa khô B tháng XII đại diện mùa mưa DANH MỤC CÁC BẢNG Bảng 2.1 Giá trị tham số sử dụng sơ đồ khép kín rối Bảng 2.2 Giá trị trung bình zo dạng đáy khác Bảng 2.3 Chỉ dẫn lựa chọn mơ hình tính tốn lan truyền sóng Bảng 5.1 Tốc độ gió trung bình tháng cực đại (m/s) vi LUAN VAN CHAT LUONG download : add luanvanchat@agmail.com Bảng 5.2 Số lượng bão (áp thấp nhiệt đới) đổ ảnh hưởng trực tiếp đến khu vực nghiên cứu Bảng 5.3 Bảng 5.4 Lưu lượng nước trung bình lớn thời kỳ quan trắc 1977 – 1985 Giá trị biên độ góc vị đặc trưng phân triều trạm Sơn Trà vii LUAN VAN CHAT LUONG download : add luanvanchat@agmail.com [m 130] 30x103 120 110 100 90 803 20x10 70 60 Sơn Trà 50 403 10x10 30 20 Cửa Đại 10 An Hòa 25 50 75 20x103 100 125 15040x103 175 200 225 60x10250 275 300 80x10325 [m ] A 130 [m ] 30x10 120 110 100 90 20x10 80 70 60 Sơn 50 Trà 10x10 40 30 20 10 Cửa Đại 0 25 50 20x10 75 100 125 An Hßa 15040x10 175 B 200 225 60x10250 275 300 80x10325 [m ] 130 [m ] 1203 30x10 110 100 90 803 20x10 70 60 50 Sơn Trà 403 10x10 30 20 10 Cửa Đại 0 25 50 20x103 75 100 125 An Hßa 40x103 150 175 200 225 60x103250 275 300 80x103[m 325] Hình 5.11 Phân bố dòng thăng tháng VII (mùa khô) (A- tầng mặt; B- tầng giữa; C- ®¸y) - 166 - LUAN VAN CHAT LUONG download : add luanvanchat@agmail.com [m ] 130 30x103 120 110 100 90 20x103 80 70 60 Sơn 50 Trà 10x10 40 30 20 10 Cửa Đại 25 50 20x10 75 25 50 75 20x103 100 125 An Hßa 150 40x10 175 200 225 15040x103 175 200 225 60x10250 A 60x10 250 275 300 275 300 80x10325 [m ] 130 [m ] 30x10 120 110 100 90 20x10 80 70 60 50 S¬n Trà 40 10x10 30 20 10 Cửa Đại 0 100 125 An Hßa B 80x10325 [m ] 130 [m ] 30x10 120 110 100 90 20x10 80 70 60 Sơn 50 Trà 10x10 40 30 20 10 Cửa Đại 25 50 20x103 75 100 125 An Hßa 15040x10 175 200 225 60x10250 275 300 80x10325 [m ] H×nh 5.12 Phân bố dòng thăng tháng XII (mùa m-a) (A- tầng mặt; B- tầng giữa; C- đáy) - 167 - LUAN VAN CHAT LUONG download : add luanvanchat@agmail.com - 168 - LUAN VAN CHAT LUONG download : add luanvanchat@agmail.com -2 -10 -20 -30 -40 -50 -60 -70 2.10 4.10 130 v w 100 [m] 50 Hìn thá B A Cử Hò Già La Cù x 250m [m/s] h Tr-ờng ng VII đại a 5.13 in o Đại dòng chảy diện mùa theo mặt khôtháng cắt XII CĐđại CLC Amùa diÖn - ; Bm-a - 169 - LUAN VAN CHAT LUONG download : add luanvanchat@agmail.com KẾT LUẬN Vùng nước ven bờ có chế độ thủy động lực vơ phức tạp, khu vực chịu tác động mạnh mẽ từ ngồi biển thủy triều, dịng chảy nước từ sơng, gió sóng gió Việc mơ chế độ hoàn lưu khu vực này tốn phức tạp có nhiều ý nghĩa thực tiễn Các kết luận án góp phần làm phong phú thêm phương pháp hữu hiệu giải toán thủy động lực ba chiều cho vùng nước nơng ven biển có cửa sơng thử nghiệm áp dụng có kết cho khu vực biển cụ thể Việt Nam Các kết khoa học thực tiễn luận án: Xây dựng chương trình tính có độ tin cậy sử dụng hệ phương trình thủy nhiệt động lực quy mô vừa đầy đủ cho phép giải số trị toán thủy động lực ba chiều hệ toạ độ Đề Các theo phương ngang toạ độ sigma theo phương thẳng đứng Mơ hình xây dựng nhằm mơ trạng thái thực hệ thống hoàn lưu vùng nước ven bờ có cửa sơng, xét đến tác động tổng cộng thủy triều, ứng suất gió bề mặt, gradien mật độ nước biến động độ muối ảnh hưởng sóng trọng lực bề mặt biển Các kết mơ hình thể quy luật diễn biến hệ thống hoàn lưu vùng nước ven bờ Mơ hình xây dựng sử dụng cơng cụ tính tốn cung cấp kết tính trường động lực chi tiết cho mơ hình tính tốn vận chuyển trầm tích, lan truyền chất nhiễm bẩn sinh thái, đồng thời mơ hình sử dụng - 169 - LUAN VAN CHAT LUONG download : add luanvanchat@agmail.com hệ mở cho phép tích hợp mơ hình khác thơng qua mođun bình lưu khuyếch tán với thuật giải có độ xác cao xây dựng Nghiên cứu, áp dụng thành công thuật giải sai phân suy giảm biến động tổng cộng (TVD) phương pháp phân tách đối xứng “bước nhỏ” Yanenko theo thời gian tính tốn q trình bình lưu - đối lưu, nhằm tăng độ ổn định bậc xác thuật giải mơ hình mơ chế độ thủy động lực khu vực với biến có gradien lớn Nghiên cứu áp dụng phương pháp đặc trưng thiết lập phương pháp đặt điều kiện cho biên lỏng biển, biên cửa sông điều kiện xáo trộn đồng phân tầng Trong điều kiện thiếu thốn số liệu đo mực nước biên lỏng việc kết hợp số liệu thực đo mực nước điểm biên cứng điều kiện lý thuyết biến dạng sóng tiến Kenvin áp dụng cho biên lỏng vng góc với bờ góp phần xác định tổ hợp biên phù hợp sử dụng cho mô hình vùng nước nơng ven bờ với tỉ lệ biên lỏng lớn Tập hợp sử dụng thành công sơ đồ khép kín rối từ biểu thức đại số đến phương pháp khép kín bậc cao k-epsilon, qua trình thử nghiệm số lựa chọn sơ đồ khép kín rối phù hợp với trình cần mơ Kết tính tốn ứng dụng mơ hình đưa sơ đồ chi tiết hoàn lưu nước ba chiều vùng biển Hội An - Cù Lao Chàm, Quảng Nam, góp phần cung cấp thơng tin phục vụ mục đích nghiên cứu q trình vận chuyển trầm tích, lan truyền chất ô nhiễm phục vụ mục đích kinh tế công tác quy hoạch vùng nuôi trồng thủy sản, quy hoạch kinh tế du lịch sinh thái Những vấn đề cần tiếp tục nghiên cứu + Tiếp tục phát triển mơ hình theo hướng tích hợp tồn tác động sóng trọng lực bề mặt lên trường dịng chảy ba chiều tiếp cận vào vùng nước sát bờ với độ phân giải cao phục vụ nghiên cứu trình vận chuyển trâm tích - 170 - LUAN VAN CHAT LUONG download : add luanvanchat@agmail.com + Xây dựng hoàn thiện phương án đo đạc thu thập số liệu có định hướng phục vụ hiệu chỉnh nghiệm chứng mơ hình số trị thủy động lực vận chuyển trầm tích áp dụng cho vùng nước ven bờ CÁC CƠNG TRÌNH Đà CƠNG BỐ CỦA TÁC GIẢ Phạm Văn Huấn, Nguyễn Tài Hợi, Nguyễn Minh Huấn Ứng dụng phương pháp bình phương nhỏ vào phân tích thủy triều dịng triều Khí tượng thủy văn - 171 - LUAN VAN CHAT LUONG download : add luanvanchat@agmail.com Biển Đơng Tổng cục Khí tượng Thủy văn, Trung tâm Khí tượng Thủy văn biển, NXB Thống kê, Hà Nội, 2000, tr 196-203 Nguyen Minh Huan A three-dimentional simulation of the tidally modulated plume in the river entrance region Vietnam National University, Hanoi Journal of Science, Natural sciences and Technology, T XIX, No1, 2003, pp 30-42 Phạm Văn Huấn, Đinh Văn Ưu, Nguyễn Minh Huấn, Đoàn Văn Bộ Các đặc trưng rối lớp biên sóng – dịng sát đáy vùng biển ven bờ Đại học Quốc Gia Hà Nội, Tạp chí Khoa học, Khoa học Tự nhiên Cơng nghệ, Phụ trương ngành Khí tượng Thủy văn Hải dương học T XIX , No1 , 2003, tr 39-46 - 172 - LUAN VAN CHAT LUONG download : add luanvanchat@agmail.com TÀI LIỆU THAM KHẢO [1] Công ty tư vấn xây dựng cảng đường thủy (TEDIport), 1999 Báo cáo khảo sát dòng chảy, nhiệt độ, độ mặn gió ven biển miền Trung (1998 – 1999) Tổng cơng ty tư vấn thiết kế G.T.V.T Hà Nội [2] Nguyễn Minh Huấn nnk., 2002 Báo cáo chuyên đề kết khảo sát mực nước, dòng chảy, chất lượng nước khu vực Hội An - Cù Lao Chàm Thuộc đề tài cấp Nhà nước mã số KC 09 12 Luận chứng khoa học mơ hình phát triển kinh tế – sinh thái số đảo, cụm đảo lựa chọn vùng biển ven bờ Việt Nam, chủ nhiệm GS TS Lê Đức Tố Đại học Quốc gia Hà Nội [3] Phạm Văn Huấn, Nguyễn Tài Hợi, Nguyễn Minh Huấn, 2000 Ứng dụng phương pháp bình phương nhỏ vào phân tích thủy triều dịng triều Khí tượng thủy văn Biển Đơng Tổng cục Khí tượng Thủy văn, Trung tâm Khí tượng Thủy văn biển, NXB Thống kê, Hà Nội, tr 196-203 [4] Bùi Hồng Long, 1999 Báo cáo chuyên đề kết nghiên cứu điều kiện động lực, trao đổi nước vật liệu vùng biển cửa Đại – Hội An hai đợt khảo sát 7/1997 – 5/1998 Thuộc đề tài cấp Nhà nước mã số KHCN 06 08 Nghiên cứu quy luật dự đốn xu bồi tụ – xói lở vùng ven biển cửa sông Việt Nam, chủ nhiệm đề tài TS Lê Phước Trình Trung tâm Khoa học Cơng nghệ quốc gia [5] Nguyễn Thị Bình Minh nnk., 2002 Báo cáo chuyên đề chế độ khí hậu, khí tượng khu vực đảo Lý Sơn Đà Nẵng Thuộc đề tài cấp Nhà nước mã số KC 09 12 Luận chứng khoa học mô hình phát triển kinh tế – sinh thái số đảo, cụm đảo lựa chọn vùng biển ven bờ Việt Nam, chủ nhiệm GS TS Lê Đức Tố Đại học Quốc gia Hà Nội [6] Nguyễn Hữu Nhân., 1999 Báo cáo kết nghiên cứu hoàn lưu nước vịnh Thái Lan Thuộc Đề tài cấp nhà nước mã số KHCN 06 03 Điều tra bổ sung vùng biển vịnh Thái Lan, chủ nhiệm đề tài TS Phan Văn Hoặc Tổng cục Khí tượng Thủy văn [7] Số liệu Khí tượng Thủy văn Việt Nam, 1989 Tập Số liệu thủy văn Chương trình tiến khoa học kỹ thuật cấp Nhà nước 42A - 173 - LUAN VAN CHAT LUONG download : add luanvanchat@agmail.com [8] Trung tâm Khí tượng Thủy văn biển Báo cáo tổng kết (1993 – 1994) Đề tài Điều tra nghiên cứu quy hoạch bảo vệ thổ cư, môi trường khai thác tiềm nông ngư nghiệp khu vực Cửa Đại – Hội An Hà Nội - Đà Nẵng [9] Trung tâm Khí tượng Thủy văn biển Báo cáo tổng kết (1994 – 1995) Đề tài Điều tra nghiên cứu quy hoạch bảo vệ thổ cư, môi trường khai thác tiềm nông ngư nghiệp khu vực Cửa Đại – Hội An Hà Nội - Đà Nẵng [10] Đinh Văn Ưu nnk., 2000 Báo cáo chun đề mơ hình chiều cấu trúc hồn lưu nhiệt muối Biển Đơng Thuộc đề tài cấp nhà nước KHCN 06 02 Nghiên cứu cấu trúc ba chiều (3D) thuỷ nhiệt động lực học Biển Đông ứng dụng chúng, chủ nhiệm đề tài PGS TS Đinh Văn Ưu Đại học Quốc gia Hà Nội [11] Baumert H., Burchard H and Kleine E., 1997 On second-moment closures for marine turbulence: A review Presented at the 29th Li`ege Colloquium on Marine Hydrodynamics Turbulence Revisited, Li`ege, May 5–9 1997 [12] Blackadar A.K., 1962 The vertical distribution of wind and turbulent exchange in a neutral atmosphere Journal of Geophysical Research, 67, 3095–3102 [13] Blumberg A.F and Mellor G.L., 1987 A description of a three-dimensional coastal ocean circulation model In : N.S Heaps (Editor), Three-dimensional Coastal Ocean Models Coastal and Estuarine Sciences, Vol 4, American Geophysical Union, Washington D.C., pp 1–16 [14] Bowden K.F., Fairbairn L.A and Hughes P., 1959 The distribution of shearing stresses in a tidal current Geophysical Journal of the Royal Astronomical Society, 2, 288–305 [15] Coastal Engineering Manual 2001 Part II, Ch – Ch U.S Army Corps of Engineers Washington - 174 - LUAN VAN CHAT LUONG download : add luanvanchat@agmail.com [16] Davies A.M., 1993 A bottom boundary layer-resolving three-dimensional tidal model : A sensitivity study of eddy viscosity formulation Journal of Physical Oceanography, 23, 1437–1453 [17] Deleersnijder E., Beckers J.M., Campin J.M., El Mohajir M., Fichefet T and Luyten P., 1997 Some mathematical problems associated with the development and use of marine models In: J.I Diaz (Editor), The mathematics of models for climatology and environment Springer Verlag, Heidelberg, 39–86 [18] Gill A.E., 1982 Atmosphere-Ocean Dynamics International Geophysics Series, Vol.30 Academic Press, Orlando [19] Grant W.D and Madsen O.S., 1979 Combined wave and current interaction with a rough bottom Journal of Geophysical Research, 84, 1797–1808 [20] Hassid S and Galperin B., 1983 A turbulent energy model for geophysical flows Boundary Layer Meteorology, 26, 397–412 [21] Heathershaw A.D., 1981 Comparisons of measured and predicted sediment transport rates in tidal currents Marine Geology, 42, 75–104 [22] Hedstrom G.W., 1979 Nonreflecting boundary conditions for nonlinear hyperbolic systems Journal of Computational Physics, 30, 222–237 [23] Hinze, J., 1975 Turbulence, McGraw Hill Book Co., New York, N.Y [24] Hirsch C., 1988 Numerical computation of internal and external flows Volume 1: Fundamentals of numerical discretisation Wiley, New York, 515 pp [25] Hirsch C., 1990 Numerical computation of internal and external flows Volume 2: Computational methods for inviscid and viscous flows Wiley, New York, 691 pp [26] Huan, N M., 2003 A three-dimentional simulation of the tidally modulated plume in the river entrance region Vietnam National University, Hanoi Journal of Science, Natural sciences and Technology, T XIX, No1, pp 30-42 - 175 - LUAN VAN CHAT LUONG download : add luanvanchat@agmail.com [27] Kowalik, Z., Murty, T S., 1993 Numerical modeling of ocean dynamics World Scientific Publishing, London [28] Kirby, J T., and Dalrymple, R A 1994 Combined Refraction/Diffraction Model REF/DIF 1, Version 2.5 Res Report CACR-94-22, Center for Applied Coastal Research, University of Dalaware [29] Launder B.E and Spalding D.B., 1974 The numerical computation of turbulent flows Computer Methods in Applied Mechanics and Engineering, 3, 269–289 [30] LeVeque R J., Conservation Law Package 4.0 Department of Applied Mathematics University of Washington [31] LeVeque R J., 2002 Finite Volume Methods for Hyperbolic Problems, Cambridge University Press [32] Lions J L., 1968 Quelques méthodes de resolution des problèmes aux limites nolinear Dunod, Paris [33] Luyten P.J., Deleersnijder E., Ozer J and Ruddick K.G., 1996 Presentation of a family of turbulence closure models for stratified shallow water flows and preliminary application to the Rhine out flow region Continental Shelf Research, 16, 101–130 [34] Manh, D V and T Yanagi, 1997 A three-dimentional numerical model of tides anf tidal currents in the Gulf of Tongking La mer, 35 [35] Manh, D V and T Yanagi, 2000 A study on residual flow in the Gulf of Tongking Journal of Oceanography, Vol 56 [36] Mellor G.L and Yamada T., 1974 A hierarchy of turbulence closure models for planetary boundary layers Journal of the Atmospheric Sciences, 31, 1791–1806 [37] Mellor G.L and Yamada T., 1982 Development of a turbulence closure model for geophysical fluid problems Reviews of Geophysics and Space Physics, 20, 851–875 - 176 - LUAN VAN CHAT LUONG download : add luanvanchat@agmail.com [38] Monin, A., and Yaglom, A., 1975 Statistical fluit mechanics, Vols and MIT Press, Cambridge [39] Munk W.H and Anderson E.R., 1948 Notes on a theory of the thermocline Journal of Marine Research, VII(3), 276–295 [40] Oey L.-Y and Chen P., 1992 A model simulation of circulation in the Northeast Atlantic shelves and seas Journal of Geophysical Research, 97, 20087–20115 [41] Pacanowski R.C and Philander S.G.H., 1981 Parameterization of vertical mixing in numerical models of tropical oceans Journal of Physical Oceanography, 11, 1443–1451 [42] Phillips N.A., 1957 A coordinate system having some special advantages for numerical forecasting Journal of Meteorology, 14, 184–185 [43] Polman, T., 1987 A three dimentional circulation model of the South China sea in: Three-dimentional models of marine and estuarine dynamics, J Nihoul and B Jamart, Eds., Elsevier Publ Co., Amsterdam [44] Press W.H., Flannery B.P., Teukolsky S.A and Vetterling W.T., 1989 Numerical Recipes The art of scientific computing Cambridge University Press, Cambridge, 702 pp [45] Rodi W., 1984 Turbulence models and their application in hydraulics International Association for Hydraulic Research, 2nd edition, Delft, Netherlands, 104 pp [46] Roed L.P and Cooper C.K., 1987 A study of various open boundary conditions for wind-forced barotropic numerical ocean models In: J.C.J Nihoul and B.M Jamart (Editors), Three-dimensional models of marine and estuarine dynamics Elsevier, Amsterdam, pp 305–335 [47] Rouse, H., 1976 Advanced mechanics of fluids, Robert E Krieger Publishers, Huntington, N Y - 177 - LUAN VAN CHAT LUONG download : add luanvanchat@agmail.com [48] Semtner A.J and Chervin R.M., 1988 A simulation of the global ocean circulation with residual eddies Journal of Geophysical Research, 93, 15502–15522 [49] Shaw, P T and S Y Chao 1994 Surface circulation in the South China sea Deep-Sea Res., 41 [50] Signell R.P., Beardsley R.C., Graber H.C and Capotondi A., 1990 Effect of wave-current interaction on wind-driven circulation in narrow, shallow embayments Journal of Geophysical Research, 95, 9671–9678 [51] Simpson J.H., Crawford W.R., Rippeth T.P., Campbell A.R and Cheok J.V.S., 1996 The vertical structure of turbulent dissipation in shelf seas Journal of Physical Oceanography, 26, 1579–1590 [52] Smagorinsky J., 1963 General circulation experiments with the primitive quations – I The basic experiment Monthly Weather Review, 91, 99–165 [53] Smith, J M., Sherlock, A R., Resio D T 2001 STWAVE: Steady-state spectral wave model User’s manual for STWAVE, version 3.0 Report ERDC/CHL SR-01-1, Coastal and Hydraulics Laboratory, US Army Corps of Engineers [54] Soulsby R.L., 1983 The bottom boundary layers of shelf seas In: B Johns (Editor), Physical Oceanography of Coastal and Shelf Seas Elsevier Oceanography Series, Vol 35, 189–266 [55] Soulsby R.L., 1997 Dynamics of Marine sands Thomas Telford Publications, London [56] Steinhorn I., 1991 Salt flux and evaporation Journal of Physical Oceanography, 21, 1681–1683 [57] Temman R., 1979 Navier-Stokes equations Theory and numerical analysis North Holland - 178 - LUAN VAN CHAT LUONG download : add luanvanchat@agmail.com [58] UNESCO, 1981 Tenth report of the joint panel on oceanographic tables and standards UNESCO Technical Papers in Marine Science No 36, UNESCO, Paris [59] Van Rjin Leo, 1989 Principles of fluid flow and surface in rivers, estuaries, seas and oceans, Aqua Publications, pp 197-201, 216-219 [60] Xing J and Davies A.M., 1996 Application of turbulence energy models to the computation of tidal currents and mixing intensities in shelf edge regions Journal of Physical Oceanography, 26, 417–447 [61]   , 1982      -            -   [62]   .,    1992         .1992  32  [63]    1991 -            .1991  31  [64]    1992               .1992  32  [65]   ., 1973           ,  - 179 - LUAN VAN CHAT LUONG download : add luanvanchat@agmail.com ... QUAN TÌNH HÌNH NGHIÊN CỨU DỊNG CHẢY BA CHIỀU VÙNG BIỂN VEN BỜ 1.1 Các phƣơng pháp nghiên cứu dịng chảy ba chiều 1.2 Tình hình nghiên cứu giới nƣớc 20 1.3 Các vấn đề nghiên cứu luận án 23 CHƢƠNG... chiều vùng biển ven bờ trình bày cách tổng quát phương pháp, kết nghiên cứu giới nước cấu trúc ba chiều dòng chảy biển vùng nước ven bờ Chương II Cơ sở tốn học mơ hình thuỷ động lực ba chiều. .. tốc theo chiều thẳng đứng lớn 1.1.3.2 Hệ phương trình thủy động lực ba chiều áp dụng cho vùng nước nơng ven bờ Hệ phương trình mơ hình thủy động lực ba chiều áp dụng cho vùng nước nông ven bờ sử