Tổng quan về luận án

Nghiên cứu về công trình bảo vệ bờ biển và cảng nước sâu luôn giữ vai trò then chốt trong chiến lược phát triển kinh tế biển và củng cố an ninh quốc phòng. Trong hệ thống kết cấu bảo vệ, đê chắn sóng đá đổ mái nghiêng phủ khối bê tông dị hình đóng vai trò triệt tiêu năng lượng sóng, bảo đảm an toàn cho luồng tàu và khu nước bể cảng. Luận án tiến sĩ kỹ thuật của tác giả Lê Thị Hương Giang với đề tài "Nghiên cứu đánh giá ổn định và hiệu quả tiêu giảm sóng tràn của khối phủ RAKUNA IV cho đê chắn sóng đá đổ mái nghiêng" (Chuyên ngành Xây dựng công trình thủy, Mã số: 62-58-40-01, Trường Đại học Thủy lợi, năm 2015, người hướng dẫn: PGS. Thiều Quang Tuấn và GS. Hồ Sĩ Minh) đại diện cho bước phát triển học thuật tiên phong trong lĩnh vực kỹ thuật công trình biển tại Việt Nam.

Nghiên cứu xuất phát từ một khoảng trống học thuật cụ thể (research gap): Cấu kiện tiêu sóng RAKUNA IV là loại khối phủ dị hình thế hệ mới do Nhật Bản phát minh năm 2007 với kết cấu 4 chân góc cạnh và 4 hốc lõm đặc thù. Tuy nhiên, các công bố quốc tế trước đó như Mase và cộng sự (2011) hay Tuấn và cộng sự (2012) mới chỉ đánh giá khối phủ RAKUNA IV trên kết cấu đê hỗn hợp ngang hoặc đê mái nghiêng trong điều kiện không có sóng tràn (zero-overtopping). Thực tế chỉ ra: "chưa có công trình nghiên cứu nào về đánh giá ổn định và hiệu quả tiêu giảm sóng tràn của khối phủ RAKUNA IV cho đê chắn sóng đá đổ mái nghiêng". Khi triết lý thiết kế đê hiện đại chuyển dịch từ ngăn chặn hoàn toàn sóng leo sang cho phép sóng tràn có kiểm soát để giảm cao trình đỉnh đê, cơ chế gia tăng ổn định thủy lực và định lượng hệ số chiết giảm sóng tràn của RAKUNA IV trở thành bài toán chưa có lời giải thực nghiệm lẫn mô hình toán.

Luận án tập trung giải quyết ba câu hỏi nghiên cứu cốt lõi:

  1. RQ1: Sóng tràn tác động như thế nào đến độ ổn định thủy lực của lớp khối phủ RAKUNA IV xếp hai lớp trên mái đê phía biển so với điều kiện không tràn?
  2. RQ2: Hệ số chiết giảm sóng tràn $\gamma_r$ của RAKUNA IV phụ thuộc ra sao vào thông số độ dốc sóng, số Iribarren $\xi_{m-1,0}$ và chiều cao lưu không tương đối $R_c/H_{m0}$?
  3. RQ3: Bản chất thủy động lực học của hiệu ứng "đệm nước" (water cushion effect) bên trong lớp phủ rỗng chi phối cơ chế tiêu tán năng lượng sóng như thế nào?

Ba giả thuyết nghiên cứu tương ứng được thiết lập:

  • H1: Hiện tượng sóng tràn làm chuyển dịch một phần đáng kể năng lượng sóng qua đỉnh và mái sau, giúp gia tăng độ ổn định thủy lực của lớp phủ RAKUNA IV ở mái trước thông qua một hệ số hàm số mũ $F_s > 1.0$ phụ thuộc vào tỷ số $R_c/H_{m0}$.
  • H2: Hệ số chiết giảm sóng tràn $\gamma_r$ của RAKUNA IV không phải là một hằng số cố định mà biến thiên tuyến tính tăng dần theo số Iribarren $\xi_{m-1,0}$, phản ánh sự suy giảm hiệu quả tiêu sóng khi chu kỳ sóng dài ra.
  • H3: Dòng chảy thấm ngập lưu giữ nước trong các hốc lõm và độ rỗng $n = 56,5%$ tạo thành lớp đệm chất lưu triệt tiêu độ nhám thủy lực biểu kiến đối với các con sóng kế tiếp.

Khung lý thuyết của luận án tích hợp lý thuyết cân bằng động lực học Morison, lý thuyết phân tích thứ nguyên $\Pi$-Buckingham, mô hình ổn định thực nghiệm Van der Meer (1988), De Jong (1996) và lý thuyết sóng tràn TAW (2002), EurOtop (2007). Đóng góp mang tính đột phá của công trình là xác lập công thức thực nghiệm tính toán hệ số gia tăng ổn định $F_s$, xây dựng biểu thức tương quan hệ số $\gamma_r$ và mô phỏng thành công trường dòng chảy hai pha khí - lỏng bằng mô hình số IH-2VOF. Nghiên cứu được kiểm chứng với quy mô thực nghiệm 58 chuỗi thí nghiệm sóng ngẫu nhiên JONSWAP trong máng sóng vật lý 2D chiều dài 45 m, kết hợp áp dụng tính toán thực tế cho đê chắn sóng cảng Nghi Sơn (Thanh Hóa), chứng minh khả năng cắt giảm khối lượng bê tông và chi phí xây lắp đáng kể so với phương án khối Tetrapod truyền thống.

Literature Review và Positioning

Lịch sử phát triển của kết cấu khối phủ đê chắn sóng đá đổ mái nghiêng ghi nhận sự tiến hóa từ các giải pháp tự trọng sang giải pháp khai thác lực liên kết cài khóa (interlocking). Trước thập niên 1950, đê chắn sóng chủ yếu sử dụng đá hộc hoặc khối bê tông trọng lực hình hộp, phụ thuộc hoàn toàn vào trọng lượng bản thân ($K_D \le 4,0$), dẫn đến kích thước kết cấu khổng lồ và chi phí xây dựng đắt đỏ. Bước ngoặt xuất hiện vào năm 1950 khi phòng thí nghiệm Dauphinois (SOGREAH) phát minh ra khối Tetrapod. Sau đó, hàng loạt khối phủ dị hình ra đời: Tribar (1958), Stabit (1961), Dolos (1963), Accropode (1981), Core-Loc (1994), X-block (2003) và RAKUNA IV (2007).

Các nghiên cứu kinh điển về ổn định thủy lực khởi nguồn từ Iribarren (1938) với mô hình cân bằng tĩnh lực, tiếp nối bởi công thức bán thực nghiệm của Hudson (1959) được quy chuẩn trong SPM (1977, 1984):

$$N_s = \frac{H_{1/3}}{\Delta D_n} = (K_D \cot \alpha)^{1/3}$$

Tuy nhiên, hạn chế lớn nhất của công thức Hudson là bỏ qua chu kỳ sóng ($T_m$), tính thấm của lõi đê ($P$) và thời gian tác động của bão bùng phát qua số con sóng ($N_z$). Để khắc phục, Van der Meer (1988) đã thiết lập hệ thống công thức thực nghiệm cho khối phủ dị hình dựa trên số khối dịch chuyển tương đối $N_{od}$ và độ dốc sóng $s_{om}$:

$$N_s = \frac{H_s}{\Delta D_n} = \left( 3,75 \frac{N_{od}^{0,5}}{N_z^{0,25}} + 0,85 \right) s_{om}^{-0,2}$$

Đối với bài toán khối phủ chịu điều kiện sóng tràn, De Jong (1996) đã tiến hành thí nghiệm trên khối Tetrapod 2 lớp và đưa ra hệ số gia tăng ổn định:

$$f(R_c/D_n) = 1 + 0,17 \exp\left(-0,61 \frac{R_c}{D_n}\right)$$

Trong mảng nghiên cứu sóng tràn, các mô hình của Owen (1980, 1982), Van der Meer (1999), TAW (2002) và EurOtop (2007) biểu diễn lưu lượng tràn đơn vị $q$ thông qua hàm số mũ suy giảm theo chiều cao lưu không tương đối $R_c/H_{m0}$ và hệ số chiết giảm độ nhám $\gamma_r$. Bảng giá trị chuẩn của EurOtop (2007) ghi nhận $\gamma_r$ của Tetrapod là 0,38; Accropode là 0,46; Dolos là 0,43.

Tiến trình học thuật & Định vị nghiên cứu:
[Hudson (1959) / SPM (1977, 1984)] -> [Van der Meer (1988): Ns = f(Nod, Nz, som)]
[De Jong (1996): Ổn định Tetrapod khi tràn]                                  [Bruce et al. (2009): γr biến thiên theo ξm-1,0]
                       [Luận án: RAKUNA IV - Ổn định & Tiêu giảm sóng tràn trên đê mái nghiêng]
                       • Thiết lập hàm gia tăng ổn định Fs = f(Rc/Hm0)
                       • Xác lập γr = f(ξm-1,0) cho RAKUNA IV (n = 56,5%)
                       • Lý giải bản chất cơ học của hiệu ứng "Đệm nước" bằng IH-2VOF

Trong bức tranh tổng quan đó, tồn tại hai luồng tranh luận học thuật lớn:

  1. Tranh luận về độ bền vững kết cấu: Khối phủ 1 lớp đối lập khối phủ 2 lớp. Khối 1 lớp (như Accropode) có độ ổn định ban đầu cao nhưng phá hoại dạng giòn (brittle failure), không có khả năng tự hàn gắn và yêu cầu nghiêm ngặt cấm sóng tràn. Ngược lại, khối 2 lớp (như Tetrapod, RAKUNA IV) có tính tự lấp chỗ (self-healing), cho phép hấp thụ sóng tràn an toàn hơn dù tiêu tốn diện tích bề mặt.
  2. Tranh luận về giá trị hệ số $\gamma_r$: Hằng số cố định đối lập hàm biến thiên. TAW (2002) áp dụng giá trị $\gamma_r$ bất biến, trong khi Bruce và cộng sự (2009) qua 13 loại khối phủ khẳng định $\gamma_r$ tăng tuyến tính với số Iribarren $\xi_{m-1,0}$.

Luận án của Lê Thị Hương Giang đã định vị chính xác tại giao điểm của hai tranh luận trên. Bằng việc so sánh trực tiếp với nghiên cứu của De Jong (1996) trên Tetrapod và Bruce và cộng sự (2009) trên các khối phủ cải tiến, luận án đã chứng minh RAKUNA IV có tính năng tiêu sóng vượt trội nhờ cấu trúc 4 hốc lõm với độ rỗng đạt 56,5% (so với 50,0% của Tetrapod), đồng thời lượng hóa được sự biến thiên của $\gamma_r$ từ 0,38 đến 0,55 theo phổ sóng thực tế.

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

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

Luận án đã mở rộng và hoàn thiện lý thuyết ổn định thủy lực của Van der Meer (1988) và De Jong (1996) cho thế hệ khối phủ dị hình rỗng cao. Bằng việc kiểm chứng quy luật truyền sóng qua đỉnh công trình, nghiên cứu đã chứng minh rằng: "Khi đê cho phép sóng tràn thì ổn định của khối phủ sẽ được gia tăng do một phần năng lượng sóng truyền qua đỉnh đê và mái phía trong, trọng lượng khối phủ và cao trình đỉnh đê giảm đáng kể".

Mô hình lý thuyết tổng quát xác lập quan hệ hàm số giữa chỉ số ổn định khi có sóng tràn $N_{s,ovt}$ và chỉ số ổn định khi không tràn $N_{s,non}$:

$$N_{s,ovt} = N_{s,non} \cdot F_s$$

Trong đó, hệ số gia tăng ổn định $F_s$ được xác lập dưới dạng hàm số thực nghiệm phụ thuộc trực tiếp vào tham số chiều cao lưu không tương đối $R_c/H_{m0}$:

$$F_s = 1 + a \cdot \exp\left( -b \frac{R_c}{H_{m0}} \right)$$

Chuyển dịch mô thức (paradigm shift) quan trọng ở đây là bác bỏ quan điểm thiết kế bảo thủ coi sóng tràn là trạng thái nguy hiểm cần triệt tiêu bằng cách tôn cao đỉnh đê. Luận án cung cấp cơ sở giải tích chứng minh rằng việc cho phép sóng tràn có kiểm soát ($0,63 \le R_c/H_{m0} \le 1,45$) tạo ra hiệu ứng giải tỏa áp lực thủy động ở mái trước, giúp lớp phủ RAKUNA IV chịu được các đỉnh sóng cao hơn mà không bị dịch chuyển phá hoại ($N_{od} < 0,5$).

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

Khung phân tích của luận án tích hợp ba trụ cột cơ học chất lưu và công trình thủy:

  1. Lý thuyết phân tích thứ nguyên $\Pi$-Buckingham: Thiết lập mối quan hệ phi thứ nguyên giữa 6 biến số chi phối $(N_s, H_s, R_c, s_{om}, N_{od}, N_z)$ dẫn xuất đến hàm trạng thái:

$$\Phi\left( N_s, \frac{N_{od}}{N_z^{0,5}}, s_{om}, \frac{R_c}{H_{m0}} \right) = 0$$

  1. Lý thuyết tương tác sóng - môi trường rỗng: Sử dụng phương trình kháng trở Darcy-Forchheimer kết hợp phương trình Navier-Stokes trung bình Reynolds (RANS) mô tả trường áp lực và dòng thấm qua lớp lõi và lớp đệm.
  2. Mô hình toán thể tích chất lưu IH-2VOF: Giải phương trình bảo toàn khối lượng và động lượng hai pha để mô phỏng chính xác quá trình biến dạng sóng, tách sóng vỡ và tái lập trường xoáy bên trong các hốc lõm của khối RAKUNA IV.

Điều kiện biên xác lập của khung phân tích: Mái dốc đồng nhất $\cot \alpha = 1,5$; lớp phủ gồm 2 lớp RAKUNA IV xếp rối ngẫu nhiên với mật độ $\Phi$; chiều dày lớp phủ $2D_n$; chiều rộng thềm đỉnh $G_c = 3D_n$; bề rộng đỉnh đê $B = 3D_n$; đáy máng nằm ngang không xuất hiện sóng vỡ sớm trước chân công trình ($h \ge 2,5 H_{m0}$).

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

Thiết kế nghiên cứu

Nghiên cứu tuân thủ thế giới quan thực chứng (positivism), vận dụng phương pháp hỗn hợp (mixed-methods) kết hợp chặt chẽ giữa mô hình vật lý thu nhỏ và mô hình toán thủy động lực học chuyên sâu.

Thiết kế thí nghiệm tuân thủ nghiêm ngặt tiêu chuẩn tương tự Froude với tỷ lệ hình học $\lambda_l = 1:50$. Nhằm mô phỏng chính xác trường sóng Biển Đông tác động lên công trình biển thực tế ($H_s = 5 \div 10\text{ m}$; $T_p = 5 \div 12\text{ s}$), các đại lượng vật lý trong mô hình được chuyển đổi chuẩn xác:

  • Tỷ lệ thời gian: $\lambda_t = \lambda_l^{1/2} = \sqrt{50} \approx 7,07$ (1 giây trong mô hình ứng với 7,07 giây nguyên hình).
  • Tỷ lệ lưu lượng sóng tràn đơn vị: $\lambda_q = \lambda_l^{3/2} = 50^{1,5} \approx 353,5$ (1 l/s/m trong mô hình tương ứng 353,5 l/s/m nguyên hình).
  • Tỷ lệ lực tác động: $\lambda_P = \lambda_l^3 = 50^3 = 125.000$.

Đặc biệt, để khắc phục hiệu ứng kích thước (scale effect) của dòng thấm trong mô hình Froude, luận án đã áp dụng phương pháp của Burcharth và cộng sự (1999): lớp đệm và lớp lõi đê được thiết kế theo một tỷ lệ riêng biệt $\lambda_{l,core} = 9$ với đường kính danh nghĩa $D_{n50} = 2,3\text{ cm}$ để bảo đảm tính tương tự của hệ số thấm Reynolds, tránh hiện tượng áp lực nước dâng giả tạo bên trong thân đê làm sai lệch độ ổn định của khối phủ mặt ngoài.

Quy trình nghiên cứu rigorous

Hệ thống thiết bị và quy trình thực nghiệm được chuẩn hóa tại Phòng thí nghiệm Thủy lực Tổng hợp – Trường Đại học Thủy lợi:

  • Máng sóng Hà Lan: Kích thước dài 45 m, cao 1,2 m, rộng 1,0 m; trang bị máy tạo sóng điện - thủy lực tích hợp hệ thống triệt tiêu sóng phản xạ chủ động (Active Reflection Compensation - ARC).
  • Thiết bị đo đạc: Hệ thống 3 đầu đo sóng kiểu điện dung đặt trước công trình tại các cự ly $0,81\text{ m}$ và $0,22\text{ m}$. Thuật toán phân tách sóng tới và sóng phản xạ của Zelt và Skjelbreia (1992) được lập trình xử lý tín hiệu trực tiếp.
  • Hệ thống ghi nhận phá hoại: Khối phủ RAKUNA IV chế tạo bằng vữa xi măng cường độ cao ($D_n = 4,11\text{ cm}$) được sơn thành 4 dải màu phân biệt (Trắng – Xanh – Vàng – Đỏ) theo cao trình để nhận diện tọa độ không gian. Hai camera độ phân giải cao ghi hình đồng bộ liên tục từ mặt chính diện và mặt cắt kính thành máng.
  • Quy trình gia tải sóng: Mỗi kịch bản ổn định vận hành kéo dài đúng $N_z = 3000$ con sóng ngẫu nhiên JONSWAP hoặc dừng khi diện tích lõi bị hở lộ vượt quá $2D_n^2$. Lưu lượng sóng tràn được thu gom bằng máng chuyên dụng và đo thể tích chính xác qua $N_z = 1000$ con sóng đầu tiên.
Sơ đồ phân bố thiết bị trong máng sóng 45m:
                                                                                   [Thùng đo thể tích Vovt]

Data và phân tích

Ma trận thực nghiệm gồm 58 bài toán thí nghiệm (trong đó có 36 thí nghiệm chính và các thí nghiệm lặp cho RAKUNA IV, cùng 8 thí nghiệm đối chứng cho Tetrapod). Các thông số đầu vào kiểm soát bao gồm:

  • Chiều cao sóng mô hình: $H_{m0} = 0,145 \div 0,214\text{ m}$.
  • Chu kỳ đỉnh: $T_p = 1,50 \div 2,60\text{ s}$.
  • Độ dốc sóng: $s_{om} = 0,015 \div 0,050$ (chia thành 9 dải phổ).
  • Mực nước trong máng: 4 mức độ sâu tĩnh ($52,5\text{ cm}$; $55,0\text{ cm}$; $57,5\text{ cm}$; $60,0\text{ cm}$).
  • Chiều cao lưu không: $R_c = 0,136 \div 0,211\text{ m}$, tương ứng dải tỷ số $R_c/H_{m0} = 0,63 \div 1,45$.

Mức độ hư hỏng được xác định thông qua số khối dịch chuyển tương đối $N_{od}$ trên dải bề rộng $B = D_n$:

$$N_{od} = \frac{N_{omov}}{B / D_n}$$

Trong đó $N_{omov}$ tính tổng toàn bộ các khối bị bật khỏi mái, trượt quá cự ly $0,5D_n$, xoay lộn góc lớn hoặc xuất hiện hiện tượng rung lắc/đá nhảy (rocking).

Mô hình số IH-2VOF được hiệu chỉnh với các hệ số cản dòng thấm: hệ số tuyến tính $\alpha = 200 \div 800$ và hệ số cản phi tuyến $\beta = 0,6 \div 1,2$, cho phép tái hiện đường quá trình mặt thoáng sóng và phổ năng lượng với sai số tương quan $R^2 > 0,92$.

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

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

  1. Quy luật gia tăng ổn định thủy lực dưới tác động của sóng tràn: Dữ liệu thực nghiệm chứng minh khi tỷ số lưu không tương đối $R_c/H_{m0}$ giảm từ 1,45 xuống 0,63, chỉ số ổn định $N_s$ của khối RAKUNA IV tăng từ 2,35 lên vượt mức 3,10 mà kết cấu vẫn nằm trong ngưỡng hư hỏng cho phép ($N_{od} \le 0,5$). Biểu thức hệ số gia tăng ổn định thực nghiệm $F_s$ được xác lập:

    $$F_s = 1 + 0,22 \exp\left( -0,85 \frac{R_c}{H_{m0}} \right)$$

    Kết quả này khẳng định mái đê phía biển trở nên ổn định hơn từ 10% đến 22% so với điều kiện không cho phép sóng tràn.

  2. Đặc tính biến thiên phi tuyến của hệ số chiết giảm sóng tràn $\gamma_r$: Khối phủ RAKUNA IV thể hiện khả năng tiêu giảm sóng vượt trội nhờ các hốc lõm tạo xoáy cục bộ. Hệ số chiết giảm sóng tràn thực đo dao động trong khoảng $\gamma_r = 0,38 \div 0,55$, thấp hơn đáng kể so với đá hộc ($\gamma_r = 0,55$) và khối Accropode ($\gamma_r = 0,46$). Tương quan giải tích giữa $\gamma_r$ và số Iribarren được xác lập:

    $$\gamma_r = 0,085 \cdot \xi_{m-1,0} + 0,15 \quad (\text{với } 2,8 \le \xi_{m-1,0} \le 5,2)$$

    Điều này chứng minh rằng sóng có chu kỳ càng dài (hoặc mái đê càng dốc) thì hiệu quả cảm nhận độ nhám càng suy giảm.

Hệ số chiết giảm sóng tràn γr theo số Iribarren ξm-1,0:
γr
            2.5     3.0     3.5     4.0     4.5     5.0

  1. Lý giải bản chất vật lý của hiệu ứng "Đệm nước" (Water Cushion Effect): Kết quả mô phỏng trường véc-tơ lưu tốc từ mô hình IH-2VOF đã làm sáng tỏ hiện tượng đệm nước. Tác giả luận án chỉ rõ: "Do hiệu ứng 'đệm nước' này mà con sóng tiếp theo leo lên mái đê không cảm nhận được nhiều độ nhám của các khối phủ nữa và gần như tác động lên một mái dốc nhẵn dẫn đến sóng leo, sóng tràn lớn hơn". Khi sóng rút trên mái dốc $1:1,5$, độ rỗng lớn $56,5%$ và cấu tạo 4 hốc lõm giữ lại một lượng nước chậm thoát, tạo thành lớp màng thủy động bôi trơn làm suy giảm ma sát cục bộ đối với đợt sóng tràn kế tiếp.

  2. Ưu thế cài khóa ngăn chặn hiện tượng phá hoại giòn: Khác với khối Tetrapod dễ xuất hiện hiện tượng xoay trượt dây chuyền khi chân đê bị xói, khối RAKUNA IV nhờ các cạnh vát và chân cài chặt vào hốc lõm lân cận tạo nên liên kết đa điểm vững chắc. Hiện tượng "đá nhảy" (rocking) giảm hơn 40% so với Tetrapod ở cùng cấp tải trọng $H_{m0} = 0,20\text{ m}$.

Implications đa chiều

  • Về mặt lý thuyết: Bổ sung vào kho tàng cơ học công trình thủy hệ thống phương trình trạng thái mới mô tả ổn định của khối phủ rỗng cao trong điều kiện sóng tràn, mở rộng biên độ ứng dụng của các tiêu chuẩn quốc tế như EurOtop và TAW.
  • Về mặt phương pháp luận: Chuẩn hóa quy trình kết hợp mô hình vật lý máng sóng ARC với mô hình số IH-2VOF, xác lập phương pháp đo đạc phân tách sóng ngẫu nhiên chuẩn xác cho các phòng thí nghiệm thủy lực ven biển.
  • Về mặt thực tiễn kỹ thuật: Cung cấp thông số thiết kế trực tiếp cho các kỹ sư công trình thủy. Tại dự án đê chắn sóng Cảng Nghi Sơn (Thanh Hóa), việc áp dụng khối phủ RAKUNA IV có tính đến hệ số gia tăng ổn định $F_s$ và cho phép sóng tràn hợp lý đã giúp hạ thấp cao trình đỉnh đê từ $+9,0\text{ m}$ xuống $+6,5\text{ m}$, giảm trọng lượng cấu kiện từ 15 tấn (Tetrapod) xuống 10 tấn (RAKUNA IV), tiết kiệm 18,5% khối lượng bê tông và giảm hơn 15% tổng chi phí xây lắp trên 100 m chiều dài đê.
  • Về mặt chính sách: Cung cấp cơ sở khoa học vững chắc để các cơ quan quản lý (Bộ Nông nghiệp & PTNT, Bộ Giao thông Vận tải) cập nhật các tiêu chuẩn thiết kế đê biển Việt Nam theo hướng thích ứng biến đổi khí hậu và tối ưu hóa chi phí đầu tư công.

Limitations và Future Research

Mặc dù đạt được những kết quả đột phá, luận án thẳng thắn thừa nhận các giới hạn học thuật:

  • Giới hạn không gian 2D: Thí nghiệm thực hiện trong máng sóng phẳng 2 chiều với góc sóng tới vuông góc ($\beta = 0^\circ$), chưa phản ánh hết tác động của sóng xiên, sóng nhiễu xạ và dòng chảy ven bờ phức tạp trong không gian 3 chiều.
  • Giới hạn mái dốc và cấu trúc xếp khối: Nghiên cứu chỉ tập trung vào một độ dốc mái đê duy nhất là $\cot \alpha = 1,5$ và phương pháp xếp rối 2 lớp, chưa mở rộng cho các mái dốc thoải hơn ($1:2; 1:3$) hoặc kiểu xếp định hướng 1 lớp.
  • Giới hạn về tải trọng gãy vỡ vật liệu: Thí nghiệm mô hình thủy lực Froude chưa mô phỏng được ứng suất uốn gãy bên trong thân cấu kiện bê tông khi xảy ra va đập động năng lớn ở tỷ lệ nguyên hình.

Chương trình nghiên cứu tương lai cần tập trung vào các hướng mở rộng:

  1. Thực hiện thí nghiệm mô hình 3D trong bể sóng rộng (wave basin) để đánh giá ảnh hưởng của góc sóng tới $\beta = 15^\circ \div 60^\circ$ và phổ sóng đa hướng.
  2. Ứng dụng mô hình số liên hợp FEM-DEM (Finite/Discrete Element Method) để phân tích ứng suất tập trung tại 4 hốc lõm nhằm ngăn ngừa rủi ro nứt gãy chân cấu kiện.
  3. Nghiên cứu tích hợp yếu tố nước biển dâng và triều cường cực đoan do biến đổi khí hậu tác động đến năng lực tiêu giảm sóng tràn của đê trong tầm nhìn 50 - 100 năm.
  4. Mở rộng thử nghiệm cấu kiện RAKUNA IV sử dụng bê tông tính năng siêu cao (UHPC) hoặc bê tông địa polyme thân thiện môi trường biển.

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

Công trình nghiên cứu tạo ra những tác động đa tầng sâu rộng:

  • Tác động học thuật: Mở ra hướng tiếp cận mới trong việc nghiên cứu thủy lực công trình bảo vệ bờ; cung cấp dữ liệu thực nghiệm chuẩn xác phục vụ hiệu chỉnh các mô hình toán quốc tế; dự kiến thu hút lượng trích dẫn cao trong các tạp chí chuyên ngành như Coastal Engineering, Ocean Engineering và các hội nghị quốc tế COPEDEC, ICCE.
  • Tác động ngành và công nghiệp xây dựng thủy: Chuyển giao công nghệ khối phủ dị hình hiện đại từ Nhật Bản vào thực tiễn Việt Nam, phá vỡ thế độc quyền của các dạng khối truyền thống, thúc đẩy ngành công nghiệp đúc cấu kiện bê tông đúc sẵn đạt tiêu chuẩn quốc tế.
  • Tác động chính sách và xã hội: Tối ưu hóa vốn đầu tư xây dựng cơ bản cho hệ thống cảng biển quốc gia; nâng cao năng lực chống chịu thiên tai của các khu kinh tế ven biển trọng điểm (Nghi Sơn, Dung Quất, Lạch Huyện); bảo vệ an toàn dân sinh trước các biến cố bão biển cực đoan.

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

  • Nghiên cứu sinh và học viên sau đại học: Tiếp cận phương pháp luận phân tích thứ nguyên $\Pi$-Buckingham chuẩn mực, quy trình vận hành máng sóng ARC và kỹ thuật mô phỏng IH-2VOF.
  • Các nhà khoa học và giảng viên chuyên ngành Công trình Thủy/Biển: Sở hữu nguồn số liệu thực nghiệm và hệ thống công thức giải tích chuẩn xác để phục vụ giảng dạy và phát triển các đề tài nghiên cứu chuyên sâu.
  • Các kỹ sư thiết kế và đơn vị tư vấn xây dựng công trình biển (TEDI, VNC, Portcoast): Ứng dụng trực tiếp hệ số $F_s$ và $\gamma_r$ của khối RAKUNA IV để tối ưu hóa kích thước mặt cắt đê, giảm giá thành xây dựng các dự án cảng nước sâu.
  • Các nhà hoạch định chính sách và chủ đầu tư hạ tầng cảng biển: Có cơ sở dữ liệu kinh tế - kỹ thuật tin cậy để thẩm định phương án đầu tư, lựa chọn giải pháp kết cấu đê chắn sóng bền vững và tiết kiệm ngân sách nhà nước.

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? Trả lời: Đóng góp độc đáo nhất là việc xây dựng thành công công thức thực nghiệm xác định hệ số gia tăng ổn định $F_s = f(R_c/H_{m0})$ cho khối phủ RAKUNA IV xếp 2 lớp trên mái nghiêng $1:1,5$. Luận án đã mở rộng lý thuyết ổn định của Van der Meer (1988) và De Jong (1996) từ phạm vi khối đá/Tetrapod truyền thống sang thế hệ khối phủ dị hình có hốc lõm với độ rỗng cao $n = 56,5%$ trong điều kiện có sóng tràn.

2. Điểm mới trong phương pháp nghiên cứu so với các công bố quốc tế trước đây là gì? Trả lời: So với nghiên cứu của Mase và cộng sự (2011) (chỉ làm việc trên đê hỗn hợp ngang không tràn) và Bruce và cộng sự (2009) (chỉ đo sóng tràn mà không đánh giá phá hoại ổn định), luận án đã tích hợp đồng thời hai bài toán tương hỗ: đánh giá chuyển động phá hoại $N_{od}$ qua $N_z = 3000$ con sóng và đo lưu lượng sóng tràn $q$ qua $N_z = 1000$ con sóng trên cùng một hệ mô hình vật lý Froude $1:50$, kết hợp chia dải 4 màu cấu kiện và mô phỏng kiểm chứng chi tiết trường dòng chảy bằng mô hình toán IH-2VOF.

3. Phát hiện bất ngờ nhất về mặt thủy lực học trong luận án là gì và dữ liệu nào chứng minh? Trả lời: Phát hiện bất ngờ nhất là khi số Iribarren $\xi_{m-1,0}$ tăng cao (sóng dài, chu kỳ lớn), hệ số chiết giảm sóng tràn $\gamma_r$ của RAKUNA IV tăng mạnh từ 0,38 lên 0,55 (tức khả năng giảm sóng tràn bị suy giảm). Dữ liệu mô phỏng IH-2VOF chứng minh hiện tượng này bắt nguồn từ hiệu ứng "đệm nước": lượng nước đọng trong 4 hốc lõm không kịp thoát ra đã tạo thành một màng chất lưu bôi trơn làm cho con sóng kế tiếp trượt qua như trên bề mặt mái nhẵn.

4. Quy trình thực nghiệm trong luận án có bảo đảm khả năng tái lập (replication protocol) không? Trả lời: Hoàn toàn bảo đảm. Luận án cung cấp chi tiết mọi thông số kỹ thuật: kích thước hình học mô hình khối RAKUNA IV ($D_n = 4,11\text{ cm}$), kích thước máng sóng ($45 \times 1,2 \times 1,0\text{ m}$), vị trí 3 đầu đo sóng Zelt & Skjelbreia ($0,81\text{ m}$ và $0,22\text{ m}$), tỷ lệ lõi riêng biệt theo Burcharth $\lambda_{l,core} = 9$, phổ sóng đầu vào JONSWAP với 9 độ dốc sóng, và 58 kịch bản thí nghiệm lặp được chuẩn hóa quy trình từng bước.

5. Luận án đã vạch ra chương trình nghiên cứu 10 năm tiếp theo như thế nào? Trả lời: Luận án định hướng lộ trình nghiên cứu tập trung vào: (1) Thử nghiệm mô hình 3D trong bể sóng lớn chịu tác động của sóng xiên đa hướng; (2) Mô phỏng số cơ học phá hủy liên hợp FEM-DEM kiểm tra độ bền mỏi và nứt gãy của 4 chân khối phủ; (3) Xây dựng cẩm nang kỹ thuật và phần mềm tự động hóa thiết kế đê chắn sóng mái nghiêng cho phép sóng tràn ứng dụng khối phủ RAKUNA IV tại bờ biển Việt Nam.

Kết luận

Luận án tiến sĩ của tác giả Lê Thị Hương Giang đã hoàn thành xuất sắc các mục tiêu nghiên cứu với 5 đóng góp cốt lõi:

  1. Thiết lập công thức giải tích thực nghiệm tính toán định lượng mức độ gia tăng ổn định $F_s$ của khối phủ RAKUNA IV trên đê chắn sóng mái nghiêng khi có sóng tràn.
  2. Xác định chính xác hệ số chiết giảm sóng tràn $\gamma_r$ và quy luật phụ thuộc tuyến tính của nó vào số Iribarren $\xi_{m-1,0}$ đối với cấu kiện có 4 hốc lõm.
  3. Lý giải thành công bản chất vật lý của hiệu ứng "đệm nước" và cơ chế tiêu tán năng lượng thủy động thông qua mô hình toán chất lưu hai pha IH-2VOF.
  4. Chuẩn hóa quy trình mô hình vật lý máng sóng kết hợp xử lý hiệu ứng kích thước dòng thấm lõi đê theo tiêu chuẩn quốc tế tiên tiến.
  5. Chứng minh hiệu quả kinh tế - kỹ thuật vượt trội thông qua bài toán ứng dụng thực tế tại đê chắn sóng Cảng Nghi Sơn, tiết kiệm hơn 15% chi phí đầu tư.

Công trình đã tạo ra bước tiến quan trọng trong việc chuyển dịch mô thức thiết kế công trình thủy ven biển tại Việt Nam, mở ra ba nhánh nghiên cứu học thuật mới về thủy động lực học công trình rỗng cao, đồng thời để lại giá trị ứng dụng thực tiễn lâu dài cho sự nghiệp phát triển bền vững hạ tầng kinh tế biển quốc gia.