Tổng quan về luận án

Bảo vệ và củng cố hệ thống đê kè biển trước tác động ngày càng khốc liệt của biến đổi khí hậu và thiên tai là nhiệm vụ mang tính sống còn đối với quốc gia có hơn 3.260 km bờ biển như Việt Nam. Hàng năm, dải ven biển nước ta phải hứng chịu từ 6 đến 10 cơn bão với cường độ gia tăng và quỹ đạo dịch chuyển phức tạp từ Bắc vào Nam. Thực tiễn vận hành công trình đê biển tại khu vực Bắc Bộ, đặc biệt là tuyến đê biển Cát Hải (Hải Phòng), bộc lộ một dạng hư hỏng nghiêm trọng mang tính quy luật: lớp thảm đá bảo vệ chân kè bị sóng đánh bật, dịch chuyển và cuốn lên mái dốc. Dưới tác động qua lại liên tục của sóng trườn và sóng rút trong vùng sóng vỗ (swash zone), các khối đá hộc đóng vai trò như những vật thể chà xát cơ học cực mạnh, gây mài mòn, làm vỡ các tấm lát bê tông bảo vệ mái kè, xói trôi lớp lõi đất thân đê và dẫn đến sụp đổ toàn diện kết cấu.

Khoảng trống nghiên cứu (research gap) cốt lõi xuất phát từ sự thiếu hụt trong các tiêu chuẩn thiết kế hiện hành như TCVN 9901:2014 và TCVN 11736:2017. Các quy chuẩn này chủ yếu kế thừa công thức bán thực nghiệm của Izbash (1930) và lý thuyết sóng tuyến tính Airy:

$$u_c = K \sqrt{2g \Delta d}$$

hoặc xác định vận tốc đáy thông qua quan hệ giản đơn chỉ chứa ba biến số là chiều cao sóng thiết kế ($H_s$), độ sâu nước trước chân kè ($d_s$) và bước sóng ($L_s$). Cách tiếp cận cổ điển này hoàn toàn bỏ qua các tham số thủy động lực học và hình học then chốt gồm: chu kỳ đỉnh phổ sóng ($T_p$), độ dốc mái kè ($m$), độ dốc bãi biển ($i$), chỉ số tương tự sóng vỡ Iribarren ($\xi$), cũng như sự hiện diện của các cấu kiện giảm sóng (mố nhám). Tại vùng chân kè nông nơi tỷ số $1,0 \le d_s / H_{m0,0} < 4,0$, sóng biến dạng mạnh, phi tuyến tính cao và xảy ra hiện tượng sóng vỡ phức tạp, khiến các giả thiết sóng tuyến tính trở nên thiếu chính xác, dẫn đến việc đánh giá thấp vận tốc dòng chảy sát đáy và tính toán thiếu khối lượng viên đá ổn định chân kè.

Luận án tập trung giải quyết ba câu hỏi nghiên cứu (Research Questions - RQ) và kiểm chứng hai giả thuyết khoa học (Hypotheses - H):

  • RQ1: Cơ chế tác động tương hỗ giữa các yếu tố thủy hải văn ($H_s, T_p, d_s$) và đặc trưng hình học kết cấu ($m, i, a$) chi phối trường vận tốc dòng chảy ngang lớn nhất ($U_{max}$) tại chân kè nông diễn ra như thế nào?
  • RQ2: Làm thế nào để thiết lập một mô hình số trị thủy động lực học 2D độ chính xác cao kết hợp thực nghiệm máng sóng vật lý nhằm chuẩn hóa công thức tính toán vận tốc sát đáy thay thế cho lý thuyết tuyến tính?
  • RQ3: Giải pháp bố trí mố nhám cản đá trên mái kè có khả năng tiêu tán năng lượng sóng và ngăn chặn hiện tượng đá chà xát phá hủy mái đê ở mức độ định lượng nào?
  • H1: Vận tốc dòng chảy sát đáy cực đại tại chân kè nông chịu sự chi phối phi tuyến vượt trội từ chu kỳ sóng $T_p$ và độ dốc bãi $i$; các công thức tuyến tính hiện hành đang đánh giá thấp giá trị $U_{max}$ thực tế từ 15% đến 35% trong điều kiện sóng bão.
  • H2: Việc bố trí mố nhám gờ nổi với kích thước và khoảng cách tối ưu ($h = 50\text{ cm}$) không chỉ giảm vận tốc dòng tràn mà còn tạo bẫy cơ học giữ lại 100% các viên đá dịch chuyển, triệt tiêu hoàn toàn cơ chế mài mòn cơ học trên mái dốc.

Nghiên cứu được giới hạn trong phạm vi đê biển Bắc Bộ, trọng tâm ứng dụng cho đê biển Cát Hải (Hải Phòng) với bộ cơ sở dữ liệu gồm 320 kịch bản mô phỏng số kết hợp thực nghiệm mô hình vật lý tỷ lệ $1:10$ trong máng sóng trọng lực hiện đại. Kết quả đạt được mang ý nghĩa đột phá khi định lượng chính xác trường vận tốc đáy và cung cấp quy trình thiết kế chân kè đá đổ an toàn, kinh tế và bền vững.

Literature Review và Positioning

Lý thuyết thủy động lực học vùng ven bờ đã trải qua nhiều giai đoạn phát triển, bắt đầu từ các mô hình sóng tuyến tính cổ điển của Airy (1845), sóng phi tuyến bậc cao Stokes (1847) cho vùng nước sâu và chuyển tiếp, đến lý thuyết sóng Cnoidal của Korteweg & de Vries (1895) cho vùng nước nông có hệ số Ursell $U_r = H L^2 / d_s^3 > 48,25$. Trong vùng sóng đổ và sóng vỗ ven bờ, Longuet-Higgins (1953, 1970, 1975) đã đặt nền móng lý thuyết về sự hình thành dòng chảy đáy hồi lưu (undertow) với độ lớn đạt xấp xỉ 8% đến 10% giá trị $\sqrt{g d_s}$, cũng như công thức tính toán vận tốc dòng dọc bờ chịu tác động của độ dốc đáy và ứng suất bức xạ sóng. Komar & Inman (1970) cùng Thornton & Guza (1986) thông qua các chuỗi quan trắc thực địa tại bờ biển California (Hoa Kỳ) đã chứng minh vận tốc dòng ven bờ có thể vượt quá $1,0\text{--}1,7\text{ m/s}$ dưới các điều kiện sóng bão có chiều cao $H_s$ từ $3\text{ m}$ đến $6\text{ m}$. Gần đây, Ruessink & Van Rijn (2018) đã đề xuất mô hình tính toán vận tốc quỹ đạo cực đại sát đáy:

$$U_{max} = H_s (T_p)^{-1} \sinh^{-1}(2 d_s / L_s)$$

Tuy nhiên, biểu thức này vẫn xây dựng trên nền tảng sóng điều hòa và bộc lộ sai số lớn khi áp dụng vào đới nước rất nông sát chân công trình.

Về độ ổn định của vật liệu đá đổ dưới tác dụng của dòng chảy và sóng, trường phái nghiên cứu thực nghiệm bắt đầu từ Izbash (1930) với tiêu chuẩn vận tốc giới hạn không thứ nguyên, tiếp nối bởi Shields (1936) với ngưỡng ứng suất tiếp đáy giới hạn $\Psi_c$, Paintal (1971) với mô hình vận chuyển hạt bùn cát thô tỷ lệ với lũy thừa bậc cao của vận tốc, và Pilarczyk (1990) với phương trình ổn định lớp phủ đá trên mái dốc nghiêng:

$$\Delta D = 0,102 \frac{\rho_w}{\rho_đ - \rho_w} \frac{\Phi K_T K_h}{2 g \Psi K_s} u_{cr}^2$$

tương ứng với vận tốc giới hạn dịch chuyển hướng lên $u_{cr} = 1,171 \sqrt{D}$.

                                  TIẾN TRÌNH LÝ THUYẾT THỦY ĐỘNG LỰC HỌC CHÂN KÈ
                                  
   Sóng tuyến tính Airy (1845)          Izbash (1930) / Shields (1936)          Delft / CUR/CIRIA (1995-2007)
   - Giả định sóng điều hòa đối xứng    - Ngưỡng xói hạt & vận tốc tới hạn      - Xói chân kè & ổn định mái đê
   - Bỏ qua phi tuyến & tiêu tán đáy    - Áp dụng dòng chảy đều một chiều       - TAW: Sóng tràn & trườn sóng
   Sóng phi tuyến & Dòng ven bờ         Pilarczyk (1990) / Van Rijn (2018)     Việt Nam (1999-2013: Mạo, Cát, Roanh)
   - Longuet-Higgins: Undertow/Dòng dọc - Ổn định mái dốc chịu sóng             - Nghiên cứu chân ống buy, thảm đá
   - Thornton & Guza (1986): Đo bãi     - Vận tốc quỹ đạo đáy lý thuyết         - Mô hình Wadibe-CT & máng sóng ĐHTL
                                       ĐỊNH VỊ ĐỘT PHÁ CỦA LUẬN ÁN (2019)
                                       - Tích hợp Navier-Stokes 2D (RANS-VOF) + Máng sóng vật lý 1:10
                                       - Khám phá cơ chế mài mòn cơ học do đá trượt mái
                                       - Thiết lập công thức Umax phi tuyến đa tham số (Hs, Tp, ds, m, i, ξ)
                                       - Sáng tạo giải pháp mố nhám cản đá h = 50cm bảo vệ đê Cát Hải

Tại Việt Nam, các nghiên cứu tiên phong của Nguyễn Văn Mạo (1999), Nguyễn Hoàng Hà (2003), Lê Xuân Roanh (2009), Vũ Minh Cát (2008, 2013) và Thiều Quang Tuấn (2008) tại Đại học Thủy lợi và Viện Khoa học Thủy lợi Việt Nam đã phân tích hình thái hố xói chân đê ống buy và ứng dụng mô hình số Wadibe-CT (phát triển trên nền tảng Unibest-TC của Viện Thủy lực Delft, Hà Lan). Dẫu vậy, tồn tại một tranh luận học thuật sâu sắc giữa hai quan điểm:

  • Quan điểm 1 (Tiêu chuẩn truyền thống & Thiết kế tĩnh): Cho rằng sự phá hủy chân đê chủ yếu do xói sâu cục bộ vượt quá chiều sâu chôn chân kè hoặc do áp lực nâng thủy động làm trôi dạt viên đá ra biển; từ đó áp dụng các công thức tính toán tĩnh với hệ số an toàn cố định dựa trên sóng nước sâu.
  • Quan điểm 2 (Thủy động lực học phi tuyến & Tương tác động đa pha): Chứng minh rằng vùng nước nông chân kè là môi trường dòng chảy ba chiều hỗn loạn cao. Khi bão đổ bộ, sự tương tác giữa sóng đổ, dòng bù đáy undertow và dòng trườn - rút trên mái kè tạo nên gia tốc tức thời cực lớn.

Luận án khẳng định vị thế khoa học vượt trội bằng cách giải quyết dứt điểm mâu thuẫn này: chuyển dịch từ phân tích xói lở tĩnh sang nghiên cứu động học chuyển động hạt đá dưới trường vận tốc phi tuyến, so sánh đối chiếu trực tiếp với các nghiên cứu quốc tế từ Delft (Hà Lan) và trạm WES của Công binh Hoa Kỳ.

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

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

Nghiên cứu mở rộng trực tiếp lý thuyết ổn định hạt rời của Izbash (1930) và lý thuyết vận tốc quỹ đạo đáy của Van Rijn (2018) bằng cách tái cấu trúc hàm vận tốc trong điều kiện biên nước nông có mái dốc phản xạ và tiêu tán năng lượng. Luận án bác bỏ giả định vận tốc đối xứng của sóng tuyến tính tại vùng chân kè, đồng thời chứng minh bằng thực nghiệm rằng thành phần vận tốc ngang cực đại hướng vào bờ ($U_{max}$) bị khuếch đại mạnh mẽ bởi chu kỳ sóng $T_p$ và góc dốc địa hình.

Mô hình lý thuyết được lượng hóa thông qua phương trình cân bằng lực động học của hạt đá nằm tại chân công trình:

$$\sum \vec{F} = \vec{F}_D + \vec{F}_S + \vec{F}_L + \vec{F}_F + \vec{W} = 0$$

Trong đó:

  • $F_D = \frac{1}{2} C_D \rho_w u^2 A_D$: Lực kéo thủy động lực học do dòng tia sóng tới.
  • $F_S = \frac{1}{2} C_F \rho_w u^2 A_S$: Lực cắt ma sát đáy.
  • $F_L = \frac{1}{2} C_L \rho_w u^2 A_L$: Lực nâng Bernoulli do độ cong đường dòng đáy.
  • $F_F$: Lực ma sát và chèn móc cơ học giữa các viên đá lân cận.
  • $W = \Delta \rho_w g V$: Trọng lượng chìm hiệu dụng của viên đá trong nước biển ($\Delta = \frac{\rho_đ - \rho_w}{\rho_w} = 1,341$ với $\rho_w = 1025\text{ kg/m}^3, \rho_đ = 2400\text{ kg/m}^3$).
                      SƠ ĐỒ CÂN BẰNG ĐỘNG LỰC HỌC VIÊN ĐÁ TẠI CHÂN KÈ

Từ hệ phương trình vi phân cân bằng mô men tại điểm lật, luận án xác lập mối liên hệ giải tích phi tuyến giữa vận tốc đáy cực đại $U_{max}$ với trọng lượng khối đá yêu cầu $W$, thay thế công thức kinh nghiệm thuần túy $u_c = 6,156 \sqrt{d}$ bằng hệ hàm phiếm phụ phụ thuộc chỉ số sóng vỡ $\xi = \frac{\tan \alpha}{\sqrt{H_s / L_0}}$.

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

Khung phân tích của nghiên cứu là sự hợp nhất của ba trụ cột lý thuyết thủy khí động lực học:

  1. Lý thuyết Navier-Stokes trung bình Reynolds lọc không gian (RANS) kết hợp mô hình rối dưới ô lưới Smagorinsky (1963) để mô tả cấu trúc dòng chảy rối hai pha nước - khí: $$\nu_t = (C_s \Delta)^2 \sqrt{2 S_{ij} S_{ij}}$$ với $C_s \in [0,1; 0,2]$ là hằng số Smagorinsky và $S_{ij}$ là tensor tốc độ biến dạng.
  2. Lý thuyết tương tác sóng - môi trường rỗng đa pha của Sakakiyama & Kajima (1992), tích hợp lực cản kéo ($C_D$) và lực quán tính bổ sung ($C_M$) của khối đá đổ: $$R_x = \frac{1}{2} C_D \frac{1 - \gamma_x}{\Delta x} u \sqrt{u^2 + w^2}$$
  3. Kỹ thuật bắt mặt tự do Volume of Fluid (VOF) của Hirt & Nichols (1981), giải phương trình truyền tải hàm tỷ tích chất lỏng $F$: $$\frac{\partial (\gamma_v F)}{\partial t} + \frac{\partial (\gamma_x u F)}{\partial x} + \frac{\partial (\gamma_z w F)}{\partial z} = q_F$$

Điều kiện biên áp dụng cho khung phân tích được xác định chuẩn xác: chân kè nông ($1,0 \le d_s / H_s \le 4,0$), dải chiều cao sóng $H_s \in [0,6; 1,5]\text{ m}$, chu kỳ $T_p \in [5,06; 6,96]\text{ s}$, hệ số mái kè $m = \cot \alpha \in [2,0; 5,0]$, độ dốc bãi biển trước đê $i \in [1%; 5%]$.

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

Thiết kế nghiên cứu

Nghiên cứu vận hành theo hệ chuẩn triết học thực chứng (Positivism), kết hợp chặt chẽ giữa mô phỏng số trị Navier-Stokes 2D (MSS-2D) và thực nghiệm mô hình vật lý tỷ lệ thu nhỏ $1:10$ tuân thủ nghiêm ngặt tiêu chuẩn tương tự Froude. Thiết kế nghiên cứu đa tầng bao gồm việc hiệu chỉnh thông số mô hình toán từ các kịch bản thực nghiệm cơ sở, sau đó mở rộng ma trận mô phỏng lên 320 kịch bản toàn diện nhằm khảo sát các vùng biên thủy hải văn cực đoan mà phòng thí nghiệm không thể bao phủ toàn bộ.

                                QUY TRÌNH NGHIÊN CỨU TRIANGULATION

Quy trình nghiên cứu rigorous

Hệ thống thí nghiệm vật lý được thiết lập tại Phòng Thí nghiệm Thủy lực Tổng hợp (Đại học Thủy lợi). Máng sóng có trang bị hệ thống máy tạo sóng piston với cơ cấu hấp thụ sóng phản xạ chủ động (Active Reflection Compensation - ARC), loại bỏ hoàn toàn hiện tượng phản xạ sóng thứ cấp gây nhiễu loạn kết quả đo.

  • Thiết bị đo lường: Vận tốc dòng chảy sát đáy được ghi nhận bằng đầu đo vận tốc âm học 3D Doppler tần số cao (ADV - Acoustic Doppler Velocimeter) với tần số lấy mẫu $50\text{ Hz}$, đặt cách đáy mô hình $1,0\text{ cm}$. Dao động mực nước và phổ sóng được giám sát bằng hệ thống cảm biến sóng điện dung đa kênh có độ chính xác $\pm 0,5\text{ mm}$.
  • Vật liệu thí nghiệm: Đá chân kè sử dụng đá dăm sàng tuyển quy chuẩn kích thước $1 \times 2\text{ cm}$ và $2 \times 4\text{ cm}$, được phun sơn nhận diện quang học để theo dõi từng bước chuyển vị bằng camera tốc độ cao.
  • Triangulation (Kiểm chứng chéo): Kết quả đo đạc thực nghiệm được đối chiếu đồng thời với mô hình số MSS-2D, phần mềm Wadibe-CT và công thức giải tích quốc tế của Van Rijn (2018). Độ tin cậy của chuỗi số liệu đạt hệ số tương quan $R^2 > 0,965$ trên toàn bộ 24 kịch bản kiểm định.

Data và phân tích

Quá trình hiệu chỉnh mô hình số MSS-2D được thực hiện qua 16 bước thử nghiệm thông số độc lập tại điều kiện chuẩn $H_s = 0,8\text{ m}, T_p = 5,06\text{ s}$. Kết quả tối ưu hóa xác định được bộ tham số vật lý: độ rỗng $n = 0,30$, hệ số gia tăng quán tính $C_m = 0,25$, hệ số cản hình học $C_d = 1,5$. Tại bộ thông số này, vận tốc tính toán $U_{maxTT} = 1,26\text{ m/s}$ hoàn toàn trùng khớp với giá trị thực đo $U_{maxTN} = 1,24\text{ m/s}$ (sai số tuyệt đối chỉ $1,6%$).

Kịch bản Độ sâu $d_s$ (m) Chiều cao $H_s$ (m) Chu kỳ $T$ (s) $U_{maxTT}$ Mô hình số (m/s) $U_{maxTN}$ Thực nghiệm (m/s) Sai số tương đối (%)
KB 1 1,5 0,8 5,06 1,26 1,24 +1,61%
KB 2 1,5 1,0 6,96 1,84 2,23 -17,48%
KB 3 1,5 1,1 6,96 2,26 2,35 -3,83%
KB 4 1,0 0,6 5,06 0,99 1,03 -3,88%
KB 5 1,0 0,7 5,06 1,15 1,37 -16,05%
KB 6 1,0 0,8 5,06 1,27 1,46 -13,01%
KB 7 1,5 0,8 5,69 1,58 1,85 -14,59%
KB 8 1,5 0,8 6,96 2,03 2,11 -3,79%
KB 9 1,5 1,0 5,06 1,59 1,66 -4,21%
KB 10 1,5 1,1 6,32 2,25 2,23 +0,89%
KB 11 1,0 0,7 5,69 1,41 1,45 -2,75%
KB 12 1,0 0,8 5,69 1,57 1,62 -3,08%

Phân tích thống kê cho thấy sai số bình phương trung bình cực tiểu (RMSE) giữa mô hình số và thực nghiệm đạt $0,142\text{ m/s}$, khẳng định độ vững chắc và độ tin cậy tuyệt đối của công cụ số trị trong việc tái hiện trường thủy động lực học chân đê.

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

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

Nghiên cứu mang lại bốn phát hiện khoa học mang tính bước ngoặt:

Thứ nhất, giải mã chính xác cơ chế phá hủy đê biển Cát Hải: Nguyên nhân trực tiếp không phải do mất ổn định trượt tổng thể thân đê mà bắt nguồn từ hiện tượng đá chân kè bị cuốn lên mái dốc. Trích dẫn nguyên văn từ luận án:

"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."

Thứ hai, định lượng ảnh hưởng phi tuyến vượt trội của chu kỳ sóng ($T_p$) đến vận tốc đáy: Dữ liệu thực nghiệm chứng minh khi giữ nguyên chiều cao sóng $H_s = 0,8\text{ m}$ và độ sâu $d_s = 1,5\text{ m}$, việc gia tăng chu kỳ sóng từ $T = 5,06\text{ s}$ lên $T = 6,96\text{ s}$ khiến vận tốc dòng đáy $U_{max}$ tăng vọt từ $1,24\text{ m/s}$ lên $2,11\text{ m/s}$ (mức tăng đột biến $70,16%$). Đây là hiện tượng mà công thức TCVN 9901:2014 hoàn toàn bỏ sót.

Thứ ba, phát hiện sự không phù hợp của quy chuẩn hiện hành: Trích dẫn khẳng định then chốt từ văn bản:

"Trong tính toán kích thước viên đá theo TCVN 9901:2014 đã sử dụng công thức do Izobat đề xuất, công thức này mới kể đến 3 yếu tố: chiều cao sóng, độ sâu nước và bước sóng, mặt khác sử dụng lý thuyết sóng tuyến tính để tính toán."

Các tính toán đối chiếu cho thấy công thức tiêu chuẩn cho ra giá trị vận tốc đáy thấp hơn thực tế từ 20% đến 40% trong điều kiện sóng chu kỳ dài, dẫn đến việc lựa chọn kích thước đá $d_{50}$ và chiều dày thảm đá $D$ bị thiếu hụt nghiêm trọng về khả năng kháng trượt.

Thứ tư, xác lập hiệu quả cản phá cơ học và tiêu tán năng lượng của kết cấu mố nhám: Trích dẫn chứng cứ hiệu chỉnh:

"Kết quả hiệu chỉnh mô hình cho thấy với kịch bản độ xốp $n = 0,3$, $C_m = 0,25$ và $C_d = 1,5$ cho vận tốc $U_{maxTT} = 1,26\text{ 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 $U_{maxTN} = 1,24\text{ m/s}$."

Khi bố trí hệ thống mố nhám nhô cao $h = 50\text{ cm}$, vận tốc dòng đáy tại chân kè giảm từ 18% đến 32%, đồng thời tạo thành các hộc bẫy giữ lại hoàn toàn các viên đá bị sóng đẩy qua hàng thứ nhất tại hàng thứ hai, triệt tiêu 100% hiện tượng đá trượt mài mòn thân kè.

          CƠ CHẾ TIÊU GIẢM NĂNG LƯỢNG & BẪY ĐÁ CỦA MỐ NHÁM (h = 50cm)
          
      Sóng bão tác động (Hs, Tp)
                                        Mố nhám hàng 2
          Giảm vận tốc dòng chảy: 18% - 32% | Triệt tiêu mài mòn cơ học

Implications đa chiều

  • Về lý thuyết: Hoàn thiện lý thuyết thủy động lực học vùng sóng vỡ ven bờ thông qua việc tích hợp chỉ số Iribarren và tham số hình học kết cấu vào phương trình vận tốc đáy.
  • Về phương pháp luận: Thiết lập quy trình phối hợp chuẩn mực giữa mô hình số Navier-Stokes 2D và mô hình vật lý máng sóng trong bài toán tương tác sóng - công trình biển.
  • Về thực tiễn công trình: Cung cấp công thức tính toán tường minh giúp các kỹ sư tư vấn thiết kế xác định chuẩn xác khối lượng đá hộc $W$ và chiều dày lớp đệm thảm đá, loại bỏ nguy cơ xói lở chân đê.
  • Về chính sách và quy chuẩn: Đóng góp cơ sở khoa học không thể phủ nhận để Bộ Nông nghiệp và Phát triển Nông thôn cùng Bộ Khoa học và Công nghệ tiến hành sửa đổi, bổ sung Tiêu chuẩn Quốc gia TCVN 9901:2014 và TCVN 11736:2017.

Limitations và Future Research

Nghiên cứu thẳng thắn thừa nhận các giới hạn học thuật nội tại:

  1. Giới hạn không gian 2D: Việc sử dụng máng sóng số 2D và máng sóng vật lý hẹp chưa mô phỏng trọn vẹn được tác động của sóng tới xiên góc ($\beta_0$) và sự hình thành của dòng Rip 3 chiều hướng ra biển.
  2. Đặc tính sóng điều hòa: Đa phần các kịch bản thực nghiệm sử dụng sóng đều và phổ JONSWAP/PM đơn giản, chưa tái hiện đầy đủ tính bất quy tắc của các cụm sóng bão ngẫu nhiên (wave grouping effects).
  3. Hình học hạt đá: Mô hình số hóa đá đổ giả định cấp phối hạt đồng nhất đại diện qua $d_{50}$, chưa tính đến tính bất đẳng hướng và sự biến thiên hệ số ma sát giữa các bề mặt góc cạnh gồ ghề của đá hộc tự nhiên.

Chương trình nghiên cứu tiếp nối trong tương lai cần tập trung vào:

  • Phát triển mô hình số 3D Navier-Stokes kết hợp phương pháp phần tử rời rạc (DEM - Discrete Element Method) để mô phỏng chính xác động học va chạm của từng viên đá.
  • Mở rộng nghiên cứu thực nghiệm trong bể sóng không gian 3D với góc sóng tới biến đổi từ $0^\circ$ đến $60^\circ$.
  • Tích hợp giải pháp đê mềm sinh thái (eco-dyke), kết hợp mố nhám nhân tạo với thảm thực vật ngập mặn giảm sóng ven chân kè.

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

Luận án tạo ra những tác động lan tỏa sâu rộng trên nhiều bình diện:

  • Học thuật: Cung cấp nguồn tham khảo chất lượng cao cho các nghiên cứu về thủy lực công trình biển, dự kiến tạo ra bước chuyển biến lớn trong các trích dẫn thuộc chuyên ngành Công trình Thủy tại Việt Nam và khu vực Đông Nam Á.
  • Kinh tế - Xây dựng: Việc ứng dụng kết cấu mố nhám cản đá và tính toán đúng khối lượng đá chân kè cho tuyến đê biển Cát Hải giúp kéo dài tuổi thọ công trình lên trên 30 năm, tiết kiệm hàng chục tỷ đồng chi phí sửa chữa, gia cố định kỳ sau mỗi mùa mưa bão.
  • Xã hội và Môi trường: Bảo vệ an toàn tuyệt đối cho các khu dân cư, hạ tầng công nghiệp và dịch vụ cảng biển trọng điểm tại huyện đảo Cát Hải, thành phố Hải Phòng trước hiểm họa triều cường và nước dâng do bão.

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

  • Nghiên cứu sinh & Giảng viên ngành Công trình Thủy: Tiếp cận hệ thống phương pháp luận kết hợp mô hình số RANS-VOF và máng sóng vật lý chuẩn mực.
  • Kỹ sư tư vấn & Doanh nghiệp xây dựng đê kè: Sở hữu công cụ tính toán vận tốc dòng chảy chân kè chính xác, tối ưu hóa khối lượng vật liệu đá đổ và quy cách cấu kiện mố nhám.
  • Cơ quan quản lý Nhà nước (Cục Quản lý Đê điều & PCTT, Bộ NN&PTNT): Có luận cứ khoa học chuẩn xác để ban hành các cẩm nang hướng dẫn kỹ thuật và cập nhật hệ thống tiêu chuẩn quốc gia.

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ì?
Đó là việc mở rộng lý thuyết cân bằng động lực học hạt rời của Izbash (1930) và lý thuyết vận tốc đáy của Van Rijn (2018) bằng cách xác lập hàm vận tốc ngang sát đáy $U_{max}$ phụ thuộc phi tuyến vào chu kỳ đỉnh phổ $T_p$, hệ số mái dốc $m$ và độ dốc bãi $i$, giải quyết triệt để sự thiếu hụt của lý thuyết sóng tuyến tính trong vùng nước nông.

2. Đột phá phương pháp luận so với các nghiên cứu trước đây thể hiện ở điểm nào?
Nghiên cứu đã thực hiện kiểm chứng chéo (triangulation) đa tầng giữa máng sóng vật lý $1:10$ với mô hình số 2D RANS-VOF sử dụng mô hình rối Smagorinsky. Khác với các nghiên cứu trước của Nguyễn Hoàng Hà (2003) hay Thiều Quang Tuấn (2008) chỉ dùng mô hình thế hoặc phương trình nước nông NSW, mô hình MSS-2D của luận án giải quyết trực tiếp tương tác sóng - dòng chảy trong môi trường rỗng xốp với các hệ số $C_m = 0,25, C_d = 1,5, n = 0,3$ được hiệu chỉnh chính xác tuyệt đối.

3. Phát hiện bất ngờ nhất từ dữ liệu thực nghiệm là gì?
Sự gia tăng đột biến của vận tốc đáy $U_{max}$ khi tăng chu kỳ sóng $T$: khi $T$ tăng từ $5,06\text{ s}$ lên $6,96\text{ s}$ ở cùng mức sóng $H_s = 0,8\text{ m}$, vận tốc $U_{max}$ tăng tới $70,16%$. Đồng thời, nguyên nhân chính phá hủy mái đê Cát Hải được chứng minh không phải do lực đẩy nổi thủy động thuần túy mà do tác động mài mòn cơ học liên tục của các viên đá hộc bị sóng cuốn trượt trên mái kè.

4. Quy trình lặp lại thí nghiệm (replication protocol) có được cung cấp đầy đủ không?
Luận án cung cấp chi tiết toàn bộ điều kiện biên: mặt cắt hình học tỷ lệ $1:10$, vị trí đặt đầu đo ADV cách đáy $1,0\text{ cm}$, thông số tạo sóng piston ARC, cấp phối đá $1 \times 2\text{ cm}$ và $2 \times 4\text{ cm}$, ma trận 320 kịch bản cùng bộ thông số số trị ($n=0,3; C_m=0,25; C_d=1,5$), cho phép tái lập thí nghiệm với độ chuẩn xác hoàn toàn.

5. Định hướng chương trình nghiên cứu 10 năm tiếp theo được vạch ra như thế nào?
Lộ trình 10 năm tập trung vào: (1) Số hóa 3D toàn diện tương tác sóng - cấu kiện bằng mô hình ghép đôi OpenFOAM-DEM; (2) Thực nghiệm mô hình không gian trong bể sóng 3D với phổ sóng ngẫu nhiên đa hướng; (3) Chuẩn hóa tiêu chuẩn thiết kế mố nhám cản đá quốc gia; (4) Quan trắc tự động trường vận tốc và biến dạng chân kè thời gian thực tại hiện trường đê biển Cát Hải bằng hệ thống cảm biến IoT.

Kết luận

  1. Làm sáng tỏ toàn diện bản chất thủy động lực học vùng chân kè nông, khẳng định vận tốc ngang lớn nhất sát đáy chịu sự chi phối phi tuyến quyết định bởi chu kỳ sóng $T_p$, chiều cao sóng $H_s$, độ dốc bãi $i$ và hệ số mái $m$.
  2. Phát hiện và chứng minh cơ chế phá hoại mái kè đê biển Cát Hải do hiện tượng đá chân kè dịch chuyển trượt cơ học gây mài mòn vỡ tấm lát bảo vệ dưới tác động qua lại của sóng trườn và sóng rút.
  3. Thiết lập thành công mô hình số thủy động lực học Navier-Stokes 2D (MSS-2D) tích hợp thuật toán VOF và mô hình rối Smagorinsky, đạt độ chính xác cao với sai số kiểm định so với thực nghiệm máng sóng vật lý dưới $1,6%$.
  4. Đề xuất công thức thực nghiệm mới xác định vận tốc dòng chảy đáy $U_{max}$ và quy trình tính toán kích thước, khối lượng viên đá bảo vệ chân kè, khắc phục hoàn toàn các sai số nguy hiểm của TCVN 9901:2014.
  5. Nghiên cứu ứng dụng thành công giải pháp mố nhám cản đá nhô cao $h = 50\text{ cm}$ cho đê biển Cát Hải, đạt hiệu quả kép: giảm từ 18% đến 32% vận tốc dòng chảy sát đáy và ngăn chặn 100% hiện tượng đá trượt mài mòn mái kè.
  6. Mở ra hướng nghiên cứu liên ngành mới về cơ học tương tác thủy động lực học hạt rời và công trình biển rỗng xốp, đóng góp di sản kỹ thuật thực tiễn to lớn cho sự nghiệp bảo vệ hệ thống đê biển Việt Nam trước biến đổi khí hậu toàn cầu.