I. Khám phá mô hình hai chiều tính dòng chảy do sóng ven bờ
Việc nghiên cứu và dự báo các quá trình thủy động lực học tại khu vực ven bờ là một nhiệm vụ cốt lõi trong kỹ thuật biển và quản lý tài nguyên ven biển. Khu vực này chịu tác động phức tạp của nhiều yếu tố, trong đó sóng biển đóng vai trò chủ đạo. Dòng chảy do sóng là một hiện tượng quan trọng, ảnh hưởng trực tiếp đến sự vận chuyển bùn cát, gây ra biến đổi đường bờ, và tác động đến sự ổn định của các công trình ven biển. Hiểu rõ và mô phỏng chính xác hiện tượng này là yêu cầu cấp thiết. Luận án tiến sĩ về mô hình hai chiều tính dòng chảy do sóng đã ra đời để giải quyết bài toán này. Mô hình cung cấp một công cụ mạnh mẽ để phân tích và dự báo các đặc trưng dòng chảy trong không gian hai chiều (2D), bao gồm cả dòng chảy dọc bờ và dòng chảy ngang bờ. Mục tiêu chính của nghiên cứu là xây dựng một mô hình số ổn định, chính xác và hiệu quả, có khả năng mô phỏng sự hình thành và phát triển của dòng chảy dưới tác động của trường sóng biến đổi. Khác với các mô hình một chiều chỉ xem xét được dòng chảy dọc bờ, mô hình hai chiều tính dòng chảy do sóng cho phép khảo sát cấu trúc phức tạp của dòng chảy, đặc biệt là sự hình thành các dòng rip (dòng xa bờ) nguy hiểm. Luận án này tập trung vào việc áp dụng các phương trình thủy động lực học tiên tiến và các thuật toán số hiệu quả để đảm bảo kết quả mô phỏng phản ánh đúng thực tế vật lý. Kết quả nghiên cứu không chỉ có ý nghĩa khoa học sâu sắc trong lĩnh vực cơ học chất lỏng mà còn mang lại giá trị ứng dụng to lớn, hỗ trợ các nhà quy hoạch và kỹ sư trong việc thiết kế công trình bảo vệ bờ, dự báo xói lở và quản lý vùng bờ một cách bền vững. Việc phát triển thành công mô hình này là một bước tiến quan trọng trong việc làm chủ công nghệ dự báo các quá trình động lực ven bờ tại Việt Nam.
1.1. Tầm quan trọng của thủy động lực học khu vực gần bờ
Khu vực gần bờ là vùng chuyển tiếp động lực giữa biển và đất liền. Mọi hoạt động kinh tế - xã hội tại đây đều chịu ảnh hưởng trực tiếp từ các quá trình thủy động lực học. Việc hiểu rõ các quy luật này là nền tảng để quy hoạch cảng biển, khu du lịch, khu nuôi trồng thủy sản và xây dựng các công trình phòng chống thiên tai. Dòng chảy do sóng là tác nhân chính gây ra vận chuyển bùn cát, dẫn đến hiện tượng bồi lấp luồng lạch hoặc xói lở bờ biển nghiêm trọng. Một mô hình hai chiều tính dòng chảy do sóng chính xác sẽ giúp dự báo các khu vực có nguy cơ xói lở cao, từ đó đề xuất giải pháp công trình phù hợp. Hơn nữa, dòng chảy còn ảnh hưởng đến việc phát tán chất ô nhiễm từ các cửa sông ra biển, là yếu tố quan trọng trong các bài toán quản lý môi trường biển.
1.2. Tổng quan các phương pháp mô phỏng dòng chảy do sóng
Lịch sử nghiên cứu đã ghi nhận nhiều phương pháp mô phỏng dòng chảy do sóng. Ban đầu là các mô hình vật lý trong phòng thí nghiệm, cho kết quả trực quan nhưng tốn kém và khó mở rộng quy mô. Tiếp theo là các mô hình giải tích, chỉ áp dụng được cho các trường hợp địa hình và điều kiện sóng đơn giản. Sự phát triển của máy tính đã mở ra kỷ nguyên của các mô hình số. Các mô hình một chiều (1D) ra đời sớm, tập trung tính toán dòng chảy ven bờ (longshore current) nhưng bỏ qua các thành phần dòng chảy ngang bờ. Để khắc phục hạn chế này, các nhà khoa học đã phát triển mô hình hai chiều tính dòng chảy do sóng, có khả năng mô phỏng toàn diện hơn trường dòng chảy phức tạp trên mặt bằng, mang lại cái nhìn tổng thể và chính xác hơn về động lực học vùng ven bờ.
II. Thách thức trong việc tính dòng chảy do sóng chính xác
Việc mô phỏng chính xác dòng chảy do sóng trong khu vực gần bờ phải đối mặt với nhiều thách thức lớn về mặt khoa học và kỹ thuật. Thách thức đầu tiên đến từ sự phức tạp của bản chất vật lý. Quá trình tương tác giữa sóng, dòng chảy và địa hình đáy biển là một hệ thống phi tuyến tính cao. Khi sóng truyền vào bờ, nó bị biến dạng do khúc xạ, nhiễu xạ, và tiêu hao năng lượng do ma sát đáy và hiện tượng sóng vỡ. Chính sự thay đổi không đều của năng lượng sóng đã tạo ra một gradient của ứng suất bức xạ, là nguyên nhân gốc rễ hình thành dòng chảy. Việc tính toán chính xác sự biến đổi của ứng suất bức xạ đòi hỏi một mô hình sóng phức tạp và đáng tin cậy. Thách thức thứ hai là việc mô tả đúng hiện tượng sóng vỡ. Sóng vỡ là một quá trình hỗn loạn, giải phóng năng lượng mạnh mẽ và là yếu tố quyết định cường độ của dòng chảy. Các mô hình hiện có thường sử dụng các tham số hóa thực nghiệm để mô tả sóng vỡ, điều này có thể đưa vào các sai số đáng kể. Thách thức thứ ba liên quan đến các yếu tố số học trong mô hình số. Việc lựa chọn sơ đồ số, điều kiện biên, và lưới tính toán đều ảnh hưởng đến độ ổn định và chính xác của mô hình. Một lưới tính quá thô sẽ không thể nắm bắt được các cấu trúc dòng chảy quy mô nhỏ, trong khi lưới quá mịn lại đòi hỏi tài nguyên tính toán khổng lồ. Việc xây dựng một mô hình hai chiều tính dòng chảy do sóng hiệu quả đòi hỏi phải cân bằng giữa độ chính xác vật lý và hiệu suất tính toán. Luận án này đã tập trung giải quyết các thách thức trên bằng cách kết hợp một mô hình sóng tiên tiến với một mô hình dòng chảy 2D, sử dụng các thuật toán số ổn định để đảm bảo kết quả đáng tin cậy.
2.1. Phân tích sự phức tạp của tương tác sóng và địa hình
Tương tác giữa sóng và địa hình là một yếu tố then chốt. Địa hình đáy biển không bằng phẳng, có các dải cát, lạch triều, hay các công trình nhân tạo. Khi sóng truyền qua các địa hình này, chiều cao và hướng sóng thay đổi mạnh mẽ. Sự thay đổi này tạo ra sự phân bố không đồng đều của năng lượng sóng, dẫn đến sự hình thành các hệ thống dòng chảy phức tạp như dòng tuần hoàn (circulation cells) và dòng rip. Việc mô tả chính xác địa hình đáy biển và sự phản ứng của trường sóng với nó là điều kiện tiên quyết cho một mô hình thủy động lực học thành công.
2.2. Hạn chế của các mô hình truyền thống và giả tĩnh
Các mô hình truyền thống thường tách rời việc tính toán biến đổi sóng và tính toán dòng chảy. Chúng thường tính toán trường sóng trước, sau đó sử dụng kết quả (lực gây chảy) để tính dòng chảy. Cách tiếp cận này, gọi là giả tĩnh (quasi-steady), bỏ qua sự tương tác ngược lại của dòng chảy lên sóng. Trong thực tế, dòng chảy có thể làm thay đổi chiều dài và chiều cao sóng, ảnh hưởng trở lại quá trình biến đổi sóng. Các mô hình một chiều thì hoàn toàn bỏ qua các thành phần dòng chảy ngang bờ. Mô hình hai chiều tính dòng chảy do sóng trong luận án này hướng đến việc khắc phục những hạn chế đó bằng cách giải đồng thời hoặc liên kết chặt chẽ hai mô đun sóng và dòng chảy.
III. Hướng dẫn xây dựng nền tảng lý thuyết mô hình hai chiều
Nền tảng của một mô hình hai chiều tính dòng chảy do sóng vững chắc dựa trên các nguyên lý cơ bản của cơ học chất lỏng và thủy động lực học ven bờ. Luận án đã xây dựng mô hình dựa trên hai trụ cột lý thuyết chính: lý thuyết về ứng suất bức xạ và hệ phương trình nước nông. Lý thuyết về ứng suất bức xạ, được tiên phong bởi Longuet-Higgins và Stewart (1964), là khái niệm trung tâm. Nó định nghĩa ứng suất bức xạ là thông lượng động lượng dư thừa do sự tồn tại của sóng. Sự thay đổi trong không gian của các thành phần ứng suất này tạo ra một lực trung bình tác động lên cột nước, gọi là lực gây chảy do sóng. Lực này là động lực chính gây ra dòng chảy ven bờ và các dòng tuần hoàn. Luận án đã trình bày chi tiết cách tính toán các thành phần của ten-xơ ứng suất bức xạ từ các đặc trưng của trường sóng như chiều cao sóng, hướng sóng và độ sâu mực nước. Trụ cột thứ hai là hệ phương trình nước nông (Shallow Water Equations) có xét đến lực gây chảy do sóng. Hệ phương trình này bao gồm phương trình liên tục (bảo toàn khối lượng) và hai phương trình động lượng theo hai phương ngang (bảo toàn động lượng). Trong các phương trình động lượng, các số hạng nguồn quan trọng được bổ sung, bao gồm gradient của ứng suất bức xạ, sức cản đáy do ma sát, và sức cản do rối. Việc lựa chọn công thức tính sức cản đáy và hệ số rối phù hợp là rất quan trọng để mô hình có thể tái hiện chính xác cường độ và sự phân bố của dòng chảy. Luận án đã kế thừa và cải tiến các công thức kinh điển để phù hợp hơn với điều kiện thực tế, tạo nên một cơ sở lý thuyết hoàn chỉnh cho việc phát triển mô hình số.
3.1. Lý thuyết ứng suất bức xạ của Longuet Higgins
Đây là khái niệm đột phá trong việc giải thích cơ chế hình thành dòng chảy do sóng. Trước đây, người ta khó có thể giải thích tại sao sóng (một chuyển động dao động) lại có thể tạo ra một dòng chảy ổn định (chuyển động tịnh tiến). Luận án trích dẫn công trình nền tảng của Longuet-Higgins và Stewart (1964), chỉ ra rằng sóng không chỉ truyền năng lượng mà còn truyền cả động lượng. Khi sóng vỡ, động lượng này được truyền vào cột nước, tạo ra lực đẩy. Ứng suất bức xạ về bản chất là một ten-xơ, có các thành phần khác nhau theo phương truyền sóng và phương vuông góc. Gradient của các thành phần này là nguồn gốc của lực gây chảy. Việc tính toán chính xác chúng là bước đầu tiên và quan trọng nhất trong mô hình hai chiều tính dòng chảy do sóng.
3.2. Hệ phương trình Boussinesq cho biến đổi sóng
Để tính toán được ứng suất bức xạ, trước hết cần có một mô hình sóng chính xác. Các mô hình dựa trên phổ năng lượng thường được sử dụng. Tuy nhiên, đối với các bài toán yêu cầu độ chi tiết cao, luận án có thể đã tham khảo hoặc sử dụng các mô hình dựa trên phương trình Boussinesq. Hệ phương trình này có khả năng mô tả tốt hơn các hiệu ứng phi tuyến và phân tán tần số của sóng khi truyền vào vùng nước nông. Mặc dù đòi hỏi tính toán phức tạp hơn, việc sử dụng một mô hình sóng dựa trên phương trình Boussinesq sẽ cung cấp một trường sóng đầu vào chi tiết, từ đó giúp tính toán lực gây chảy và dòng chảy do sóng một cách chính xác hơn, đặc biệt là trong các khu vực có địa hình phức tạp.
IV. Phương pháp số trong mô hình tính dòng chảy do sóng
Từ nền tảng lý thuyết, việc chuyển hóa các phương trình toán học thành một công cụ tính toán hiệu quả là nhiệm vụ của các phương pháp số. Luận án đã phát triển một mô hình số để giải hệ phương trình nước nông hai chiều. Quá trình này bao gồm ba bước chính: rời rạc hóa không gian và thời gian, xây dựng sơ đồ số để giải các phương trình rời rạc, và xử lý các điều kiện biên. Về rời rạc hóa, miền tính toán được chia thành một lưới các ô chữ nhật hoặc tam giác. Các biến số như độ sâu, vận tốc dòng chảy được xác định tại các điểm nút hoặc tâm ô lưới. Bước thời gian cũng được chia thành các khoảng nhỏ. Sự lựa chọn mật độ lưới và bước thời gian ảnh hưởng trực tiếp đến độ chính xác và tính ổn định của mô hình. Về sơ đồ số, luận án có thể đã sử dụng phương pháp sai phân hữu hạn (Finite Difference Method) hoặc phương pháp phần tử hữu hạn (Finite Element Method). Các phương pháp này chuyển đổi các phương trình vi phân thành một hệ các phương trình đại số tuyến tính có thể giải được bằng máy tính. Sự ổn định của sơ đồ số là yếu tố cực kỳ quan trọng, đảm bảo rằng các sai số nhỏ không bị khuếch đại trong quá trình tính toán. Luận án đã trình bày chi tiết thuật toán giải, bao gồm cả các kỹ thuật xử lý số hạng phi tuyến và đảm bảo tính bảo toàn khối lượng. Cuối cùng, việc xử lý điều kiện biên là rất cần thiết để mô hình hoạt động đúng. Các điều kiện biên xác định cách mô hình tương tác với môi trường bên ngoài, ví dụ như biên ngoài khơi (nơi sóng truyền vào) và biên đường bờ. Việc áp đặt các điều kiện biên hợp lý, như biên bức xạ hay biên tường cứng, đảm bảo rằng sóng và dòng chảy được truyền đi một cách tự nhiên mà không gây ra phản xạ số học giả tạo. Quá trình này hoàn thiện việc xây dựng mô hình hai chiều tính dòng chảy do sóng.
4.1. Rời rạc hóa phương trình bằng phương pháp sai phân
Phương pháp sai phân hữu hạn (FDM) là một lựa chọn phổ biến do tính đơn giản và hiệu quả. Phương pháp này xấp xỉ các đạo hàm trong phương trình vi phân bằng các biểu thức sai phân tại các điểm trên lưới tính. Luận án có thể đã sử dụng một sơ đồ sai phân xen kẽ (staggered grid), nơi các biến vận tốc và mực nước được tính tại các vị trí khác nhau trong một ô lưới. Sơ đồ này giúp tăng cường tính ổn định và tránh được các dao động số học giả. Việc lựa chọn sơ đồ sai phân (tiến, lùi, hay trung tâm) cho các số hạng khác nhau trong phương trình cũng được cân nhắc kỹ lưỡng để tối ưu hóa độ chính xác và hiệu quả.
4.2. Thiết lập điều kiện biên và thuật toán tính toán
Điều kiện biên quyết định sự thành công của một mô phỏng. Tại biên ngoài khơi, các đặc trưng của sóng tới (chiều cao, chu kỳ, hướng) phải được cung cấp. Tại biên đường bờ, điều kiện không thấm (vận tốc pháp tuyến bằng không) thường được áp dụng. Đối với các biên mở, điều kiện biên bức xạ (radiation boundary condition) được sử dụng để cho phép sóng và dòng chảy đi ra khỏi miền tính toán mà không bị phản xạ ngược trở lại. Thuật toán tổng thể thường là một vòng lặp theo thời gian, trong đó ở mỗi bước thời gian, mô hình sẽ cập nhật các giá trị vận tốc và mực nước trên toàn bộ lưới tính dựa trên giá trị của bước thời gian trước đó và các lực tác động.
V. Kết quả ứng dụng mô hình hai chiều tính dòng chảy thực tế
Giá trị của một mô hình số nằm ở khả năng tái hiện các hiện tượng thực tế và đưa ra các dự báo đáng tin cậy. Luận án đã dành một phần quan trọng để kiểm nghiệm và ứng dụng mô hình hai chiều tính dòng chảy do sóng đã xây dựng. Quá trình kiểm nghiệm (validation) được thực hiện bằng cách so sánh kết quả mô phỏng với các bộ dữ liệu đáng tin cậy. Các bộ dữ liệu này có thể đến từ hai nguồn: các thí nghiệm trong phòng thí nghiệm (physical model) và các đợt đo đạc thực địa (field measurement). Việc so sánh với các thí nghiệm trong máng sóng cho phép kiểm tra khả năng của mô hình trong việc mô phỏng các quá trình vật lý cơ bản dưới các điều kiện được kiểm soát chặt chẽ. Việc so sánh với dữ liệu đo đạc tại một khu vực ven biển cụ thể giúp đánh giá hiệu năng của mô hình trong điều kiện thực tế phức tạp. Các chỉ số thống kê như sai số trung bình bình phương (RMSE) hay hệ số tương quan được sử dụng để lượng hóa độ chính xác của mô hình. Sau khi được kiểm nghiệm, mô hình được áp dụng để giải quyết các bài toán thực tiễn. Một trong những ứng dụng quan trọng nhất là dự báo vận chuyển bùn cát và biến đổi đường bờ. Bằng cách kết hợp trường dòng chảy tính toán từ mô hình với một công thức vận chuyển bùn cát, có thể dự báo được các khu vực bồi tụ và xói lở. Kết quả này cung cấp cơ sở khoa học cho việc quy hoạch các giải pháp bảo vệ bờ, như xây dựng đê chắn sóng, mỏ hàn, hay nuôi bãi nhân tạo. Ngoài ra, mô hình còn được sử dụng để nghiên cứu sự hình thành và hoạt động của các dòng rip, một hiện tượng nguy hiểm cho người tắm biển. Việc xác định được các vị trí và điều kiện hình thành dòng rip giúp đưa ra các cảnh báo kịp thời, nâng cao an toàn cho các hoạt động du lịch biển.
5.1. Kiểm nghiệm mô hình qua các bài toán thực nghiệm
Luận án có thể đã sử dụng các bộ dữ liệu thực nghiệm kinh điển, ví dụ như thí nghiệm về dòng chảy ven bờ trên bãi biển phẳng của Kraus và cộng sự, hoặc thí nghiệm về dòng rip trong bể sóng của trường Đại học Delaware. Việc mô phỏng lại thành công các kịch bản này chứng tỏ rằng các quá trình vật lý cốt lõi, như sự biến đổi ứng suất bức xạ và ma sát đáy, đã được mô hình hóa một cách chính xác. Kết quả so sánh thường được trình bày dưới dạng các biểu đồ đối chiếu giữa giá trị tính toán và đo đạc cho vận tốc dòng chảy và mực nước dâng.
5.2. Dự báo vận chuyển bùn cát và biến đổi đường bờ
Vận chuyển bùn cát là hệ quả trực tiếp của dòng chảy do sóng. Mô hình dòng chảy cung cấp trường vận tốc gần đáy, là đầu vào cho các công thức tính toán lượng bùn cát bị dịch chuyển. Luận án có thể đã áp dụng các công thức nổi tiếng như CERC hoặc Van Rijn. Bằng cách tính toán tốc độ vận chuyển bùn cát tại mỗi vị trí, mô hình có thể dự báo sự thay đổi cao độ đáy biển theo thời gian. Sự thay đổi này chính là nguyên nhân gây ra biến đổi đường bờ. Đây là ứng dụng thực tiễn giá trị nhất, giúp đưa ra các kịch bản xói lở/bồi tụ dưới các điều kiện sóng khác nhau.
VI. Đóng góp và hướng phát triển cho mô hình dòng chảy do sóng
Luận án tiến sĩ về mô hình hai chiều tính dòng chảy do sóng đã mang lại những đóng góp khoa học và thực tiễn quan trọng. Về mặt khoa học, luận án đã xây dựng thành công một công cụ tính toán hoàn chỉnh, tích hợp các lý thuyết tiên tiến về thủy động lực học ven bờ và các phương pháp số hiệu quả. Nghiên cứu đã làm rõ hơn về cơ chế hình thành và cấu trúc của dòng chảy do sóng trong các điều kiện địa hình phức tạp, góp phần làm phong phú thêm kho tàng tri thức trong lĩnh vực cơ học chất lỏng ứng dụng. Việc kiểm nghiệm thành công mô hình với các bộ dữ liệu thực nghiệm và thực địa đã khẳng định tính đúng đắn và độ tin cậy của phương pháp luận được áp dụng. Về mặt thực tiễn, mô hình này là một công cụ hữu ích cho các hoạt động quản lý tổng hợp vùng bờ tại Việt Nam. Nó có thể được ứng dụng trong việc đánh giá tác động của biến đổi khí hậu đến xói lở bờ biển, quy hoạch các công trình ven biển, tối ưu hóa thiết kế cảng, và dự báo các rủi ro thiên tai như dòng rip. Sự ra đời của mô hình góp phần nâng cao năng lực tự chủ về công nghệ, giảm sự phụ thuộc vào các phần mềm thương mại đắt đỏ. Tuy nhiên, nghiên cứu khoa học là một quá trình liên tục. Luận án cũng đã chỉ ra những hướng phát triển trong tương lai. Một hướng đi tiềm năng là mở rộng mô hình thành ba chiều (3D) để có thể khảo sát cấu trúc của dòng chảy theo chiều sâu, đặc biệt quan trọng cho việc nghiên cứu dòng chảy dưới đáy (undertow). Một hướng khác là tích hợp mô hình với các yếu tố khác như dòng triều và gió để có một bức tranh động lực tổng thể hơn. Việc áp dụng các kỹ thuật học máy và trí tuệ nhân tạo (AI) để hiệu chỉnh tham số và cải thiện tốc độ tính toán cũng là một triển vọng đầy hứa hẹn. Những kết quả của luận án là nền tảng vững chắc cho các nghiên cứu tiếp theo, hướng tới một hệ thống dự báo động lực ven bờ ngày càng hoàn thiện và chính xác.
6.1. Tổng kết những đóng góp mới của luận án tiến sĩ
Đóng góp chính của luận án là việc xây dựng và kiểm nghiệm thành công một mô hình hai chiều tính dòng chảy do sóng phù hợp với điều kiện Việt Nam. Nghiên cứu đã cải tiến một số tham số hóa trong mô hình, chẳng hạn như hệ số ma sát đáy hoặc mô hình rối, để kết quả mô phỏng gần với thực tế hơn. Luận án cũng đã cung cấp một bộ công cụ có thể được chuyển giao và áp dụng bởi các cơ quan quản lý, các viện nghiên cứu và các công ty tư vấn trong lĩnh vực kỹ thuật biển, tạo ra tác động xã hội và kinh tế trực tiếp.
6.2. Triển vọng tích hợp mô hình 3D và trí tuệ nhân tạo
Tương lai của lĩnh vực mô hình hóa hướng tới sự phức tạp và chính xác cao hơn. Việc phát triển mô hình 3D sẽ cho phép nghiên cứu các hiện tượng không thể mô tả trong không gian 2D, như sự phân bố vận tốc theo phương thẳng đứng. Bên cạnh đó, trí tuệ nhân tạo (AI) và học máy (Machine Learning) có thể được sử dụng để tạo ra các mô hình thay thế (surrogate models) có tốc độ tính toán nhanh hơn hàng nghìn lần, phục vụ cho các hệ thống cảnh báo sớm thời gian thực. Việc kết hợp mô hình số vật lý với AI là hướng đi tất yếu để tạo ra những công cụ dự báo thế hệ mới.
