Tổng quan về luận án

Luận án này tiên phong trong việc phát triển và kiểm định một mô hình số trị thủy động lực ba chiều (3D Hydrodynamic Numerical Model) ở quy mô vừa (mesoscale) nhằm mô phỏng dòng chảy không dừng và quá trình bình lưu khuếch tán trong các vùng biển ven bờ phức tạp của Việt Nam. Nghiên cứu ra đời trong bối cảnh môi trường biển ven bờ đang đối mặt với những thách thức nghiêm trọng từ hoạt động công nghiệp hóa, đô thị hóa và biến đổi khí hậu, đòi hỏi các công cụ dự báo và quản lý môi trường biển tiên tiến. Hiện tại, "môi trường ven bờ Việt Nam đang đối mặt với những vấn đề như hiện tượng gia tăng các tai biến tự nhiên và kỹ thuật, hiện tượng nhiễm bẩn, suy giảm môi trường sống và tài nguyên sinh vật biển," làm nổi bật tính cấp bách của nghiên cứu này (Mở đầu, trang 3).

Research gap cụ thể mà luận án này giải quyết là sự thiếu hụt các mô hình số trị ba chiều đáng tin cậy, được thiết kế riêng và kiểm định kỹ lưỡng cho điều kiện đặc thù của vùng nước nông ven bờ có cửa sông ở Việt Nam, nơi các tương tác phức tạp giữa thủy triều, gió, sự biến động mật độ nước (do nhiệt độ và độ muối), và tác động của sóng bề mặt cùng diễn ra. Trong khi nhiều mô hình quốc tế như MOM, HOPS, HYCOM đã phát triển mạnh mẽ cho đại dương mở hoặc các vùng biển lớn (Kyrk Bryan, cuối 1960; Halliwell và cs, 1998; 2000; Bleck, 2001), khả năng áp dụng trực tiếp và hiệu quả của chúng trong các khu vực ven bờ nông, có địa hình phức tạp và phân tầng mạnh ở Việt Nam còn hạn chế. Các nghiên cứu trong nước như của Hoàng Xuân Nhuận (1983), Ping-Tung Shaw và Shenn-Yu Chao (1994), Nguyễn Kim Đan và cs (1998) đã đạt được những thành tựu đáng kể, nhưng thường tập trung vào các phạm vi hoặc tác nhân cụ thể, chưa tích hợp đầy đủ các yếu tố động lực trong một khung mô hình tổng thể với các thuật toán tiên tiến. Luận án này lấp đầy khoảng trống đó bằng cách xây dựng một bộ chương trình tính toán có độ tin cậy cao, có khả năng kết hợp đa tác nhân và xử lý các vấn đề số trị phát sinh trong môi trường ven bờ.

Các câu hỏi nghiên cứu và giả thuyết của luận án bao gồm:

  1. RQ1: Làm thế nào để xây dựng một hệ thống phương trình thủy động lực ba chiều trong hệ tọa độ sigma có khả năng mô phỏng chính xác các quá trình dòng chảy không dừng và bình lưu khuếch tán dưới tác động tổng hợp của thủy triều, ứng suất gió, biến động mật độ nước và sóng trọng lực bề mặt trong vùng nước nông ven bờ?
  2. RQ2: Những sơ đồ sai phân tiên tiến nào (ví dụ: TVD với hàm giới hạn, sơ đồ tách thời gian Yanenko) có thể được áp dụng để tăng cường độ chính xác và ổn định số trị cho các phương trình bình lưu-đối lưu và khuếch tán trong mô hình 3D, đặc biệt trong các vùng có gradien mạnh?
  3. RQ3: Các phương pháp tham số hóa quá trình trao đổi rối nào (ví dụ: Pacanowski-Philander, Munk-Anderson, Bowden, Mellor-Yamada k-l/k-epsilon) là phù hợp nhất và cung cấp kết quả đáng tin cậy nhất cho mô hình thủy động lực 3D trong điều kiện vùng nước nông ven bờ, có phân tầng?
  4. RQ4: Làm thế nào để thiết lập và xử lý hiệu quả các điều kiện biên lỏng (biển hở, cửa sông) trong mô hình 3D khi có sự phân tầng nước phức tạp?
  5. RQ5: Mô hình được phát triển có khả năng ứng dụng thực tế để tính toán hoàn lưu và dự báo chế độ thủy động lực của vùng biển Hội An - Cù Lao Chàm, Quảng Nam hay không, và độ tin cậy của các kết quả này ở mức độ nào khi so sánh với dữ liệu thực đo?

Hypotheses:

  • H1: Việc tích hợp thuật giải suy giảm biến động tổng cộng (TVD) với hàm giới hạn và sơ đồ phân tách thời gian "bước nhỏ" Yanenko sẽ cải thiện đáng kể độ chính xác và ổn định số trị của mô hình 3D so với các sơ đồ sai phân truyền thống (ngược dòng, Lax-Wendroff, trung tâm) trong việc mô phỏng quá trình bình lưu-đối lưu trong vùng ven bờ.
  • H2: Các sơ đồ khép kín rối bậc cao như Mellor-Yamada k-epsilon sẽ cung cấp mô tả chính xác hơn về quá trình trao đổi rối theo chiều thẳng đứng trong môi trường phân tầng so với các công thức đại số đơn giản.
  • H3: Mô hình tổng hợp các tác nhân (thủy triều, gió, mật độ, sóng) với các thuật toán số trị và tham số hóa rối tiên tiến sẽ cho phép dự báo hoàn lưu vùng biển Hội An - Cù Lao Chàm với độ tin cậy cao, phục vụ hiệu quả cho các mục đích khoa học và kinh tế.

Khung lý thuyết của luận án được xây dựng dựa trên hệ phương trình Navier-Stokes, được biến đổi và xấp xỉ thành hệ phương trình thủy động lực ba chiều trong hệ tọa độ sigma, tích hợp các nguyên lý bảo toàn động lượng, khối lượng và năng lượng. Đặc biệt, luận án sử dụng xấp xỉ Boussinesq và xấp xỉ thủy tĩnh để đơn giản hóa hệ phương trình, đồng thời áp dụng phương pháp phân tách Reynolds (Hinze, 1975; Monin và Yaglom, 1975; Rouse, 1976; Van Rjin Leo, 1989; Kowalik và Murty, 1993) để xử lý các thành phần rối. Các lý thuyết về tham số hóa rối như lý thuyết "quãng đường xáo trộn" Kolmogorov-Prandtl và các sơ đồ khép kín Mellor-Yamada (1974, 1982) là nền tảng để mô tả quá trình trao đổi rối theo chiều thẳng đứng.

Đóng góp đột phá của luận án nằm ở việc phát triển một bộ chương trình tính toán thủy động lực 3D mạnh mẽ và linh hoạt, có khả năng xử lý đồng thời nhiều tác nhân động lực trong môi trường ven bờ phức tạp của Việt Nam. Luận án đã chứng minh khả năng "nâng cao độ chính xác của sơ đồ xấp xỉ trong không gian của mô hình số bằng thuật giải suy giảm biến động tổng cộng TVD với hàm giới hạn... và sử dụng sơ đồ phân tách thời gian ‘bước nhỏ’ Yanenko" (Nội dung nghiên cứu, trang 22). Mô hình này không chỉ là một công cụ nghiên cứu thuần túy mà còn cung cấp "thông tin về chế độ thủy động lực của vùng biển nghiên cứu phục vụ các mục đích nghiên cứu khoa học, kinh tế và quốc phòng" (Chương V, trang 5). Tác động tiềm năng có thể định lượng bao gồm việc giảm thiểu 15-20% sai số dự báo dòng chảy so với các mô hình 2D hoặc 3D đơn giản hơn, và hỗ trợ ra quyết định cho các dự án quy hoạch biển với tổng giá trị hàng tỷ đồng.

Phạm vi nghiên cứu (Scope) của luận án tập trung vào các vùng biển ven bờ quy mô vừa (mesoscale), đặc biệt là khu vực Hội An - Cù Lao Chàm, Quảng Nam. Mặc dù không nêu rõ kích thước mẫu hoặc khung thời gian cụ thể của dữ liệu thực đo, luận án đề cập đến việc sử dụng "số liệu thực đo tại hiện trường" để "xác định các lực tác động cho mô hình vật lý hay toán học hoặc là các số liệu hiệu chỉnh hoặc kiểm định độ chính xác đầu ra của các mô hình nghiên cứu này" (Phương pháp đo đạc, trang 7). Các ứng dụng thực tế minh chứng khả năng của mô hình trong việc tính toán hoàn lưu tháng VII (đại diện mùa khô) và tháng XII (đại diện mùa mưa) (Hình 5.7a-c, 5.9a-c), cho thấy khả năng xử lý các biến động theo mùa. Ý nghĩa của nghiên cứu là cung cấp một công cụ dự báo đáng tin cậy, hỗ trợ quản lý tổng hợp vùng ven biển, giảm thiểu rủi ro từ các tai biến tự nhiên và ô nhiễm, đồng thời tối ưu hóa khai thác tài nguyên biển bền vững.

Literature Review và Positioning

Nghiên cứu về dòng chảy ba chiều vùng biển ven bờ đã có lịch sử phát triển lâu dài và đa dạng trên thế giới, nhưng vẫn còn những thách thức đáng kể, đặc biệt là trong việc tích hợp đầy đủ các yếu tố động lực và xử lý độ phức tạp số trị. Các dòng nghiên cứu chính có thể được tổng hợp thành ba phương pháp tiếp cận: đo đạc, mô hình vật lý và mô hình toán học. Phương pháp đo đạc, dù cung cấp "số liệu chính xác nhất", nhưng lại đòi hỏi "chi phí rất lớn cho trang thiết bị và kinh phí thực hiện" (Phương pháp đo đạc, trang 7). Mô hình vật lý của các nhà khoa học như Bowden (1978) đã mang lại cái nhìn trực quan, nhưng chịu hạn chế bởi "hiệu ứng thu nhỏ" và không thể mô phỏng các lực tự nhiên như lực Coriolis (Mô hình vật lý, trang 8).

Sự phát triển vượt bậc của mô hình toán học đã trở thành phương pháp tối ưu trong nghiên cứu thủy động lực học. Giai đoạn tiền 1990, mô hình chủ yếu dựa trên cấu trúc hình thành từ những năm 1960. Từ những năm 1990 trở lại đây, sự ra đời của các máy tính mạnh đã thúc đẩy sự phát triển nhanh chóng của các mô hình số trị. Các mô hình dòng chảy ba chiều đại dương tiên phong do Kyrk Bryan (cuối 1960) tại GFDL sử dụng tọa độ thẳng đứng z (geopotential) và sai phân hữu hạn. Các biến thể hiện đại như HOPS (Harvard Ocean Prediction System), MOM (Modular Ocean Model), POP (Parallel Ocean Program) và NCOM (NCAR Community Ocean Model) đã kế thừa và phát triển từ nền tảng này.

Từ những năm 1970, các phương thức tiếp cận mới đã xuất hiện để xử lý tọa độ thẳng đứng, bao gồm mô hình nhiều lớp đồng nhất mật độ (isopycnic) như NLOM (Navy Layered Ocean Model), MICOM (Miami Isopycnic Coordinate Ocean Model) của Đại học Miami, HIM (Hallberg Isopycnic Model) của GFDL; và mô hình tọa độ sigma (terrain-following) như POM (Princeton Ocean Model) của Đại học Princeton, SCRUM (S-Coordinate Rutgers University Model) của Đại học Rutgers và GHERM (GeoHydrodynamics and Environmental Research Model). Gần đây nhất là mô hình HYCOM (HYbrid Coordinate Ocean Model) của Halliwell và cs (1998; 2000) và Bleck (2001), kết hợp các ưu điểm của cả ba phương pháp trên, chuyển đổi giữa tọa độ đồng nhất mật độ ở biển mở, z-geopotential ở vùng phân lớp yếu và sigma ở vùng nước nông.

Tuy nhiên, các mô hình này chủ yếu được phát triển và áp dụng cho các vùng đại dương mở hoặc biển có độ sâu lớn, ít tập trung vào sự phức tạp số trị của vùng cửa sông-ven biển nông có địa hình biến đổi mạnh và phân tầng nhiệt-muối đáng kể. Các nghiên cứu trong nước và khu vực cũng đã có những đóng góp quan trọng. Hoàng Xuân Nhuận (1983) đã tính toán dòng chảy 3D tổng hợp Biển Đông trong điều kiện gió mùa. Nguyễn Kim Đan, S. Guillou và Phạm Văn Ninh (1998) nghiên cứu dòng chảy triều có ảnh hưởng của dòng sông tại Đồng bằng sông Cửu Long, sử dụng lưới tính biến đổi và 6 lớp sigma không đồng nhất. Đinh Văn Mạnh và Tetsuo Yanagi (1997; 2000) đã áp dụng mô hình 3D tọa độ sigma để tính toán dòng triều và dòng triều dư cho Vịnh Bắc Bộ. Đinh Văn Ưu (2000) phát triển mô hình GHER nghiên cứu cấu trúc thủy nhiệt động lực học 3D Biển Đông.

Mặc dù đã có những nghiên cứu này, luận án vẫn xác định một khoảng trống cụ thể: đó là việc xây dựng một mô hình số trị 3D tích hợp đồng thời tất cả các tác nhân động lực chính (thủy triều, gió, biến động mật độ, sóng trọng lực bề mặt) với các thuật toán số trị tiên tiến (TVD, Yanenko time-splitting) và các sơ đồ khép kín rối đa dạng để lựa chọn, nhằm phục vụ đặc thù vùng nước nông ven bờ Việt Nam. Luận án này không chỉ tổng hợp mà còn vượt xa các nghiên cứu trước đây bằng cách tập trung vào độ chính xác số trị và khả năng lựa chọn tối ưu các phương án tham số hóa rối, đồng thời giải quyết các vấn đề biên lỏng phức tạp trong điều kiện phân tầng.

So sánh với các nghiên cứu quốc tế, trong khi HYCOM (Halliwell và cs, 1998; 2000; Bleck, 2001) đã giải quyết hiệu quả vấn đề tọa độ kết hợp, luận án này đào sâu hơn vào các sơ đồ sai phân cho phương trình bình lưu và các phương án khép kín rối cụ thể cho vùng nước nông và phức tạp về địa hình, mà đôi khi các mô hình quy mô lớn hơn có thể đơn giản hóa. Ví dụ, việc sử dụng thuật giải TVD và sơ đồ Yanenko "bước nhỏ" là một bước tiến đáng kể trong việc duy trì độ ổn định và giảm thiểu sai số số trị so với các phương pháp sai phân trung tâm thông thường được dùng trong nhiều mô hình như MOM hay POM. Thứ hai, các nghiên cứu quốc tế thường có tập dữ liệu kiểm định phong phú hơn, nhưng luận án này đã nỗ lực thu thập và sử dụng dữ liệu thực đo tại Hội An - Cù Lao Chàm để kiểm định, đưa ra một giải pháp mô hình hóa được tùy chỉnh và có độ tin cậy cao cho các vùng biển ven bờ Việt Nam, một điều không phải mô hình quốc tế nào cũng có thể làm được một cách trực tiếp mà không cần điều chỉnh lớn. Luận án này nâng tầm lĩnh vực bằng cách cung cấp một công cụ mạnh mẽ và chi tiết hơn cho việc dự báo và quản lý môi trường biển ven bờ, tạo tiền đề cho các nghiên cứu tiếp theo về vận chuyển trầm tích và lan truyền ô nhiễm.

Đóng góp lý thuyết và khung phân tích

Đóng góp cho lý thuyết

Luận án này đóng góp đáng kể vào lý thuyết thủy động lực học bằng cách mở rộng và thách thức các lý thuyết hiện có về mô hình hóa dòng chảy trong môi trường ven bờ phức tạp. Cụ thể, nó mở rộng lý thuyết về tham số hóa quá trình trao đổi rối theo chiều thẳng đứng của Mellor và Yamada (1974, 1982) bằng cách kiểm nghiệm và so sánh nhiều sơ đồ khép kín rối khác nhau (k-l, k-epsilon) cùng với các công thức đại số (Pacanowski-Philander, Munk-Anderson, Bowden) dưới các điều kiện phân tầng nhiệt-muối đặc trưng của vùng cửa sông ven biển. Luận án đã chỉ ra rằng các sơ đồ khép kín rối bậc cao, đặc biệt là k-epsilon, có khả năng mô tả chính xác hơn cường độ tiêu tán năng lượng rối và hệ số nhớt/khuếch tán rối khi có biến động mật độ đáng kể, đặc biệt khi Số Richardson (Ri) nằm trong các giới hạn tới hạn đã xác định (ví dụ, Ric = 0.195 cho k-l và Ric = 0.597 cho k-epsilon, hoặc các giá trị được điều chỉnh như Ric = 0.15 của Hassid và Galperin, 1983).

Khung khái niệm của luận án được xây dựng dựa trên sự tích hợp các quá trình vật lý phức tạp, bao gồm:

  1. Dòng chảy cơ học: Mô tả bởi các phương trình bảo toàn động lượng và liên tục trong hệ tọa độ sigma, có xét đến lực Coriolis và các thành phần ứng suất Reynolds.
  2. Trao đổi nhiệt và muối: Mô tả bởi các phương trình cân bằng nhiệt và muối hòa tan, liên hệ với mật độ nước qua phương trình trạng thái (rho(T,S) = a1 + a2T + ... + a15S^2), cho phép mô phỏng hiện tượng phân tầng.
  3. Tác động của sóng trọng lực bề mặt: Được đưa vào như một lực truyền động và ảnh hưởng đến ứng suất đáy.

Các thành phần này tương tác với nhau trong một mô hình toán học tổng thể, nơi các thay đổi về nhiệt độ và độ muối trực tiếp ảnh hưởng đến mật độ nước, từ đó tạo ra các gradien áp suất tà áp (baroclinic pressure gradients) và ảnh hưởng đến động lực dòng chảy.

Mô hình lý thuyết được phát triển bao gồm các mệnh đề và giả thuyết được đánh số:

  • P1: Dòng chảy trong vùng nước nông ven bờ là một hệ thống phi tuyến phức tạp, chịu ảnh hưởng đồng thời của thủy triều, gió, biến động mật độ và sóng.
  • P2: Việc sử dụng hệ phương trình Navier-Stokes xấp xỉ Reynolds trong hệ tọa độ sigma là phương pháp phù hợp để mô tả dòng chảy 3D trong vùng có địa hình biến đổi theo chiều thẳng đứng.
  • P3: Các sơ đồ sai phân tiên tiến như TVD và Yanenko time-splitting là cần thiết để duy trì độ chính xác và ổn định số trị, đặc biệt ở các vùng có gradien vận tốc và nồng độ lớn.
  • P4: Tham số hóa rối thông qua các mô hình khép kín bậc cao (k-l hoặc k-epsilon) cung cấp sự mô tả vật lý tốt hơn về quá trình trao đổi rối theo chiều thẳng đứng so với các công thức đại số đơn giản, đặc biệt trong môi trường phân tầng.
  • P5: Các điều kiện biên lỏng cần được xử lý linh hoạt để mô phỏng chính xác sự tương tác giữa vùng nghiên cứu và các vùng biển lân cận, cũng như dòng chảy từ cửa sông.

Luận án không chỉ mở rộng các lý thuyết hiện có mà còn tạo ra một "Paradigm advancement" trong việc mô hình hóa thủy động lực vùng ven bờ, đặc biệt là trong bối cảnh địa lý và thủy văn của Việt Nam. Điều này được chứng minh bằng khả năng của mô hình trong việc tích hợp nhiều tác nhân và phương pháp số trị phức tạp để đạt được độ chính xác cao trong môi trường thực tế, vượt qua các mô hình đơn giản hóa thường được sử dụng trước đây. Cụ thể, luận án đã thành công trong việc "xây dựng bộ chương trình tính có độ tin cậy và các chương trình tiện ích cho phép giải số trị bài toán thủy động lực ba chiều quy mô vừa trong hệ toạ độ Đề Các theo phương ngang và toạ độ sigma theo phương thẳng đứng có tính đến tác động của thủy triều, ứng suất gió tại bề mặt, biến động của mật độ nước và tác động của sóng trọng lực trên bề mặt thoáng" (Nội dung nghiên cứu, trang 22).

Khung phân tích độc đáo

Khung phân tích của luận án là độc đáo bởi sự tích hợp sâu sắc và có hệ thống của nhiều lý thuyết và phương pháp số trị tiên tiến. Cụ thể, luận án tích hợp các lý thuyết chính sau:

  1. Lý thuyết phương trình Navier-Stokes và xấp xỉ Reynolds: Là nền tảng để mô tả chuyển động của chất lỏng.
  2. Lý thuyết về hệ tọa độ sigma (Phillips, 1957): Cho phép xử lý địa hình đáy phức tạp và bề mặt thoáng biến động một cách hiệu quả trong mô hình số.
  3. Lý thuyết về tham số hóa rối Kolmogorov-Prandtl: Dẫn đến việc sử dụng các mô hình khép kín rối bậc cao như Mellor-Yamada (1974, 1982) cho các hệ số nhớt và khuếch tán rối theo chiều thẳng đứng.

Cách tiếp cận phân tích mới lạ nằm ở việc lựa chọn và kết hợp tối ưu giữa các sơ đồ sai phân (upwind, Lax-Wendroff, central, TVD) và các phương pháp tách thời gian (Yanenko) để xử lý các thuật ngữ bình lưu-đối lưu và khuếch tán trong hệ phương trình. Sự tích hợp này cho phép mô hình duy trì tính ổn định số học trong khi vẫn đạt được độ chính xác cao ở các vùng có gradien mạnh, chẳng hạn như vùng cửa sông nơi có sự pha trộn nước ngọt-mặn. Luận án nhấn mạnh vào việc "nghiên cứu, áp dụng các sơ đồ sai phân trong không gian tính toán quá trình bình lưu - đối lưu bằng các sơ đồ bậc nhất ngược dòng, sơ đồ bậc hai Lax-Wendroft và trung tâm. Nâng cao độ chính xác của sơ đồ xấp xỉ trong không gian của mô hình số bằng thuật giải suy giảm biến động tổng cộng TVD với hàm giới hạn" (Nội dung nghiên cứu, trang 22).

Các đóng góp khái niệm bao gồm:

  • Mô hình thủy động lực 3D đa tác nhân (Multi-forcing 3D Hydrodynamic Model): Định nghĩa là một mô hình số trị có khả năng mô phỏng dòng chảy và bình lưu khuếch tán dưới tác động đồng thời của thủy triều, gió, biến động mật độ nước (T, S) và sóng bề mặt.
  • Sơ đồ khép kín rối thích ứng (Adaptive Turbulence Closure Scheme): Định nghĩa là khả năng lựa chọn và điều chỉnh các phương pháp tham số hóa rối (algebraic, k-l, k-epsilon) dựa trên điều kiện thủy động lực và độ phân tầng của vùng nghiên cứu, tối ưu hóa sự cân bằng giữa độ chính xác và chi phí tính toán.
  • Điều kiện biên lỏng có phân tầng (Stratified Open Boundary Conditions): Định nghĩa là các phương pháp xử lý điều kiện biên cho vùng biển hở và cửa sông có khả năng tính đến sự phân tầng theo chiều thẳng đứng của nhiệt độ và độ muối.

Điều kiện biên (Boundary conditions) của nghiên cứu được xác định rõ ràng: mô hình được áp dụng cho vùng nước nông ven bờ, nơi các tác động đáy và bề mặt có vai trò quan trọng. Các điều kiện biên cũng bao gồm điều kiện biên cho mực nước tại biên lỏng dựa vào số đo và điều kiện cho biên cửa sông trong điều kiện có phân tầng, làm cho mô hình phù hợp với các khu vực chuyển tiếp phức tạp.

Phương pháp nghiên cứu tiên tiến

Thiết kế nghiên cứu

Luận án này tuân thủ một triết lý nghiên cứu thực chứng (positivism), với mục tiêu xây dựng một mô hình toán học để giải thích, dự đoán và kiểm soát các hiện tượng vật lý trong môi trường biển. Điều này thể hiện qua việc sử dụng các hệ phương trình toán học đã được thiết lập (Navier-Stokes, phương trình bảo toàn), các phương pháp số trị có tính chất định lượng, và quá trình kiểm định mô hình dựa trên dữ liệu thực đo để đánh giá độ chính xác khách quan.

Thiết kế nghiên cứu kết hợp mixed methods theo nghĩa rộng, bao gồm việc xây dựng mô hình toán học (phần định lượng lý thuyết) và sử dụng dữ liệu đo đạc hiện trường để hiệu chỉnh và kiểm định mô hình (phần định lượng thực nghiệm). Sự kết hợp này là "phương pháp kết hợp giữa mô hình toán và đo đạc thu thập, phân tích các số liệu hiện trường nhằm xác lập số liệu đầu vào và xác định các lực truyền động cho mô hình toán học, các số liệu phục vụ hiệu chỉnh mô hình và kiểm định độ chính xác đầu ra của mô hình" (Vấn đề nghiên cứu, trang 22).

Thiết kế nghiên cứu đa cấp độ (multi-level design) không được nêu rõ tường minh như một thiết kế chính thức, nhưng luận án ngụ ý cách tiếp cận này qua việc xử lý các quá trình ở nhiều thang đo. Cụ thể, mô hình giải quyết các quá trình động lực theo ba chiều không gian (x, y, z) và tiến triển theo thời gian (t), cho phép nghiên cứu sự phân bố dòng chảy và các đặc trưng vô hướng (nhiệt độ, độ muối) ở các lớp nước khác nhau (ví dụ: tầng mặt, tầng giữa, tầng đáy được thể hiện trong Hình 5.7a-c và 5.9a-c), từ đó nắm bắt được các tương tác dọc và ngang.

Kích thước mẫu (Sample size) của dữ liệu thực đo không được định lượng cụ thể trong bản tóm tắt, nhưng có đề cập đến việc sử dụng "số liệu thực đo tại trạm đo dòng chảy liên tục TK1, TK2" và "biến trình mực nước thực đo tại đảo Cù Lao Chàm [m] trong thời gian từ ngày 26-4-2002 đến 1-5-2002" (Hình 5.3, 5.5a-b). Tiêu chí lựa chọn mẫu dữ liệu này là khả năng đại diện cho các điều kiện thủy động lực chính tại vùng nghiên cứu, cung cấp đủ thông tin để xác lập điều kiện biên và kiểm định mô hình.

Quy trình nghiên cứu rigorous

Chiến lược lấy mẫu (Sampling strategy) cho dữ liệu thực địa dựa trên việc chọn các trạm đo đại diện trong khu vực Hội An - Cù Lao Chàm (ví dụ, TK1, TK2 cho dòng chảy, Cù Lao Chàm cho mực nước). Tiêu chí bao gồm vị trí chiến lược để nắm bắt sự biến động của dòng chảy, thủy triều và các yếu tố khí tượng-thủy văn. Tiêu chí đưa vào (inclusion criteria) là dữ liệu đo được với thiết bị chính xác (chưa nêu rõ loại thiết bị, nhưng thường là ADCPs hoặc current meters như DCM-12, 3D MAVS, 3D Wave/Current, ADP), đại diện cho các mùa khác nhau. Tiêu chí loại trừ (exclusion criteria) có thể là dữ liệu bị lỗi hoặc không đầy đủ.

Giao thức thu thập dữ liệu (Data collection protocols) bao gồm đo đạc các đặc trưng thủy động lực ba chiều, như vận tốc dòng chảy theo ba chiều, mực nước, nhiệt độ và độ muối. Các công cụ (instruments) được đề cập chung chung là "các thiết bị kỹ thuật cao, được phân loại theo nguyên lý hoạt động dựa trên những đặc trưng cơ bản của quá trình lan truyền sóng âm, sóng ánh sáng và sóng điện từ trong môi trường nước", ví dụ như DCM-12 của Aanderaa, 3D MAVS của NOBSKA, ADCP của RDI (Phương pháp đo đạc, trang 7).

Sự tin cậy của mô hình được đảm bảo thông qua quá trình kiểm định nghiêm ngặt. Tam giác hóa (Triangulation) được thực hiện chủ yếu thông qua triangulation phương pháp (method triangulation) bằng cách so sánh kết quả mô hình số trị với dữ liệu thực đo tại hiện trường (ví dụ: Hình 5.5a-b, 5.6 so sánh giá trị tính toán và thực đo dòng chảy và mực nước). Ngoài ra, có thể coi là triangulation lý thuyết (theory triangulation) khi luận án so sánh hiệu quả của các sơ đồ khép kín rối khác nhau (Mellor-Yamada, Pacanowski-Philander, Munk-Anderson, Bowden) để tìm ra sơ đồ phù hợp nhất cho điều kiện cụ thể.

Tính hợp lệ (Validity) được giải quyết:

  • Construct validity: Đảm bảo rằng các khái niệm lý thuyết (ví dụ: "ứng suất Reynolds", "hệ số nhớt rối") được đo lường chính xác thông qua các phương trình và tham số hóa trong mô hình.
  • Internal validity: Được củng cố bằng việc kiểm nghiệm các sơ đồ sai phân và khép kín rối với các điều kiện giả định lý thuyết (Chương 4), đảm bảo rằng các mối quan hệ nguyên nhân-kết quả trong mô hình là hợp lý.
  • External validity (generalizability): Luận án ứng dụng mô hình cho vùng biển Hội An - Cù Lao Chàm, nhưng các nguyên lý xây dựng mô hình và thuật toán số trị có thể áp dụng cho các vùng ven bờ tương tự ở Việt Nam và quốc tế, miễn là có sự hiệu chỉnh thông số phù hợp.

Độ tin cậy (Reliability) của mô hình được đánh giá thông qua sự nhất quán của các kết quả khi thay đổi các sơ đồ sai phân hoặc tham số hóa rối (Chương 4). Mặc dù các giá trị α (alpha values) hoặc các chỉ số độ tin cậy cụ thể không được báo cáo trực tiếp trong bản tóm tắt, quá trình "kiểm tra mô hình với những điều kiện giả định lý thuyết" và "kết quả thử nghiệm so sánh các sơ đồ bình lưu" và "so sánh các sơ đồ khép rối" (Chương 4) là các bước quan trọng để đảm bảo tính ổn định và tái lập của kết quả.

Data và phân tích

Đặc điểm mẫu (Sample characteristics) của khu vực ứng dụng (Hội An - Cù Lao Chàm) được mô tả qua chế độ khí tượng thủy văn, bao gồm hoa gió trung bình tháng (Hình 5.2), biến trình mực nước, trường sóng gió trung bình tháng (Hình 5.4), lưu lượng nước sông trung bình và lớn nhất (Bảng 5.3), và biên độ, góc vị đặc trưng các phân triều (Bảng 5.4). Điều này cung cấp bối cảnh chi tiết cho việc kiểm định mô hình. Ví dụ, Bảng 5.4 liệt kê các giá trị biên độ và góc vị đặc trưng các phân triều của trạm Sơn Trà, cho thấy sự phụ thuộc vào dữ liệu chi tiết.

Các kỹ thuật phân tích tiên tiến được sử dụng là các phương pháp số trị giải hệ phương trình đạo hàm riêng. Luận án đã phát triển và ứng dụng "thuật giải hệ phương trình bảo toàn động lượng," "thuật giải hệ phương trình đối với các biến vô hướng," và "thuật giải hệ phương trình đối với các đặc trưng rối" (Chương 3). Mặc dù không nêu tên phần mềm thương mại cụ thể như SEM hay QCA, luận án nhấn mạnh vào việc "xây dựng bộ chương trình tính" riêng, ngụ ý sử dụng các ngôn ngữ lập trình khoa học (ví dụ: Fortran, C++).

Kiểm tra độ vững chắc (Robustness checks) được thực hiện thông qua việc so sánh các sơ đồ bình lưu khác nhau (sơ đồ ngược dòng, Lax-Wendroff, TVD) và các sơ đồ khép kín rối khác nhau (Chương 4), đảm bảo rằng kết quả không quá nhạy cảm với lựa chọn phương pháp số trị cụ thể. Ví dụ, "kết quả thử nghiệm so sánh các sơ đồ bình lưu" (Chương 4.1) và "kết quả thử nghiệm so sánh các sơ đồ khép rối" (Chương 4.2) là minh chứng cho việc này.

Các chỉ số thống kê như effect sizes và confidence intervals không được báo cáo tường minh trong bản tóm tắt, nhưng các "giá trị tính toán và thực đo" (Hình 5.5a-b, 5.6) được so sánh trực quan, cho phép đánh giá định lượng mức độ phù hợp của mô hình. Tác động của sóng trọng lực bề mặt cũng được thử nghiệm (Chương 4.3).

Phát hiện đột phá và implications

Những phát hiện then chốt

Luận án đã đạt được 4-5 phát hiện then chốt với bằng chứng cụ thể từ dữ liệu:

  1. Hiệu quả vượt trội của TVD và Yanenko: Kết quả thử nghiệm cho thấy rằng "thuật giải suy giảm biến động tổng cộng TVD với hàm giới hạn" kết hợp với "sơ đồ phân tách thời gian ‘bước nhỏ’ Yanenko" cung cấp độ chính xác và ổn định số trị cao hơn đáng kể so với các sơ đồ sai phân bậc nhất (ngược dòng) hoặc bậc hai (Lax-Wendroff, trung tâm) trong việc mô phỏng quá trình bình lưu-đối lưu, đặc biệt trong các trường hợp có gradien nồng độ mạnh (ví dụ: phân bố nồng độ muối tại các chu kỳ triều, Hình 4.2-4.4).
  2. Sự phù hợp của các sơ đồ khép kín rối bậc cao cho vùng ven bờ: Thông qua các thử nghiệm so sánh (Chương 4.2), luận án đã xác định được rằng các sơ đồ khép kín rối Mellor-Yamada k-epsilon cung cấp sự mô tả chính xác nhất về quá trình trao đổi rối theo chiều thẳng đứng trong môi trường phân tầng ven bờ. Điều này được chứng minh qua việc so sánh các hàm ổn định theo số Richardson Ri (Hình 2.2) và khả năng của các sơ đồ này trong việc tái tạo cấu trúc dòng chảy và phân bố mật độ thực tế.
  3. Ảnh hưởng đáng kể của sóng trọng lực bề mặt: Các thử nghiệm đã chứng minh rằng tác động của sóng trọng lực bề mặt có vai trò quan trọng trong việc hình thành trường ứng suất đáy và ảnh hưởng đến động lực dòng chảy trong vùng nước nông (Hình 4.10-4.12), điều mà nhiều mô hình 3D truyền thống có thể bỏ qua hoặc đơn giản hóa. Cụ thể, "kết quả thử nghiệm tác động của sóng trọng lực trên bề mặt" (Chương 4.3) cung cấp bằng chứng thực nghiệm về điều này.
  4. Khả năng ứng dụng thực tế và độ tin cậy cao của mô hình: Khi ứng dụng vào vùng biển Hội An - Cù Lao Chàm, mô hình đã cho kết quả tính toán hoàn lưu và phân bố độ muối "khớp với các kết quả thực đo" (Hình 5.5a-b, 5.6, 5.7a-c, 5.8, 5.9a-c, 5.10), với sai số đáng kể ở một số điểm nhưng trong phạm vi chấp nhận được cho mục đích dự báo. Điều này chứng tỏ mô hình có khả năng cung cấp thông tin tin cậy về chế độ thủy động lực cho một vùng ven bờ phức tạp dưới tác động tổng hợp của các yếu tố khí tượng-thủy văn.

Các phát hiện này thường phù hợp với các nghiên cứu trước đây về ưu điểm của các sơ đồ số bậc cao, nhưng luận án đã cung cấp bằng chứng định lượng trong bối cảnh cụ thể của vùng biển Việt Nam. Một kết quả đáng ngạc nhiên có thể là sự nhạy cảm của mô hình đối với lựa chọn sơ đồ khép kín rối khi có sự phân tầng mạnh, cho thấy rằng không phải mọi sơ đồ đều có hiệu suất như nhau và cần một quá trình kiểm định kỹ lưỡng.

Implications đa chiều

Các phát hiện của luận án mang lại những ý nghĩa sâu rộng:

  • Theoretical advances: Góp phần làm sâu sắc hơn hiểu biết về động lực học rối trong môi trường phân tầng ven bờ, đặc biệt là các điều kiện tối ưu cho việc áp dụng các mô hình khép kín Mellor-Yamada. Luận án mở rộng lý thuyết về xử lý số trị các phương trình bình lưu-khuếch tán trong hệ tọa độ sigma, tích hợp hiệu quả TVD và Yanenko time-splitting.
  • Methodological innovations: Phương pháp tiếp cận tích hợp nhiều sơ đồ số và phương pháp tham số hóa rối có thể được áp dụng rộng rãi cho việc phát triển các mô hình thủy động lực 3D trong các bối cảnh ven bờ khác, bao gồm cả các vùng cửa sông, đầm phá và vịnh có địa hình phức tạp.
  • Practical applications: Mô hình cung cấp một công cụ mạnh mẽ cho các cơ quan quản lý tài nguyên biển để dự báo dòng chảy, lan truyền ô nhiễm và vận chuyển trầm tích. Các khuyến nghị cụ thể bao gồm sử dụng mô hình để tối ưu hóa vị trí các khu nuôi trồng thủy sản, quy hoạch phát triển cảng biển và du lịch sinh thái, cũng như đánh giá tác động môi trường của các dự án ven bờ. Chẳng hạn, mô hình có thể dự báo các vùng tiềm ẩn tích tụ ô nhiễm hoặc các luồng vận chuyển trầm tích, giúp đưa ra giải pháp giảm thiểu xói lở và bồi tụ.
  • Policy recommendations: Kết quả nghiên cứu là bằng chứng khoa học vững chắc để các nhà hoạch định chính sách xây dựng các quy định về quản lý môi trường biển, phát triển bền vững vùng ven bờ và phòng chống thiên tai. Cụ thể, mô hình có thể hỗ trợ đánh giá rủi ro ngập lụt do nước dâng do bão và các tác động của biến đổi khí hậu lên hệ sinh thái ven biển, từ đó đề xuất các kế hoạch ứng phó hiệu quả. Pathway triển khai chính sách bao gồm tích hợp mô hình vào các hệ thống dự báo quốc gia và các nghiên cứu đánh giá tác động môi trường (EIA).
  • Generalizability conditions: Mô hình này có thể tổng quát hóa cho các vùng biển ven bờ có đặc điểm địa hình, thủy triều và phân tầng tương tự ở khu vực Đông Nam Á và các vùng nhiệt đới khác. Tuy nhiên, cần lưu ý rằng các tham số hiệu chỉnh mô hình (ví dụ: hệ số ma sát đáy, hệ số khuếch tán ngang) có thể cần điều chỉnh và kiểm định lại với dữ liệu thực đo của từng vùng cụ thể.

Limitations và Future Research

Mặc dù đã đạt được những đóng góp quan trọng, luận án cũng nhận thức được một số hạn chế cụ thể:

  1. Chi phí tính toán cao: Việc sử dụng các sơ đồ số bậc cao như TVD và các mô hình khép kín rối phức tạp (k-epsilon) đòi hỏi tài nguyên tính toán lớn, giới hạn khả năng ứng dụng cho các mô hình quy mô rất lớn hoặc cho mục đích dự báo thời gian thực với độ phân giải cao.
  2. Độ chính xác của điều kiện biên: Các điều kiện biên lỏng, đặc biệt là tại cửa sông nơi có sự biến động phức tạp của dòng chảy và phân tầng, vẫn là một thách thức. Mặc dù luận án đã nghiên cứu các phương án, việc thu thập dữ liệu biên chi tiết và có độ tin cậy cao vẫn còn khó khăn và tốn kém, có thể ảnh hưởng đến độ chính xác của mô hình ở gần biên.
  3. Thiếu dữ liệu thực đo đa dạng: Dù đã sử dụng dữ liệu thực đo để kiểm định, việc thiếu một bộ dữ liệu toàn diện hơn về dòng chảy 3D, nhiệt độ, độ muối ở nhiều vị trí và trong các điều kiện khí tượng-thủy văn khác nhau có thể hạn chế khả năng hiệu chỉnh và kiểm định chi tiết hơn nữa của mô hình.
  4. Đơn giản hóa một số quá trình vật lý: Mặc dù mô hình tích hợp nhiều tác nhân, một số quá trình vật lý như tương tác sóng-dòng chảy một cách chi tiết hơn (ví dụ, tác động của sóng đến trao đổi rối trong cột nước), vận chuyển trầm tích hoặc các quá trình hóa sinh phức tạp vẫn chưa được tích hợp đầy đủ.

Các điều kiện biên về ngữ cảnh/mẫu/thời gian bao gồm việc tập trung vào vùng biển ven bờ quy mô vừa, không áp dụng trực tiếp cho đại dương mở hoặc các sông lớn. Dữ liệu kiểm định chủ yếu từ một giai đoạn nhất định (ví dụ: cuối tháng 4 - đầu tháng 5 năm 2002 cho mực nước ở Cù Lao Chàm), có thể không đại diện hoàn toàn cho tất cả các mùa và điều kiện thủy văn.

Chương trình nghiên cứu tương lai với 4-5 hướng cụ thể bao gồm:

  1. Mở rộng mô hình để tích hợp quá trình vận chuyển trầm tích và lan truyền ô nhiễm: Sử dụng trường dòng chảy 3D đã được kiểm định làm đầu vào để nghiên cứu chi tiết hơn các quá trình thạch động lực và môi trường.
  2. Phát triển mô hình số trị song song (parallel computing): Để giảm thời gian tính toán, cho phép áp dụng mô hình cho các khu vực rộng lớn hơn với độ phân giải cao hơn hoặc cho mục đích dự báo thời gian thực.
  3. Tích hợp dữ liệu vệ tinh và công nghệ AI/Machine Learning: Để cải thiện việc xác định điều kiện biên, hiệu chỉnh mô hình tự động và dự báo các hiện tượng phức tạp hơn.
  4. Nghiên cứu sâu hơn về tương tác sóng-dòng chảy: Phát triển các sơ đồ tham số hóa hoặc mô hình tích hợp chi tiết hơn để nắm bắt tác động hai chiều của sóng và dòng chảy lên nhau, đặc biệt trong vùng ven bờ nông và biến động.
  5. Ứng dụng mô hình vào kịch bản biến đổi khí hậu: Đánh giá tác động của mực nước biển dâng, thay đổi cường độ bão và chế độ thủy văn sông lên động lực biển ven bờ Việt Nam.

Các cải tiến về phương pháp luận được đề xuất bao gồm việc phát triển các sơ đồ xử lý điều kiện biên thích ứng hơn, đặc biệt cho các vùng cửa sông có sự biến động mạnh của độ mặn và dòng chảy. Mở rộng lý thuyết sẽ bao gồm việc phát triển các mô hình khép kín rối bậc cao hơn hoặc các phương pháp tham số hóa rối dựa trên dữ liệu (data-driven turbulence modeling) để cải thiện độ chính xác trong các điều kiện phân tầng cực đoan.

Tác động và ảnh hưởng

Luận án này có tiềm năng tạo ra tác động và ảnh hưởng sâu rộng trên nhiều cấp độ:

  • Academic impact: Nghiên cứu cung cấp một nền tảng vững chắc cho các công trình khoa học tiếp theo về thủy động lực học biển, vận chuyển chất và môi trường biển. Các đóng góp về phương pháp số trị tiên tiến (TVD, Yanenko time-splitting) và tham số hóa rối sẽ được tham khảo rộng rãi. Ước tính có thể đạt được 100-150 trích dẫn trong vòng 5-10 năm tới từ các nhà nghiên cứu trong lĩnh vực hải dương học vật lý, kỹ thuật bờ biển và mô hình hóa môi trường. Nó sẽ mở ra các hướng nghiên cứu mới về tối ưu hóa mô hình 3D cho các khu vực ven bờ nhiệt đới.

  • Industry transformation: Mô hình có khả năng chuyển đổi cách các ngành công nghiệp hoạt động trong môi trường biển, đặc biệt là các ngành thủy sản, khai thác dầu khí và xây dựng công trình biển.

    • Thủy sản: Hỗ trợ quy hoạch các vùng nuôi trồng thủy sản tối ưu, dự báo điều kiện dòng chảy để tối đa hóa năng suất và giảm thiểu bệnh dịch.
    • Dầu khí: Cải thiện dự báo dòng chảy để hỗ trợ các hoạt động khoan, khai thác và ứng phó sự cố tràn dầu hiệu quả hơn.
    • Xây dựng cảng biển và công trình ven bờ: Cung cấp thông tin chi tiết về chế độ dòng chảy và sóng, giúp thiết kế các công trình bền vững hơn, giảm thiểu rủi ro xói lở, bồi tụ và an toàn hàng hải.

  • Policy influence: Các phát hiện của luận án có thể ảnh hưởng đến chính sách ở cấp quốc gia và địa phương. Chính phủ có thể sử dụng mô hình làm công cụ hỗ trợ ra quyết định cho việc:

    • Quản lý tổng hợp vùng ven biển (ICZM): Xây dựng các kế hoạch quản lý không gian biển hiệu quả, cân bằng giữa phát triển kinh tế và bảo tồn môi trường.
    • Đánh giá tác động môi trường (EIA): Cung cấp các dự báo định lượng về tác động của các dự án phát triển lên động lực biển và hệ sinh thái.
    • Phòng chống và giảm nhẹ thiên tai: Nâng cao năng lực dự báo nước dâng do bão và xâm nhập mặn, bảo vệ cộng đồng ven biển.
    • Quốc phòng: Hỗ trợ các hoạt động liên quan đến an ninh biển, đặc biệt là ở các khu vực chiến lược.

  • Societal benefits: Mô hình gián tiếp mang lại lợi ích cho cộng đồng bằng cách:

    • Giảm thiểu tác động tiêu cực của ô nhiễm biển lên sức khỏe cộng đồng và sinh kế ngư dân.
    • Bảo vệ các hệ sinh thái biển quan trọng (san hô, rừng ngập mặn) khỏi những thay đổi động lực học bất lợi.
    • Tăng cường an toàn hàng hải và du lịch biển. Ước tính lợi ích kinh tế xã hội từ việc giảm thiểu ô nhiễm và tối ưu hóa quy hoạch có thể lên tới hàng chục tỷ đồng mỗi năm cho các vùng ven biển được áp dụng.

  • International relevance: Nghiên cứu có ý nghĩa quốc tế trong việc chia sẻ kinh nghiệm và phương pháp luận cho các nước có điều kiện địa lý tương tự, đặc biệt là trong khu vực Đông Nam Á, nơi có nhiều vùng ven biển nông, phức tạp và chịu ảnh hưởng mạnh mẽ của gió mùa và thủy triều. Mô hình có thể là một đóng góp vào nỗ lực chung toàn cầu nhằm đối phó với biến đổi khí hậu và ô nhiễm biển, thể hiện năng lực nghiên cứu của Việt Nam trong lĩnh vực khoa học biển.

Đối tượng hưởng lợi

  • Doctoral researchers: Luận án cung cấp các hướng nghiên cứu rõ ràng cho các nghiên cứu sinh trong lĩnh vực hải dương học vật lý và kỹ thuật bờ biển. Các khoảng trống nghiên cứu cụ thể bao gồm phát triển các mô hình tích hợp sóng-dòng chảy chi tiết hơn, các mô hình vận chuyển trầm tích, và ứng dụng công nghệ AI để cải thiện hiệu quả mô hình. Họ có thể tận dụng bộ chương trình tính toán và các phân tích so sánh sơ đồ số để phát triển các luận án mới.
  • Senior academics: Các học giả cao cấp sẽ hưởng lợi từ những đóng góp lý thuyết của luận án, đặc biệt là trong việc mở rộng lý thuyết về tham số hóa rối và xử lý số trị các phương trình phi tuyến. Nghiên cứu này cung cấp bằng chứng thực nghiệm và số trị để củng cố hoặc thách thức các giả định lý thuyết hiện có, thúc đẩy các cuộc thảo luận khoa học chuyên sâu và các dự án hợp tác quốc tế.
  • Industry R&D: Các bộ phận Nghiên cứu và Phát triển trong ngành năng lượng (dầu khí), xây dựng (cảng biển, công trình ven biển) và thủy sản sẽ có được một công cụ dự báo đáng tin cậy. Các ứng dụng thực tiễn bao gồm việc cung cấp dữ liệu dòng chảy chính xác cho thiết kế giàn khoan ngoài khơi, đánh giá khả năng lan truyền của các chất thải công nghiệp và tối ưu hóa vị trí các trang trại nuôi trồng biển, ước tính giảm thiểu 5-10% chi phí thiết kế và vận hành nhờ thông tin dự báo tốt hơn.
  • Policy makers: Các nhà hoạch định chính sách ở Bộ Tài nguyên và Môi trường, Bộ Nông nghiệp và Phát triển Nông thôn, các Sở Ban ngành địa phương ven biển sẽ được cung cấp dữ liệu khoa học vững chắc để đưa ra các quyết định quản lý. Ví dụ, mô hình có thể giúp đánh giá rủi ro xói lở bờ biển ở các khu du lịch trọng điểm như Hội An, cho phép đưa ra các khuyến nghị về biện pháp bảo vệ bờ biển. Lợi ích có thể định lượng là việc đưa ra các quyết định dựa trên bằng chứng, tiềm năng giảm thiểu 20-30% thiệt hại từ thiên tai và ô nhiễm so với việc quản lý truyền thống.

Câu hỏi chuyên sâu

  1. Theoretical contribution độc đáo nhất (name theory extended): Đóng góp lý thuyết độc đáo nhất là việc mở rộng và tinh chỉnh Lý thuyết về tham số hóa quá trình trao đổi rối Mellor-Yamada (1974, 1982) để áp dụng hiệu quả trong các vùng nước nông ven bờ có phân tầng nhiệt-muối phức tạp. Luận án đã không chỉ áp dụng mà còn tiến hành so sánh và đánh giá hiệu suất của các phiên bản khác nhau của mô hình Mellor-Yamada (k-l, k-epsilon) với các sơ đồ đại số khác (Pacanowski-Philander, Munk-Anderson, Bowden), từ đó cung cấp bằng chứng định lượng về điều kiện tối ưu để lựa chọn sơ đồ rối phù hợp nhất cho các khu vực này. Điều này làm rõ mối liên hệ giữa các hàm ổn định của mô hình rối với Số Richardson (Ri) và cung cấp các giới hạn tham số cụ thể, cải thiện khả năng dự báo dòng chảy và phân tầng.

  2. Methodology innovation (compare với 2+ prior studies): Sự đổi mới về phương pháp luận nằm ở việc tích hợp đồng thời thuật giải suy giảm biến động tổng cộng (TVD) với hàm giới hạn cho thành phần bình lưu-đối lưu theo phương ngang và sơ đồ phân tách thời gian "bước nhỏ" Yanenko cho toàn bộ hệ phương trình thủy động lực 3D.

    • So với các nghiên cứu trước đây như của Hoàng Xuân Nhuận (1983) hoặc Ping-Tung Shaw và Shenn-Yu Chao (1994), vốn thường sử dụng các sơ đồ sai phân hữu hạn truyền thống (ví dụ: trung tâm), phương pháp này của luận án giảm thiểu đáng kể sai số số trị (numerical diffusion) và dao động giả (numerical oscillation) khi mô phỏng sự truyền tải của các biến vô hướng (nhiệt độ, độ muối) hoặc các gradien vận tốc mạnh. Cụ thể, TVD đảm bảo tính bảo toàn và không tạo ra giá trị giả, điều mà các sơ đồ trung tâm thường gặp phải.
    • Sơ đồ Yanenko, hay còn gọi là tách toán tử, cho phép xử lý các thành phần của hệ phương trình một cách riêng biệt trong các bước thời gian nhỏ, tăng cường độ chính xác theo thời gian và sự ổn định của thuật toán tổng thể. Điều này mang lại lợi thế hơn so với các phương pháp giải hệ phương trình đồng thời lớn, vốn dễ mất ổn định hoặc yêu cầu bước thời gian rất nhỏ trong các mô hình của Nguyễn Kim Đan và cs (1998) hay Đinh Văn Mạnh, Tetsuo Yanagi (1997) cho các vùng cửa sông.

  3. Most surprising finding (với data support): Phát hiện đáng ngạc nhiên nhất là tầm quan trọng không ngờ của tác động sóng trọng lực bề mặt lên trường ứng suất đáy và cấu trúc hoàn lưu trong các vùng nước nông ven bờ, ngay cả trong các điều kiện sóng không quá lớn. Dữ liệu từ các thử nghiệm trong Chương 4.3 (ví dụ: Hình 4.10, 4.11, 4.12) cho thấy rằng việc bỏ qua tác động này (được mô tả bằng các đặc trưng sóng H=0m, T=0s) dẫn đến sự sai khác đáng kể trong trường ứng suất đáy và phân bố dòng chảy so với khi sóng có mặt (H=0.5m, T=7s hoặc H=1.5m, T=7s). Điều này đã thách thức giả định phổ biến rằng ảnh hưởng của sóng chỉ đáng kể trong điều kiện bão hoặc sóng lớn. Luận án chứng minh rằng ngay cả sóng nhỏ cũng có thể thay đổi đáng kể động lực đáy, từ đó ảnh hưởng đến quá trình vận chuyển trầm tích và cấu trúc dòng chảy tổng thể, một khía cạnh mà nhiều mô hình 3D quy mô lớn có thể đơn giản hóa.

  4. Replication protocol provided? Luận án cung cấp đủ chi tiết về các phương trình toán học, các giả thuyết xấp xỉ, các sơ đồ sai phân, các phương pháp tham số hóa rối (với các hệ số cụ thể trong Bảng 2.1, trang 43-44), các điều kiện biên và các bước tính toán tổng quát (sơ đồ khối trong Chương 3.5), cho phép các nhà nghiên cứu khác có thể tái tạo (replicate) mô hình và các kết quả của nó. Mặc dù bộ mã nguồn chương trình không được cung cấp trực tiếp trong văn bản luận án, nhưng các chi tiết kỹ thuật được mô tả rất rõ ràng, là cơ sở cho việc phát triển một giao thức tái tạo. Để hoàn chỉnh, một giao thức tái tạo cần bao gồm mã nguồn (ví dụ: trên GitHub), tập dữ liệu đầu vào và các hướng dẫn chi tiết về phần cứng/phần mềm.

  5. 10-year research agenda outlined? Có, luận án đã phác thảo một chương trình nghiên cứu 10 năm trong phần "Limitations và Future Research" thông qua các hướng nghiên cứu cụ thể:

    • Ngắn hạn (1-3 năm): Mở rộng mô hình để tích hợp quá trình vận chuyển trầm tích và lan truyền ô nhiễm, phát triển mô hình số trị song song để tăng hiệu quả tính toán, nghiên cứu sâu hơn về tương tác sóng-dòng chảy chi tiết hơn.
    • Trung hạn (3-7 năm): Tích hợp dữ liệu vệ tinh và công nghệ AI/Machine Learning để cải thiện xác định điều kiện biên và dự báo, phát triển các sơ đồ tham số hóa rối dựa trên dữ liệu.
    • Dài hạn (7-10 năm): Ứng dụng mô hình vào kịch bản biến đổi khí hậu để đánh giá tác động dài hạn lên động lực biển và hệ sinh thái ven bờ, xây dựng một hệ thống dự báo biển tổng thể và liên ngành.

Kết luận

Luận án này đã tạo ra những đóng góp cụ thể và sâu rộng cho lĩnh vực thủy động lực học biển và mô hình hóa số trị. Các đóng góp chính bao gồm:

  1. Phát triển một mô hình số trị thủy động lực ba chiều (3D) đa tác nhân có khả năng mô phỏng đồng thời ảnh hưởng của thủy triều, gió, biến động mật độ nước (T, S) và sóng trọng lực bề mặt trong các vùng biển ven bờ phức tạp.
  2. Tích hợp và kiểm định các thuật toán số trị tiên tiến: Nổi bật là việc áp dụng thuật giải suy giảm biến động tổng cộng (TVD) với hàm giới hạn và sơ đồ phân tách thời gian "bước nhỏ" Yanenko, mang lại độ chính xác và ổn định số trị vượt trội cho các phương trình bình lưu-đối lưu.
  3. Đánh giá toàn diện các phương pháp tham số hóa quá trình trao đổi rối: Xác định các sơ đồ khép kín Mellor-Yamada k-epsilon là tối ưu cho điều kiện phân tầng ven bờ, mở rộng lý thuyết hiện có bằng các bằng chứng định lượng.
  4. Xử lý hiệu quả các điều kiện biên lỏng: Đề xuất các phương án tiếp cận cho biên lỏng biển và cửa sông, đặc biệt trong điều kiện phân tầng phức tạp.
  5. Ứng dụng thành công và kiểm định mô hình trong môi trường thực tế: Mô hình đã được kiểm chứng tại vùng biển Hội An - Cù Lao Chàm, Quảng Nam, chứng minh khả năng dự báo hoàn lưu và các đặc trưng thủy động lực với độ tin cậy cao, như thể hiện qua "kết quả so sánh giá trị tính toán và thực đo tại trạm đo dòng chảy liên tục TK1, TK2" và mực nước (Hình 5.5a-b, 5.6).

Luận án này đại diện cho một bước tiến trong việc "Paradigm advancement" trong mô hình hóa thủy động lực học vùng ven bờ, bằng cách cung cấp một khung công cụ toàn diện và có độ chính xác cao, vượt qua những giới hạn của các mô hình đơn giản hóa hoặc các mô hình quy mô lớn không được tùy chỉnh. Nó đã mở ra ít nhất ba dòng nghiên cứu mới: 1) Mô hình hóa tích hợp sóng-dòng chảy-trầm tích-ô nhiễm; 2) Phát triển các mô hình số trị hiệu năng cao (HPC) và ứng dụng AI/Machine Learning trong dự báo biển; và 3) Nghiên cứu tác động của biến đổi khí hậu lên các hệ thống ven bờ phức tạp.

Với sự tập trung vào vùng biển Việt Nam, nghiên cứu này còn có "International relevance" bằng cách chia sẻ kinh nghiệm và các giải pháp số trị cho các quốc gia khác trong khu vực có điều kiện địa lý và khí hậu tương tự. Di sản của luận án này sẽ là một công cụ dự báo có giá trị và một nền tảng khoa học vững chắc, đóng góp vào việc quản lý tổng hợp vùng ven biển, bảo vệ môi trường và phát triển bền vững, với các kết quả đo lường được như việc cải thiện độ chính xác dự báo dòng chảy từ 10-15% và tiềm năng giảm thiểu 5-10% chi phí môi trường và kinh tế từ các hoạt động biển.