Tổng quan về luận án

Bờ biển Việt Nam trải dài trên 3.260 km, nằm trong khu vực nhiệt đới gió mùa chịu áp lực khốc liệt từ 6 đến 10 cơn bão nhiệt đới mỗi năm với cường độ ngày càng dị thường và quỹ đạo phức tạp. Dưới tác động kết hợp của sóng bão, triều dâng và dòng chảy ven bờ, hệ thống đê kè biển tại các trọng điểm xung yếu như Cát Hải (Hải Phòng), Thịnh Long (Nam Định) liên tục chịu tổn thất nghiêm trọng. Luận án tiến sĩ kỹ thuật "Nghiên cứu vận tốc dòng chảy do sóng tại chân kè nông trong thiết kế chân kè đá đổ, áp dụng cho đê biển Cát Hải, Hải Phòng" (Chuyên ngành: Kỹ thuật xây dựng công trình thủy, Mã số: 9580202) do NCS. Phạm Văn Lập thực hiện dưới sự hướng dẫn của PGS. Lê Xuân Roanh và GS. Ngô Trí Viềng tại Trường Đại học Thủy lợi, đại diện cho công trình tiên phong giải quyết mắt xích xung yếu nhất trong chuỗi ổn định công trình bảo vệ bờ biển: động lực học dòng chảy đáy tại chân kè nông.

Nghiên cứu nhận diện một khoảng trống học thuật và thực tiễn (research gap) cốt lõi: Các tiêu chuẩn thiết kế đê biển hiện hành tại Việt Nam (như TCVN 9901:2014 và tiêu chuẩn ngành 14TCN130-2002) đang áp dụng công thức xác định vận tốc lớn nhất tại chân công trình dựa trên lý thuyết sóng tuyến tính giản lược của Izbash (1930) hoặc Van Rijn (2018):

$$U_{max} = \frac{\pi H_s}{T_p \sinh(2\pi d_s / L_s)}$$

hoặc công thức kinh nghiệm không thứ nguyên:

$$u_c = 6{,}156\sqrt{d}$$

Các công thức này hoàn toàn bỏ qua tính chất phi đối xứng của quỹ đạo hạt nước trong vùng sóng vỡ, bỏ qua độ dốc lưu vực bãi ($i$), độ dốc mái kè ($m$), chu kỳ đỉnh phổ ($T_p$) và cơ chế tiêu tán năng lượng của mố nhám ($a$). Hệ quả là các tính toán theo quy phạm đánh giá thấp ứng suất cắt đáy thực tế, dẫn đến sự mất ổn định của lớp đá đổ gia cố chân kè. Tại thực địa đê biển Cát Hải, tác giả ghi nhận hiện tượng phá hủy đặc thù:

"Thực tế cho thấy: đá hộc rải trên nền giáp chân kè bị di chuyển, hoặc bị đẩy lên mái, hoặc di động dọc bờ, hoặc rời ra phía biển. Nổi trội hiện nay là hiện tượng đá trượt lên mái, lúc lên, lúc xuống, ma sát giữa đá và mái kè đã gây mài mòn và làm hư hại kết cấu bảo vệ mái."

Để giải quyết triệt để vấn đề, nghiên cứu thiết lập hệ thống 2 câu hỏi nghiên cứu (RQ) và 2 giả thuyết khoa học (H):

  • RQ1: Các thông số thủy động lực học phi tuyến và biên hình học công trình chi phối phân bố vec-tơ vận tốc dòng chảy sát đáy ($U_{max}$) tại chân kè nông như thế nào?
  • RQ2: Quy luật tương tác giữa sóng biển, mố nhám giảm sóng và khối lượng viên đá bảo vệ chân kè được định lượng ra sao nhằm ngăn chặn hiện tượng mài mòn mái đê?
  • H1: Vận tốc dòng chảy đáy cực đại tại chân kè nông không chỉ phụ thuộc vào tỷ số $H_s/d_s$ mà tỷ lệ thuận phi tuyến với chu kỳ sóng $T_p$, độ dốc bãi $i$ và hệ số dốc mái $m$, đòi hỏi phải mô hình hóa bằng lý thuyết sóng phi tuyến bậc cao kết hợp phương trình Navier-Stokes trong môi trường rỗng.
  • H2: Việc tích hợp hệ mố nhám cản đá chuyên dụng có kích thước và khoảng cách tối ưu sẽ tiêu giảm từ 18% đến 35% vận tốc dòng chảy đáy cực đại, đồng thời tạo bẫy cơ học cô lập hoàn toàn sự di chuyển của các viên đá hộc rời rạc lên bề mặt mái kè.

Khung lý thuyết của luận án được xây dựng dựa trên sự giao thoa giữa Thủy động lực học sóng phi tuyến vùng ven bờ (Nonlinear Coastal Hydrodynamics), Cơ học vận chuyển bùn cát và chuyển vị hạt rời (Incipient Motion Mechanics) và Mô hình dòng chảy 2 pha qua môi trường rỗng (Two-Phase Porous Media Flow). Phạm vi nghiên cứu bao quát 320 kịch bản mô phỏng số trị Navier-Stokes 2 pha (MSS-2D) kết hợp hệ thống thí nghiệm mô hình vật lý tỷ lệ 1:10 trên máng sóng hiện đại tại Phòng Thí nghiệm Thủy lực Tổng hợp (Đại học Thủy lợi), áp dụng tính toán kiểm chứng cho 10 mặt cắt đặc trưng tuyến đê biển Cát Hải.

Literature Review và Positioning

Nghiên cứu ổn định công trình bảo vệ bờ biển trước tác động của sóng và dòng chảy ven bờ phát triển qua ba giai đoạn lý thuyết chính:

1. Dòng lý thuyết thủy động lực học sóng sát bờ và vùng sóng vỗ (Swash Zone Dynamics)

Khởi đầu từ lý thuyết sóng tuyến tính Airy (1845) với giả định hạt nước chuyển động theo quỹ đạo elip đối xứng khép kín, các nghiên cứu kinh điển đã nhanh chóng bộc lộ hạn chế khi tiến vào vùng nước nông ($d_s / L_s < 0{,}05$). Stokes (1847) phát triển lý thuyết sóng phi tuyến bậc cao ($n = 1 \div 5$), trong khi Korteweg và de Vries (1895) giới thiệu lý thuyết sóng Cnoidal áp dụng cho số Ursell cao ($U_r = H L^2 / d_s^3 > 48{,}25$). Longuet-Higgins (1953, 1970, 1975) đã đặt nền móng giải tích cho các cấu trúc dòng chảy phức tạp ven bờ, chứng minh sự hình thành của dòng phản hồi sát đáy (undertow) với vận tốc đạt tới $8% \div 10%$ của $\sqrt{g d_s}$ và dòng dọc bờ do gradient ứng suất bức xạ (radiation stress). Ruessink và Van Rijn (2018) bổ sung mô hình dòng leo trong vùng sóng tràn với vận tốc đáy $u_b = \sqrt{g d_s}$ dao động từ $1{,}0\text{ m/s}$ đến $3{,}0\text{ m/s}$.

2. Dòng lý thuyết ngưỡng vận động của hạt vật liệu rời (Incipient Motion of Armor Units)

Izbash (1930) đề xuất phương trình xác định vận tốc tới hạn dựa trên sự cân bằng tĩnh giữa lực cản thủy động ($F_D$), lực nâng ($F_L$) và trọng lượng chìm ($W$). Shields (1936) nâng tầm lý thuyết lớp biên bằng đồ thị thông số không thứ nguyên $\Psi_c$ liên hệ với số Reynolds hạt ($Re_*$). Tuy nhiên, Paintal (1971) chứng minh rằng đối với các hạt thô (đá hộc), ngưỡng vận động không tiệm cận về 0 như giả định của Shields mà tốc độ vận chuyển tỷ lệ thuận với lũy thừa 16 của ứng suất tiếp đáy (tương đương lũy thừa 13 của vận tốc dòng chảy). Pilarczyk (1990) mở rộng bài toán cho dòng chảy xiên trên mái dốc, thiết lập quan hệ giữa chiều dày lớp bảo vệ $D$ và bình phương vận tốc tới hạn $U_{cr}^2$:

$$D = \frac{0{,}035}{\Phi \cdot \Delta \cdot K_s \cdot 2g} U_{cr}^2 \approx 0{,}728 , U_{cr}^2$$

3. Tranh luận học thuật và định vị của luận án

Một cuộc tranh luận sâu sắc tồn tại giữa hai trường phái:

  • Trường phái giải tích tuyến tính quy chuẩn (TCVN 9901:2014, CIRIA/CUR Rock Manual 1991): Xem vận tốc đáy tại chân kè là một hàm số đối xứng của chiều cao sóng nước sâu và độ sâu nước tĩnh, bỏ qua tương tác phức tạp của cấu trúc độ nhám và dòng chảy ngược.
  • Trường phái mô hình hóa số trị phi tuyến đa pha: Đại diện bởi các nghiên cứu của Van Gent (1995) về phương trình nước nông phi tuyến (NLSW) trong môi trường xốp, và Sakakiyama & Kajima (1992) về mô hình Navier-Stokes mở rộng.

So sánh với các nghiên cứu quốc tế điển hình như nghiên cứu thực địa bãi biển California của Thornton & Guza (1986), nghiên cứu dòng tia sát bờ của Thomas Piter Muray (2004 tại Đại học Griffith, Úc), và nghiên cứu mô hình hóa xói chân kè tại Viện Thủy lực Delft (Hà Lan) của Van Rijn (2018), luận án đã định vị sự đóng góp độc đáo của mình: Không dừng lại ở việc mô phỏng hiện tượng xói bãi cát phẳng thông thường, luận án trực tiếp giải mã cơ chế tương tác sóng - dòng chảy tại vùng chân kè dốc có cấu tạo mố nhám tiêu năng, bắc cầu nối giữa cơ học thủy khí phi tuyến và kỹ thuật công trình thủy thực hành.

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

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

Luận án thách thức trực tiếp tính chuẩn xác của các công thức kinh nghiệm tuyến tính đang sử dụng trong quy chuẩn kỹ thuật quốc gia, cung cấp một bước nhảy vọt về mặt nhận thức lý luận:

  • Mở rộng lý thuyết ổn định hạt rời trên mái dốc của Pilarczyk (1990) và Izbash (1930): Bổ sung các biến số hình học công trình bao gồm độ dốc mái kè ($m$), chiều cao mố nhám ($a$), khoảng cách bố trí mố ($C$) và chỉ số tương tự sóng vỡ Iribarren: $$\xi = \frac{\tan\alpha}{\sqrt{H_s/L_0}}$$
  • Phát triển hệ phương trình toán học mô tả dòng chảy không ổn định trong môi trường xốp hai chiều: Dựa trên mô hình của Sakakiyama và Kajima (1992), kết hợp phương pháp thể tích lưu chất (VOF) và phép lọc không gian rối Smagorinsky (1963), chuẩn hóa hệ số nhớt rối: $$\nu_t = C_s^2 \Delta^2 \sqrt{2 S_{x,z} S_{x,z}}$$ với $C_s = 0{,}2 \div 0{,}25$.
  • Xác lập hệ 3 mệnh đề lý thuyết (Theoretical Propositions):
    • Mệnh đề P1: Vận tốc dòng chảy sát đáy cực đại $U_{max}$ tỷ lệ thuận với căn bậc hai của năng lượng sóng truyền và có quan hệ hàm mũ với chu kỳ sóng $T_p$, thay vì biến thiên tuyến tính đơn thuần theo chiều cao sóng $H_s$.
    • Mệnh đề P2: Góc dốc mái kè $m$ quyết định pha giao thoa giữa sóng tới và sóng phản xạ; mái kè dốc hơn ($m = 2{,}0$ so với $m = 5{,}0$) làm gia tăng ứng suất tiếp đáy cục bộ tại điểm nối tiếp chân đê lên tới 35%.
    • Mệnh đề P3: Cấu trúc mố nhám cản đá hoạt động theo nguyên lý tiêu tán năng lượng dòng xiết và biến đổi trường áp suất cục bộ, hình thành vùng bóng thủy lực làm triệt tiêu lực kéo $F_D$ tác động lên các viên đá bảo vệ chân.

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

Khung phân tích của luận án tích hợp đồng thời 3 trường phái lý thuyết trong một chuỗi tiến trình động học khép kín:

Thành phần khung phân tích Cơ sở lý thuyết nền tảng Phương pháp mô hình hóa Tham số kiểm soát chính
Pha 1: Lan truyền sóng nước nông Sóng phi tuyến Boussinesq & Stokes bậc 2 Máng sóng số MSS-2D & VOF $H_s, T_p, L_s, d_s, U_r$
Pha 2: Tương tác sóng - kết cấu chân kè RANS 2D khép kín rối Smagorinsky Mô hình rỗng Sakakiyama-Kajima Hệ số cản $C_d$, Quán tính $C_m$, Độ xốp $n$
Pha 3: Vận động hạt rời & Mài mòn mái Cân bằng giới hạn lực hạt rời Izbash & Shields Mô hình vật lý tỷ lệ 1:10 (ADV đo đạc) $U_{max}, U_{cr}, d_{50}, \Delta, \Phi, \Psi_c$
Pha 4: Tiêu tán năng lượng mố nhám Thủy lực tiêu năng đáy nhám gia cường Đối chứng thực nghiệm 24 kịch bản Chiều cao mố $a$, Khoảng cách $C$, Hiệu suất $N_o$

Điều kiện biên của khung phân tích được xác định rõ: Áp dụng cho vùng chân kè nông và rất nông ($0{,}3 < d_s / H_{m0,0} < 4{,}0$), bãi biển thoải dốc $i \le 1%$, chịu tác động của phổ sóng ngẫu nhiên JONSWAP hoặc Pierson-Moskowitz với chiều cao sóng thiết kế $H_s \le 3{,}0\text{ m}$.

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

Thiết kế nghiên cứu

Luận án tuân thủ thế giới quan thực chứng luận khoa học (Positivism) kết hợp hiện thực phê phán (Critical Realism), sử dụng thiết kế nghiên cứu hỗn hợp bậc cao (Multi-level Experimental & Numerical Design) kết hợp tương hỗ giữa mô hình số trị và mô hình vật lý.

Quy trình nghiên cứu rigorous

Quy trình thực nghiệm được triển khai nghiêm ngặt qua ba giai đoạn:

Giai đoạn 1: Thiết lập Mô hình Số MSS-2D

Giai đoạn 2: Thực nghiệm Máng sóng Vật lý Tỷ lệ 1:10

Giai đoạn 3: Kiểm chuẩn Chéo & Phân tích Đa biến

Độ tin cậy và tính hợp lệ của nghiên cứu được kiểm soát tuyệt đối: Hệ số phản xạ sóng $K_r$ được kiểm soát qua công nghệ ARC, triệt tiêu sóng phản xạ thứ cấp từ vách máng; độ lặp lại của thí nghiệm đạt mức tin cậy trên 95%.

Data và phân tích

Trong quá trình hiệu chỉnh mô hình số MSS-2D, tác giả đã chạy 16 kịch bản biến thiên đồng thời độ xốp $n$ ($0{,}20 \div 0{,}45$), hệ số gia tăng quán tính $C_m$ ($0{,}20 \div 0{,}80$) và hệ số cản $C_d$ ($0{,}5 \div 2{,}5$) dưới điều kiện sóng $H_s = 0{,}8\text{ m}$, $T_p = 5{,}06\text{ s}$.

"Kết quả hiệu chỉnh mô hình cho thấy với kịch bản độ xốp n = 0,3, Cm = 0,25 và Cd = 1,5 cho vận tốc UmaxTT = 1,26 m/s cho kết quả phù hợp với kết quả thí nghiệm vật lý với cùng điều kiện UmaxTN(m/s) = 1,24 m/s"

Dưới đây là kết quả kiểm chuẩn chéo giữa vận tốc tính toán mô hình số ($U_{maxTT}$) và vận tốc đo đạc thực nghiệm vật lý ($U_{maxTN}$) trên 12 kịch bản mái nghiêng không có mố nhám:

Kịch bản Độ sâu chân kè $d_s$ (m) Chiều cao sóng $H_s$ (m) Chu kỳ sóng $T$ (s) $U_{maxTT}$ (m/s) $U_{maxTN}$ (m/s) Sai số tuyệt đối $\Delta U$ (m/s)
KB1 1,5 0,8 5,06 1,26 1,24 +0,02
KB2 1,5 1,0 6,96 1,84 2,23 -0,39
KB3 1,5 1,1 6,96 2,26 2,35 -0,09
KB4 1,0 0,6 5,06 0,99 1,03 -0,04
KB5 1,0 0,7 5,06 1,15 1,37 -0,22
KB6 1,0 0,8 5,06 1,27 1,46 -0,19
KB7 1,5 0,8 5,69 1,58 1,85 -0,27
KB8 1,5 0,8 6,96 2,03 2,11 -0,08
KB9 1,5 1,0 5,06 1,59 1,66 -0,07
KB10 1,5 1,1 6,32 2,25 2,23 +0,02
KB11 1,0 0,7 5,69 1,41 1,45 -0,04
KB12 1,0 0,8 5,69 1,57 1,62 -0,05

Tập dữ liệu kiểm định khẳng định sự hội tụ chặt chẽ giữa hai phương pháp với hệ số tương quan tuyến tính $R^2 = 0{,}946$, sai số tương đối trung bình dưới $8{,}5%$, khẳng định độ tin cậy vượt trội của mô hình số MSS-2D trong việc mở rộng tính toán cho 320 kịch bản thủy hải văn phức tạp.

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

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

  1. Phát hiện sai số mang tính hệ thống của tiêu chuẩn hiện hành: Tiêu chuẩn TCVN 9901:2014 và công thức Izbash đánh giá thấp vận tốc dòng chảy sát đáy từ $25%$ đến $45%$ trong điều kiện sóng chu kỳ dài ($T_p > 6{,}0\text{ s}$). Khi độ sâu $d_s = 1{,}5\text{ m}$ và $H_s = 1{,}0\text{ m}$, vận tốc thực tế đo được lên tới $2{,}23\text{ m/s}$ (KB2), vượt xa ngưỡng tính toán lý thuyết tuyến tính ($1{,}42\text{ m/s}$).
  2. Xác lập vai trò thống trị của chu kỳ sóng ($T_p$) đối với $U_{max}$: Trong cùng điều kiện độ sâu $d_s = 1{,}5\text{ m}$ và chiều cao sóng $H_s = 0{,}8\text{ m}$, khi chu kỳ $T$ tăng từ $5{,}06\text{ s}$ (KB1) lên $6{,}96\text{ s}$ (KB8), vận tốc $U_{maxTN}$ tăng vọt từ $1{,}24\text{ m/s}$ lên $2{,}11\text{ m/s}$ (tăng trưởng $70{,}16%$).
  3. Làm sáng tỏ quy luật phá hủy mài mòn kết cấu do va đập hạt rời hai chiều: Luận án chứng minh cơ chế phá hủy đê Cát Hải không phải do xói ngầm đơn thuần mà do sự chuyển động lặp của các viên đá hộc chân kè. Vận tốc tới hạn khiến khối đá cấp phối bắt đầu trượt lên mái được xác định chính xác theo công thức: $$U_{cr} = 1{,}171\sqrt{D}$$ Khi sóng rút, dòng undertow kéo đá trượt ngược lại, tạo chu kỳ mài mòn cơ học phá vỡ bề mặt tấm lát bê tông M200/M300.
  4. Định lượng chiều sâu hố xói cân bằng trước chân kè:

    "Chiều sâu hố xói lớn nhất trước chân đê sau 4 giờ chịu tác dụng của sóng (thời gian mô phỏng một cơn bão) đạt giá trị khoảng 1/3, (0,35Hm0) chiều cao sóng trước hố xói"

  5. Hiệu ứng tiêu năng và cản đá đột phá của kết cấu mố nhám: Bố trí hàng mố nhám chiều cao $h = 50\text{ cm}$ làm suy giảm vận tốc ngang cực đại từ $18%$ đến $35%$, đồng thời giữ lại $100%$ các viên đá dịch chuyển tại hàng mố thứ hai.

Implications đa chiều

  • Về mặt lý thuyết: Cung cấp mô hình giải tích mở rộng giải thích trường vận tốc phi đối xứng trong vùng nước cực nông, làm phong phú thêm lý thuyết cơ học môi trường rỗng của Sakakiyama & Kajima.
  • Về mặt phương pháp luận: Chuẩn hóa quy trình phối hợp giữa mô hình số trị MSS-2D (RANS-VOF) và mô hình máng sóng vật lý tỷ lệ 1:10, giảm thiểu $75%$ chi phí thí nghiệm thực tế mà vẫn đảm bảo độ bao quát 320 kịch bản.
  • Về mặt ứng dụng công trình: Đề xuất quy trình 5 bước kiểm tra điều kiện ổn định và chọn kích thước viên đá bảo vệ chân kè, thay thế kết cấu ống buy truyền thống hay bị nứt vỡ bằng kết cấu dầm khóa kết hợp mố nhám cản đá liên hoàn.
  • Về mặt chính sách quản lý: Cung cấp luận cứ khoa học thực chứng để Bộ Nông nghiệp & PTNT, Bộ Xây dựng sửa đổi bổ sung Tiêu chuẩn Quốc gia TCVN 9901:2014 về thiết kế đê biển.

Limitations và Future Research

Luận án thẳng thắn thừa nhận các giới hạn học thuật cần tiếp tục hoàn thiện:

  1. Giới hạn không gian 2D: Nghiên cứu tập trung trên máng sóng phẳng hai chiều (2D), chưa mô phỏng trọn vẹn hiện tượng sóng tới xiên góc ($\beta_0 \ne 0$) và dòng xoáy ba chiều (3D turbulent vortex) hình thành dọc theo đường chân kè uốn khúc.
  2. Giả định hình học hạt rời: Các viên đá thí nghiệm được giả định có dạng tương đối đồng nhất, chưa bao quát đầy đủ sự phân bố cấp phối hạt ngẫu nhiên và hệ số góc cạnh phức tạp của đá hộc khai thác từ các mỏ tự nhiên.
  3. Tương tác động lực với nền cát xói lở: Luận án giới hạn bài toán tương tác trên đáy cố định và bán biến dạng, chưa tích hợp mô hình vận chuyển bùn cát lơ lửng liên tục của nền cát mịn dưới chân thảm đá.

Định hướng nghiên cứu tương lai:

  • Nâng cấp mô hình số trị lên 3D RANS-VOF sử dụng mã nguồn mở OpenFOAM hoặc phương pháp hạt mượt SPH (DualSPHysics).
  • Tích hợp mô hình phần tử rời rạc DEM (Discrete Element Method) liên hợp với CFD để mô phỏng chính xác từng va chạm cơ học giữa viên đá và tấm lát mái đê.
  • Nghiên cứu ứng dụng bê tông tính năng siêu cao (UHPC) và polyme composite có khả năng chống mài mòn chịu lực va đập cực lớn cho hệ mố nhám.

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

+-----------------------------------------------------------------------------------+
|                        MA TRẬN TÁC ĐỘNG TOÀN DIỆN                                 |
+--------------------------+--------------------------------------------------------+
| Lĩnh vực tác động        | Chỉ số Đo lường & Giá trị Thực tiễn                    |
+--------------------------+--------------------------------------------------------+
| 1. Giới Học thuật        | Trở thành tài liệu tham chiếu nòng cốt cho các NCS     |
|                          | ngành Công trình Thủy; dự kiến thu hút trích dẫn cao  |
|                          | trong các nghiên cứu sóng ven bờ tại Đông Nam Á.       |
+--------------------------+--------------------------------------------------------+
| 2. Ngành Xây dựng Thủy lợi| Ứng dụng trực tiếp cho 10 mặt cắt đê Cát Hải; tiết kiệm|
|                          | hàng chục tỷ đồng kinh phí sửa chữa duy tu sau bão.    |
+--------------------------+--------------------------------------------------------+
| 3. Chính sách Nhà nước   | Cung cấp cơ sở khoa học hiệu chỉnh TCVN 9901:2014     |
|                          | và Định hướng quy hoạch đê biển quốc gia đến 2030.     |
+--------------------------+--------------------------------------------------------+
| 4. Lợi ích Xã hội        | Bảo vệ an toàn tuyệt đối cho huyện đảo Cát Hải,        |
|                          | khu kinh tế Đình Vũ và dải dân cư ven biển Hải Phòng.  |
+--------------------------+--------------------------------------------------------+

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

  • Nghiên cứu sinh và Giới học thuật: Tiếp cận cơ sở dữ liệu thực nghiệm 320 kịch bản sóng - cấu trúc rỗng, làm chủ phương pháp luận ghép nối mô hình số Navier-Stokes và máng sóng vật lý hiện đại.
  • Kỹ sư tư vấn thiết kế công trình thủy: Sở hữu công cụ giải tích và quy trình 5 bước để xác định chính xác khối lượng đá đổ $W$, chiều dày thảm đá $D$ và kích thước mố nhám, loại bỏ rủi ro thiết kế thiếu an toàn.
  • Ban Quản lý Dự án và Doanh nghiệp Xây dựng: Nâng cao tuổi thọ công trình đê kè biển từ 15 năm lên trên 30 năm, tối ưu hóa biện pháp thi công lắp ghép mố nhám đúc sẵn ($h = 50\text{ cm}$).
  • Cơ quan Quản lý Nhà nước (Bộ Nông nghiệp & PTNT, Tổng cục Thủy lợi): Có căn cứ pháp lý và tiêu chuẩn kỹ thuật vững chắc để phê duyệt các dự án đầu tư củng cố đê kè biển toàn quốc thích ứng với biến đổi khí hậu.

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

1. Đóng góp lý thuyết độc đáo nhất của luận án là gì và đã mở rộng lý thuyết nào?

Đóng góp độc đáo nhất là việc phát hiện và toán học hóa sự phụ thuộc phi tuyến mạnh mẽ của vận tốc cực đại sát đáy $U_{max}$ vào chu kỳ sóng đỉnh phổ $T_p$, hệ số mái dốc $m$ và độ nhám công trình, mở rộng trực tiếp lý thuyết ổn định hạt rời trên mái dốc của Pilarczyk (1990) và lý thuyết sóng nước nông của Van Rijn (2018). Luận án đã phá vỡ giới hạn của công thức Izbash truyền thống vốn chỉ dựa trên độ sâu nước và chiều cao sóng.

2. Tính đổi mới trong phương pháp nghiên cứu so với các công bố trước đây thể hiện ở điểm nào?

So với các nghiên cứu trước đây (như Wadibe-CT của ĐHTL hay Unibest-TC của Delft chỉ tập trung vào xói bãi phẳng cát mịn), luận án là công trình đầu tiên tại Việt Nam thiết lập quy trình kiểm chuẩn 2 chiều giữa mô hình số trị MSS-2D giải hệ phương trình Navier-Stokes khép kín rối Smagorinsky trong môi trường rỗng (Sakakiyama & Kajima) với mô hình máng sóng vật lý 1:10 có hệ thống hấp thụ sóng phản xạ chủ động ARC và đầu đo lưu tốc âm học ADV tần số cao sát đáy.

3. Phát hiện bất ngờ nhất từ tập dữ liệu thực nghiệm là gì?

Phát hiện bất ngờ nhất là sự gia tăng đột biến của vận tốc đáy khi chu kỳ sóng kéo dài: Tại cùng độ sâu $d_s = 1{,}5\text{ m}$ và $H_s = 0{,}8\text{ m}$, chỉ cần chu kỳ $T$ tăng từ $5{,}06\text{ s}$ lên $6{,}96\text{ s}$, vận tốc đáy $U_{max}$ đo được tăng tới $70{,}16%$ (từ $1{,}24\text{ m/s}$ lên $2{,}11\text{ m/s}$), giải thích nguyên nhân các trận bão có chu kỳ sóng dài phá hủy chân kè khốc liệt dù chiều cao sóng quan trắc không đổi.

4. Luận án có cung cấp quy trình tái lập thí nghiệm (Replication Protocol) chuẩn xác không?

Có. Luận án mô tả chi tiết toàn bộ điều kiện biên: Kích thước máng sóng, tỷ lệ hình học $1:10$, thông số hiệu chỉnh mô hình số ($n = 0{,}3$, $C_m = 0{,}25$, $C_d = 1{,}5$), quy cách đá thí nghiệm ($2\times 4\text{ cm}$ và $1\times 2\text{ cm}$), và 24 kịch bản kiểm định chi tiết từ KB1 đến KB12 kèm thông số $H_s, T_p, d_s$, cho phép bất kỳ phòng thí nghiệm thủy lực nào trên thế giới tái lập hoàn toàn kết quả.

5. Khung chương trình nghiên cứu 10 năm tới được định hình ra sao?

Khung chương trình tập trung vào: (1) Mô hình hóa tương tác 3D sóng xiên - dòng triều - kết cấu chân kè; (2) Ứng dụng phương pháp số hạt rời DEM mô phỏng va đập mài mòn chi tiết; (3) Tích hợp vật liệu bê tông UHPC và hệ thống cảm biến áp lực đáy thông minh giám sát thời gian thực ổn định đê biển quốc gia.

Kết luận

  1. Làm sáng tỏ bản chất thủy động lực học phi tuyến: Chứng minh vận tốc dòng chảy sát đáy tại chân kè nông chịu sự chi phối tổng hợp của $H_s, T_p, d_s, m, i, \xi$, bác bỏ tính đơn giản hóa của các công thức sóng tuyến tính truyền thống trong TCVN 9901:2014.
  2. Khám phá cơ chế phá hủy mài mòn hai chiều do đá hộc dịch chuyển: Xác lập ngưỡng vận tốc tới hạn $U_{cr} = 1{,}171\sqrt{D}$ gây mất ổn định hạt rời chân kè, biến đá hộc thành tác nhân mài mòn phá hủy tấm lát bê tông mái đê.
  3. Phát triển và chuẩn hóa công cụ mô hình số MSS-2D: Tích hợp thành công mô hình môi trường rỗng Navier-Stokes 2 pha với bộ tham số chuẩn xác ($n = 0{,}3, C_m = 0{,}25, C_d = 1{,}5$), đạt độ tương thích cao ($R^2 = 0{,}946$) với thực nghiệm vật lý 1:10.
  4. Sáng chế giải pháp kết cấu mố nhám cản đá tiêu năng: Đột phá kỹ thuật với hệ mố nhám tiêu giảm $18% \div 35%$ vận tốc đáy và cô lập hoàn toàn sự dịch chuyển của đá hộc lên thân đê.
  5. Đề xuất quy trình thiết kế chuẩn hóa áp dụng thực tiễn: Ứng dụng thành công cho tuyến đê biển Cát Hải (Hải Phòng), mang lại giải pháp công trình bền vững, hiệu quả kinh tế - kỹ thuật cao và mở ra hướng nghiên cứu mới cho ngành kỹ thuật công trình biển Việt Nam.