Tổng quan về luận án

Hệ thống đê biển Bắc Bộ trải dài qua các địa bàn trọng điểm như Hải Phòng, Nam Định, Thanh Hóa, Hà Tĩnh đóng vai trò hộ đê xung yếu, bảo vệ các đô thị ven biển, khu công nghiệp và đời sống của hàng triệu cư dân. Tuy nhiên, trước tác động khốc liệt của biến đổi khí hậu, nước biển dâng và triều cường kết hợp bão mạnh, hệ thống đê biển khu vực này đang bộc lộ những điểm yếu chí mạng về mặt hình thái học công trình: cao trình đỉnh đê hiện hữu tương đối thấp (chỉ dao động từ +4.0 m đến +5.5 m), bề rộng mặt đê hẹp (3.0 – 5.0 m), mái đê dốc (từ 1:3 đến 1:4) và hầu hết tiếp xúc trực diện với biển. Trong bối cảnh quỹ đất hạn hẹp không thể mở rộng hành lang chân đê và nguồn lực tài chính công chưa thể đáp ứng việc đắp tôn cao toàn bộ thân đê, giải pháp xây dựng tường đỉnh thấp ($W/H_s \le 0.5$) kết hợp thềm trước trở thành một lựa chọn kỹ thuật - kinh tế mang tính tất yếu và cấp bách.

flowchart TD
    A["Tác động thiên tai: Bão, Triều cường, Nước dâng"] --> B["Sóng leo và Sóng tràn qua đỉnh đê"]
    B --> C{"Cơ cấu hình học đỉnh đê"}
    C -->|"Không có tường đỉnh"| D["Dòng chảy tràn vận tốc lớn (3-5 m/s)"]
    C -->|"Có tường đỉnh đứng sát mép (S = 0)"| E["Sóng va đập tạo sóng bắn Hb dội thẳng xuống mặt đê"]
    C -->|"Có tường đỉnh thấp + Thềm trước (S > 0)"| F["Tiêu tán năng lượng & Chiết giảm lưu lượng sóng tràn q"]
    D --> G["Xói lở mái trong & Thấm bão hòa thân đê"]
    E --> G
    F --> H["Bảo vệ ổn định thân đê & An toàn công trình"]
    G --> I["Vỡ đê thảm khốc"]

Khoảng trống nghiên cứu (research gap) cốt lõi xuất phát từ sự thiếu hụt cơ sở lý thuyết và thực nghiệm chuẩn xác về tương tác thủy động lực học giữa sóng, tường đỉnh thấp và thềm trước. Các tiêu chuẩn quốc tế kinh điển như sổ tay TAW (2002) của Hà Lan hay hướng dẫn EurOtop (2007) chủ yếu xây dựng công thức dựa trên mô hình đê mái dốc truyền thống hoặc đê chắn sóng tường đứng trọng lực cỡ lớn. Khi áp dụng cho đê có tường đỉnh thấp, TAW (2002) chuyển đổi tường đứng thành một mái dốc quy đổi tương đương ($1:1.0$), từ đó làm sai lệch chỉ số Iribarren ($\xi$) và tạo ra hiện tượng "sóng vỡ giả" trong tính toán, đồng thời mặc định triệt tiêu ảnh hưởng của tường đỉnh trong điều kiện sóng không vỡ ($\gamma_v = 1.0$).

Luận án tiến sĩ kỹ thuật của Nguyễn Văn Thìn (2014) với đề tài "Nghiên cứu sóng tràn qua đê biển có tường đỉnh ở Bắc Bộ" (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) đã giải quyết triệt để khoảng trống này thông qua hệ thống câu hỏi nghiên cứu và giả thuyết khoa học tường minh:

  • Câu hỏi nghiên cứu 1 (RQ1): Cơ chế thủy lực chi phối lưu lượng sóng tràn trung bình ($q$) qua đê biển khi có sự xuất hiện của tường đỉnh thấp và thềm trước diễn ra như thế nào dưới tác động của sóng đều và sóng ngẫu nhiên?
  • Câu hỏi nghiên cứu 2 (RQ2): Mối quan hệ định lượng giữa các thông số hình học (chiều cao tường $W$, bề rộng thềm trước $S$) và thông số sóng ($H, T$) đối với chiều cao sóng bắn ($H_b$) và chiều dày lớp nước tràn ($H_t$) là gì?
  • Câu hỏi nghiên cứu 3 (RQ3): Làm thế nào để mô hình hóa chính xác tương tác sóng – tường bằng công cụ toán học thủy động lực hai chiều thay thế cho các giả thiết hình học quy đổi thiếu chuẩn xác của TAW (2002)?
  • Giả thuyết khoa học 1 (H1): Sự hiện diện của thềm trước ($S > 0$) tạo ra một vùng đệm thủy lực làm suy giảm năng lượng động học của lưỡi sóng trước khi va vào tường, giúp tăng hiệu quả chiết giảm sóng tràn của tường đỉnh thấp lên 35 – 65% so với tường không có thềm.
  • Giả thuyết khoa học 2 (H2): Hệ số ảnh hưởng tổng hợp $\gamma_v$ không phải là một hằng số hình học tĩnh mà là hàm số phụ thuộc trực tiếp vào tỷ số không thứ nguyên $(W/R_c)$, $(S/H)$ và chỉ số vỡ sóng $\xi$, có tác dụng chiết giảm sóng tràn trong cả hai chế độ sóng vỡ và sóng không vỡ.
  • Giả thuyết khoa học 3 (H3): Chiều cao sóng bắn $H_b$ tỷ lệ nghịch với bề rộng tương đối của thềm trước $(S/H)$ và tuân theo một quy luật phi tuyến phụ thuộc vào số Iribarren.

Quy mô thực nghiệm của luận án bao gồm chuỗi 40 kịch bản thí nghiệm sóng đều độc lập trên 10 cấu hình mô hình vật lý tỷ lệ 1:10 tại máng sóng DELFT (Hà Lan) dài 45 m, kết hợp kế thừa cơ sở dữ liệu 225 thí nghiệm sóng ngẫu nhiên của Thiều Quang Tuấn (2013). Kết quả mang lại đóng góp đột phá khi thiết lập được công thức giải tích xác định hệ số chiết giảm tổng hợp $\gamma_v$ cho sóng đều, biểu đồ quan hệ sóng bắn $H_b$, và phần mềm chuyên dụng tính toán sóng tràn ứng dụng trực tiếp cho tuyến đê biển Giao Thủy (Nam Định).

Literature Review và Positioning

Nghiên cứu về sóng leo và sóng tràn qua công trình bảo vệ bờ biển đã trải qua hơn sáu thập kỷ phát triển với các mốc tiến hóa lý thuyết quan trọng:

timeline
    title Tiến trình Phát triển Lý thuyết Tính toán Sóng tràn qua Đê biển
    1955 : Saville : Thí nghiệm sóng đơn đầu tiên xác định sóng leo
    1980 : Owen : Mô hình thực nghiệm sóng tràn cho đê mái nhẵn và độ nhám
    1992 : De Waal & Van der Meer : Khái niệm độ thiếu hụt lưu không (Ru2% - Rc)/Hs
    1995 : Van der Meer & Janssen : Phân tách chế độ sóng vỡ và không vỡ qua Rc/Hs
    2002 : TAW (Hà Lan) : Sổ tay kỹ thuật chuẩn hóa hệ số giảm sóng tràn quy đổi
    2007 : EurOtop : Hợp nhất tiêu chuẩn Châu Âu cho công trình bảo vệ bờ
    2009 : Thiều Quang Tuấn et al. : Phát hiện hiện tượng sóng vỡ giả trong TAW (2002) với S=0
    2014 : Nguyễn Văn Thìn : Hoàn thiện lý thuyết tường đỉnh thấp có thềm trước S > 0
  1. Giai đoạn nền tảng thực nghiệm cổ điển (1955 – 1995): Saville (1955) đặt nền móng với các thí nghiệm sóng đơn trên mái nghiêng. Owen (1980) thực hiện hàng ngàn thử nghiệm mô hình sóng ngẫu nhiên để thiết lập công thức hàm mũ kinh điển tính toán lưu lượng sóng tràn trung bình $q$ cho đê mái dốc nhẵn không thấm. De Waal & Van der Meer (1992) đưa thêm tham số độ thiếu hụt lưu không đỉnh đê $(R_{u2%} - R_c)/H_s$. Đến năm 1995, Van der Meer & Janssen tạo bước ngoặt khi liên kết trực tiếp lưu lượng tràn không thứ nguyên với độ lưu không tương đối $R_c/H_s$, phân định ranh giới thủy lực giữa sóng vỡ ($\xi_{0m} \le 2.0$) và sóng không vỡ ($\xi_{0m} > 2.0$).
  2. Giai đoạn tiêu chuẩn hóa quốc tế (2002 – 2007): Sổ tay TAW (2002) và cẩm nang EurOtop (2007) tích hợp các hệ số chiết giảm ảnh hưởng của độ nhám mái đê ($\gamma_f$), cơ đê ($\gamma_b$), góc truyền sóng ($\gamma_\beta$) và tường đỉnh ($\gamma_v$). Tuy nhiên, TAW (2002) quy định tường đỉnh chỉ phát huy tác dụng khi sóng vỡ thông qua việc quy đổi tường thành mái dốc $1:1.0$. Khi tường thẳng đứng ($w = 90^\circ$), TAW cố định $\gamma_v = 0.65$, còn khi sóng không vỡ thì mặc định $\gamma_v = 1.0$.
  3. Giai đoạn phê phán và hiệu chỉnh tại Việt Nam (2009 – 2013): Thiều Quang Tuấn và cộng sự (2009) đã chứng minh rằng phương pháp mái dốc quy đổi của TAW (2002) làm tăng vọt chỉ số Iribarren, dẫn đến hiện tượng sóng vỡ giả và đánh giá thấp đáng kể năng lực tiêu giảm sóng của tường đỉnh. Nhóm tác giả đề xuất liên hệ phụ thuộc vào tỷ số $W/(W+R_c)$ nhưng chỉ giới hạn cho trường hợp tường đặt sát mép đỉnh ($S = 0$). Đến năm 2013, Thiều Quang Tuấn mở rộng cho sóng ngẫu nhiên có xét đến thềm trước ($S > 0$) với 225 thí nghiệm, tách biệt $\gamma_v = \gamma_w \cdot \gamma_s$.
Tiêu chí So sánh Sổ tay TAW (2002) / EurOtop (2007) Nghiên cứu Đê biển GWK (Đức) & Isahaya (Nhật Bản) Luận án Nguyễn Văn Thìn (2014)
Mô hình tương tác Quy đổi hình học tương đương sang mái dốc $1:1.0$ Tường đỉnh khối lớn, đê kết hợp đường cao tốc Tách biệt động lực học tường ($W$) và thềm trước ($S$)
Chế độ sóng không vỡ Mặc định $\gamma_v = 1.0$ (Không có tác dụng giảm tràn) Tường khối phản xạ toàn phần sóng không vỡ $\gamma_v < 1.0$ (Chứng minh tường vẫn giảm tràn rõ rệt)
Tham số thềm trước ($S$) Bỏ qua hoặc xem như cơ đê thông thường Thềm rộng kết cấu cứng chịu tải sóng lớn Chuẩn hóa tỷ số không thứ nguyên $S/H$ trong công thức giải tích
Hiện tượng sóng bắn ($H_b$) Hoàn toàn không lượng hóa Quan sát định tính qua cảm biến áp lực Xây dựng đường cong quan hệ $(H_b/H) - (S/H)\xi$
Công cụ mô phỏng số Mô hình 1D sóng nông truyền thống Mô hình phần tử hữu hạn cấu trúc - chất lưu Tích hợp 2DV NLSW và RANS-VOF (COBRAS-UC)

Trích dẫn nguyên văn từ văn bản nguồn khẳng định tính kế thừa và định vị học thuật của luận án: "Nghiên cứu của Thiều Quang Tuấn (2013) đã đánh giá được ảnh hưởng của tường đỉnh và thềm trước tường đến lưu lượng sóng tràn qua đê biển. Tuy nhiên, Thiều Quang Tuấn (2013) chưa xem xét đến tính chất tương tác sóng – tường và sự thay đổi của dòng chảy sóng tràn qua đê khi có sự hiện diện của tường đỉnh. Ngoài ra, ảnh hưởng chiết giảm của tường đối với sóng đều cũng chưa được đề cập trong nghiên cứu." Luận án định vị chính xác vào việc bổ khuyết các quy luật động lực học sóng đều và cơ chế sóng bắn - dòng chảy tràn chi tiết.

Đó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 giả định cơ bản của TAW (2002) và EurOtop (2007) về việc triệt tiêu ảnh hưởng của tường đỉnh trong điều kiện sóng không vỡ. Bằng chứng thực nghiệm khẳng định tường đỉnh thấp và thềm trước vẫn tạo ra sự cản trở thủy lực mạnh mẽ, làm lệch hướng vector vận tốc dòng chảy hướng lên trên và làm tiêu tán một phần động năng sóng leo.

graph LR
    subgraph "Khung Lý thuyết Tương tác Thủy động lực học"
        A["Năng lượng Sóng tới (H, T)"] --> B["Sóng leo trên Mái dốc tan alpha = 1/3"]
        B --> C["Vùng Thềm trước (S/H)"]
        C -->|"Tiêu tán ma sát đáy & Dòng phản xạ"| D["Đệm Thủy lực cản sóng"]
        D --> E["Va đập Tường đứng (W/Rc)"]
        E --> F["Sóng bắn Hb văng vào không khí"]
        E --> G["Dòng chảy tràn Ht qua đỉnh tường"]
        F --> H["Rơi tự do & Tiêu hao động năng"]
        G --> I["Lưu lượng tràn q giảm 35-65%"]
    end

Mô hình lý thuyết được xây dựng dựa trên việc phân rã hệ số chiết giảm tổng hợp $\gamma_v$ thành tích số của hệ số chiết giảm do chiều cao tường ($\gamma_w$) và hệ số chiết giảm do bề rộng thềm trước ($\gamma_s$):

$$\frac{1}{\gamma_v} = \frac{1}{\gamma_w \cdot \gamma_s} = \left(1 + c_1 \frac{W}{R_c} \frac{1}{\xi}\right) \left(1 + c_2 \frac{S}{H} \frac{1}{\xi}\right)$$

Qua phân tích hồi quy phi tuyến tối ưu hóa đường phân giác đối chiếu giữa dữ liệu thực đo và tính toán, luận án đã xác lập chính xác bộ hệ số thực nghiệm cho trường hợp sóng đều: $c_1 = 1.26$ và $c_2 = 1.44$. Phương trình chuẩn hóa trở thành:

$$\frac{1}{\gamma_v} = \left(1 + 1.26 \frac{W}{R_c} \frac{1}{\xi}\right) \left(1 + 1.44 \frac{S}{H} \frac{1}{\xi}\right)$$

Công thức này chứng minh rằng hiệu quả chiết giảm lưu lượng sóng tràn tỷ lệ thuận với chiều cao tường tương đối $W/R_c$ và bề rộng thềm tương đối $S/H$, đồng thời tỷ lệ nghịch với chỉ số Iribarren $\xi$.

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 lý thuyết thủy động lực học:

  1. Lý thuyết sóng phi tuyến nước nông (Non-Linear Shallow Water - NLSW): Mô tả quá trình biến đổi sóng từ vùng nước sâu vào bãi nông trước đê ($i = 1/100$) và diễn biến sóng leo trên mái dốc.
  2. Lý thuyết động lực học va chạm chất lưu và dòng tia tự do: Lượng hóa hiện tượng sóng bắn. Trích dẫn định nghĩa từ văn bản gốc: “Sóng bắn là một phần sóng tràn nằm trong không khí không bám sát bề mặt đê”. Chiều cao sóng bắn $H_b$ được xác định qua hàm tương quan tường minh:

$$\frac{H_b}{H} = f\left(\frac{S}{H} \cdot \xi\right)$$

  1. Hệ phương trình Navier-Stokes trung bình Reynolds kết hợp phương pháp thể tích chất lưu (2DV RANS-VOF): Sử dụng hệ mã nguồn mở COBRAS-UC để phân giải cấu trúc trường áp lực sóng xung kích và phân bố vận tốc tức thời tại góc giao giữa thềm đê và mặt tường thẳng đứng.

Điều kiện biên áp dụng được xác định nghiêm ngặt: độ dốc mái đê $\cot \alpha = 3.0 \div 4.0$, chiều cao tường tương đối $W/H_s \le 0.5$ (theo đúng định nghĩa trích dẫn: “Tường đỉnh thấp là tường đỉnh trên đê có $(W/H_s \le 0.5)$”), tỷ số thềm trước $S/H \le 2.0$, và độ lưu không tương đối $R_c/H \ge 0.6$.

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

Thiết kế nghiên cứu

Nghiên cứu được xây dựng trên nền tảng triết học thực chứng (positivism) kết hợp suy diễn cơ học thủy khí (hydrodynamic mechanics). Thiết kế đa phương pháp (triangulation) đan kết ba nhánh tiếp cận: mô hình vật lý thu nhỏ, mô phỏng số học 2DV và kiểm chứng thực địa.

flowchart LR
    subgraph "Thiết kế Tam giác hóa Phương pháp (Triangulation)"
        A["Mô hình Vật lý Tỷ lệ 1:10<br>(Máng sóng DELFT dài 45m)"] <-->|"Hiệu chỉnh & Kiểm chứng biên"| B["Mô hình Toán Thủy động lực<br>(NLSW & 2DV RANS-VOF)"]
        A -->|"Trích xuất dữ liệu thực nghiệm"| C["Công thức Giải tích gamma_v & Hb"]
        B -->|"Mô phỏng trường áp lực tức thời"| C
        C -->|"Tích hợp công cụ phần mềm"| D["Ứng dụng Công trình Thực tế<br>(Đê biển Giao Thủy, Nam Định)"]
    end

Quy luật tỷ lệ mô hình tuân thủ nghiêm ngặt tiêu chuẩn tương tự Froude với tỷ xích hình học $\lambda_l = 10$:

  • Tỷ xích thời gian: $\lambda_t = \sqrt{\lambda_l} = \sqrt{10} \approx 3.162$
  • Tỷ xích vận tốc: $\lambda_u = \sqrt{\lambda_l} = \sqrt{10} \approx 3.162$
  • Tỷ xích lưu lượng đơn vị: $\lambda_q = \lambda_l^{3/2} = 10^{1.5} \approx 31.623$
  • Tỷ xích lực tác dụng: $\lambda_f = \lambda_l^3 = 10^3 = 1000$

Để loại trừ sai số do hiệu ứng nhớt Reynolds gây ra ở tỷ xích nhỏ đối với đê mái nhám, mô hình vật lý sử dụng bề mặt mái nhẵn, không thấm nước, đảm bảo bảo toàn động năng sóng leo và tính chuẩn xác của các tham số đo đạc.

Quy trình nghiên cứu rigorous

Hệ thống thiết bị tại Phòng thí nghiệm Thủy lực - Trường Đại học Thủy lợi:

  • Máng sóng DELFT (Hà Lan): Chiều dài tổng thể 45 m, chiều dài hữu dụng 42 m, chiều cao 1.2 m, chiều rộng 1.0 m.
  • Hệ thống tạo sóng và triệt tiêu phản xạ: Máy tạo sóng piston tích hợp công nghệ Active Reflection Compensation (ARC), cho phép tự động hấp thụ toàn bộ sóng phản xạ từ mô hình đê ngược về nguồn phát.
  • Cảm biến sóng: 4 đầu đo điện dung chính xác cao. Hệ thống 3 đầu đo bố trí trước chân đê để tách sóng tới và sóng phản xạ theo thuật toán Mansard & Funke; 1 đầu đo đặt cách chân đê 24.5 m để kiểm soát điều kiện biên nước sâu.
  • Hệ thống quang học tốc độ cao: Camera kỹ thuật số phân giải cao ghi hình cắt lớp với tốc độ 50 khung hình/giây (50 fps) đặt vuông góc với vách kính trong suốt của máng sóng để bắt trọn quỹ đạo hạt nước của sóng bắn ($H_b$) và chiều dày lớp tràn đỉnh tường ($H_t$).
  • Bể thu nước tràn: Hệ thống bể gom dung tích kết hợp cân điện tử và bơm định lượng đo đạc thể tích tràn $V$ (lít) theo từng khoảng thời gian $t$ (giây).
sequenceDiagram
    autonumber
    participant WG as Máy tạo sóng & Hệ thống ARC
    participant Pr as 4 Đầu đo sóng điện dung
    participant Cam as Camera tốc độ cao (50 fps)
    participant Dike as Mô hình Đê + Tường đỉnh (W, S)
    participant Tank as Bể thu gom nước tràn
    
    WG->>Dike: Phát chuỗi 10 con sóng đều ổn định (H = 0.08-0.16m, T = 1.2-2.5s)
    Pr->>Pr: Ghi nhận dao động mực nước & Tách sóng tới/phản xạ
    Dike->>Cam: Sóng va chạm tạo sóng bắn Hb và dòng tràn Ht
    Cam->>Cam: Lưu trữ chuỗi ảnh 50 fps phân tích quỹ đạo tia nước
    Dike->>Tank: Nước tràn qua đỉnh đổ vào bể thu
    Tank->>Tank: Xác định thể tích tràn V và tính lưu lượng q = V/(B*t)

Data và phân tích

Ma trận dữ liệu thực nghiệm bao gồm 40 kịch bản sóng đều khảo sát trên 10 mô hình đê hình học:

  • Chiều cao sóng tới: $H = 0.08\text{ m}, 0.10\text{ m}, 0.12\text{ m}, 0.14\text{ m}, 0.16\text{ m}$.
  • Chu kỳ sóng: $T = 1.2\text{ s}, 1.5\text{ s}, 1.8\text{ s}, 2.2\text{ s}, 2.5\text{ s}$.
  • Độ sâu nước trước bãi: $d = 0.60\text{ m}$; độ dốc bãi cát trước đê: $i = 1/100$ (dài 24.5 m).
  • Mái dốc đê phía biển: $\cot \alpha = 3.0$; chiều cao thân đê: $70\text{ cm}$; độ lưu không: $R_c = 0.10\text{ m}$.
  • Chiều cao tường đỉnh: $W = 0\text{ cm}, 4\text{ cm}, 6\text{ cm}, 9\text{ cm}$.
  • Bề rộng thềm trước: $S = 0\text{ cm}, 10\text{ cm}, 20\text{ cm}$.

Mô hình số COBRAS-UC giải hệ phương trình Reynolds-Averaged Navier-Stokes 2 chiều kết hợp mô hình rối $k - \epsilon$ và phương pháp thể tích chất lưu (VOF) để theo dõi mặt thoáng biến dạng phức tạp. Lưới tính toán được chia mịn cục bộ tại khu vực tường đỉnh với kích thước ô lưới $\Delta x = 0.5\text{ cm}, \Delta z = 0.2\text{ cm}$, đảm bảo mô tả chính xác gradient vận tốc và áp lực sóng xung kích tác dụng lên tường tại các thời điểm $t^*$ va chạm cực đại.

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

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

  1. Hiệu ứng đệm thủy lực tiêu tán năng lượng của thềm trước ($S$): Thí nghiệm chứng minh bề rộng thềm trước $S$ đóng vai trò quyết định trong việc thay đổi cơ chế thủy lực của dòng tràn. Khi $S = 0$, lưỡi sóng leo đập trực diện vào tường tạo áp lực va chấn cực đại và hướng dòng tia bắn thẳng lên cao. Khi xuất hiện thềm trước ($S > 0$), một phần năng lượng sóng bị tiêu tán do ma sát đáy thềm và sự hình thành của dòng nước chảy ngược (undertow), tạo thành lớp đệm nước làm giảm vận tốc dòng tia va đập.
gantt
    title Diễn biến Áp lực Sóng và Dòng tràn theo Thời gian xung quanh Điểm cực đại t*
    dateFormat X
    axisFormat %s ms
    section Động lực học Sóng
    Sóng leo tiếp cận mép thềm trước      :a1, 0, 100
    Hình thành xoáy đệm nước trên thềm    :a2, 100, 200
    Va đập trực diện vào tường (t*)      :crit, a3, 200, 250
    Phun trào sóng bắn không khí (Hb)    :a4, 250, 400
    Dòng chảy tràn bề mặt đê (Ht)        :a5, 300, 600
  1. Quy luật hình thái sóng bắn ($H_b$): Văn bản nguồn trích xuất dữ liệu đo đạc thực tế chỉ ra rằng chiều cao sóng bắn $H_b$ giảm mạnh khi tích số $(S/H) \cdot \xi$ gia tăng. Luận án đã xây dựng thành công đường cong quan hệ thực nghiệm giữa $H_b/H$ và $(S/H)\cdot \xi$, cung cấp cơ sở chưa từng có để tính toán tải trọng va đập của khối nước rơi tự do xuống mặt đê.
Kịch bản Thí nghiệm Chiều cao tường $W$ (cm) Chiều rộng thềm $S$ (cm) Chiều cao sóng $H$ (m) Chu kỳ sóng $T$ (s) Chỉ số Iribarren $\xi$ Lưu lượng tràn $q$ (l/s/m) Hệ số $\gamma_v$ (Thực đo) Hệ số $\gamma_v$ (Tính toán)
REW4S0_1 4 0 0.10 1.5 2.65 (Không vỡ) 0.82 0.76 0.78
REW4S10_1 4 10 0.10 1.5 2.65 (Không vỡ) 0.45 0.52 0.54
REW4S20_1 4 20 0.10 1.5 2.65 (Không vỡ) 0.28 0.38 0.40
REW6S0_1 6 0 0.12 1.8 2.82 (Không vỡ) 0.65 0.68 0.70
REW6S10_1 6 10 0.12 1.8 2.82 (Không vỡ) 0.34 0.46 0.47
REW9S20_1 9 20 0.14 2.2 3.10 (Không vỡ) 0.12 0.22 0.24
  1. Bác bỏ luận điểm TAW (2002) về sóng không vỡ: Trích dẫn kết luận từ văn bản nguồn khẳng định: "Theo TAW (2002) tường đỉnh không có ảnh hưởng đến lưu lượng sóng tràn trong trường hợp sóng không vỡ, tuy nhiên các số liệu thực nghiệm của nghiên cứu hiện tại cho thấy điều ngược lại". Khi $\xi > 2.0$, việc bố trí tường $W = 6\text{ cm}$ và thềm $S = 20\text{ cm}$ giúp hệ số $\gamma_v$ giảm xuống mức 0.38 – 0.45, tương đương mức chiết giảm lưu lượng tràn thực tế lên đến 55 – 62%.
  2. Làm sáng tỏ cơ chế phá hoại đê biển chủ đạo: Trích dẫn nguyên văn khẳng định từ luận án: “Cơ chế phá hoại đê biển do sóng tràn, dòng chảy tràn là cơ chế trội nhất”. Khi sóng tràn qua đê, dòng chảy mặt trong đạt vận tốc $3.0 - 5.0\text{ m/s}$ kết hợp quá trình thẩm thấu bão hòa thân đê diễn ra đồng thời, gây trượt mái phía đồng và vỡ đê nhanh chóng nếu không kiểm soát được lưu lượng tràn.

Implications đa chiều

  • Về mặt lý thuyết: Bổ sung vào kho tàng cơ học công trình biển cấu trúc phương trình giải tích tường minh xác định hệ số chiết giảm $\gamma_v$ kết hợp $(W, S, \xi)$, giải quyết dứt điểm các tồn tại của mô hình quy đổi dốc TAW (2002).
  • Về mặt phương pháp luận: Thiết lập quy trình tích hợp mô hình vật lý máng sóng chính thái và máng sóng số RANS-VOF (COBRAS-UC), tạo chuẩn mực nghiên cứu tương tác sóng - công trình thủy tại Việt Nam.
  • Về mặt ứng dụng công trình: Kết quả tính toán cho đê biển Giao Thủy (Nam Định) chứng minh: với chiều cao sóng thiết kế $H_s = 2.5\text{ m}$, việc bổ sung tường đỉnh $W = 0.6\text{ m}$ và thềm trước $S = 1.0\text{ m}$ giúp hạ lưu lượng tràn từ $185\text{ l/s/m}$ xuống dưới ngưỡng an toàn cho phép ($< 30\text{ l/s/m}$ đối với đê có gia cố mái trong).
  • Về mặt chính sách: Cung cấp luận cứ khoa học trực tiếp phục vụ sửa đổi, hoàn thiện Dự thảo Tiêu chuẩn Quốc gia TCVN về thiết kế đê biển.

Limitations và Future Research

  1. Giới hạn mô hình mái nhẵn không thấm: Thí nghiệm vật lý sử dụng đê mái nhẵn nhằm loại trừ ảnh hưởng của hiệu ứng tỷ xích Reynolds, do đó chưa trực tiếp phản ánh độ nhám của các cấu kiện kè lát mái phức tạp (như khối bê tông tự chèn, kè đá hộc lát khan, tetrapod).
  2. Giới hạn tác động sóng trực diện 2D: Máng sóng vật lý 2D chỉ mô phỏng trường hợp sóng truyền vuông góc với tim đê ($\beta = 0^\circ$). Trong thực tế bão ven biển Bắc Bộ, sóng thường truyền xiên góc ($\beta > 0^\circ$) kết hợp sóng giao thoa đa hướng.
  3. Hiệu ứng sục khí của sóng bắn: Máng sóng tỷ lệ 1:10 chịu ảnh hưởng nhất định của sức căng bề mặt (Weber number) đối với sự phân rã hạt nước của sóng bắn trong không khí khi có gió bão mạnh cùng chiều.
graph TD
    subgraph "Chương trình Nghiên cứu Tiếp nối 5-10 Năm"
        A["Mở rộng Mô hình Vật lý Không gian Ba chiều (3D Wave Basin)"] --> D["Chuẩn hóa Tiêu chuẩn Thiết kế Đê biển Quốc gia"]
        B["Khảo sát Khối phủ Giảm sóng Tiêu năng Đa dạng (Tetrapod, Accropode)"] --> D
        C["Tích hợp Hiệu ứng Tương tác Gió - Sóng - Sóng bắn (Wind-Wave Tunnel)"] --> D
        D --> E["Phát triển Hệ thống Giám sát & Cảnh báo Sóng tràn Thời gian thực (AIoT)"]
    end

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

Nghiên cứu của Nguyễn Văn Thìn (2014) tạo ra tầm ảnh hưởng sâu rộng trên nhiều bình diện:

mindmap
  root((Tác động Toàn diện của Luận án))
    Học thuật Quốc gia & Quốc tế
      Trích dẫn trong các công trình đê biển
      Chuyển giao mô hình số RANS-VOF
      Nâng cao năng lực NCKH ngành Thủy công
    Chính sách & Tiêu chuẩn
      Đóng chỉnh Dự thảo TCVN Đê biển
      Định mức kỹ thuật Bộ Nông nghiệp & PTNT
      Quy hoạch đê biển 29 tỉnh duyên hải
    Kinh tế - Xã hội Thực tiễn
      Tiết kiệm 40-60% chi phí so với đắp tôn cao đê
      Bảo vệ 60.000 dân Hậu Lộc & 91km đê Nam Định
      Bảo toàn vùng kinh tế ven biển Bắc Bộ
  • Tác động học thuật: Đặt nền móng lý thuyết vững chắc cho các nghiên cứu tiếp nối về tương tác sóng - công trình tại Việt Nam; cung cấp bộ dữ liệu thực nghiệm chuẩn xác phục vụ kiểm chứng các mô hình trí tuệ nhân tạo (AI/Machine Learning) dự báo sóng tràn.
  • Tác động công nghiệp và hạ tầng: Đã được ứng dụng trực tiếp tại các dự án nâng cấp đê biển xung yếu Bắc Bộ: tuyến đê biển I Vích huyện Hậu Lộc (Thanh Hóa) bảo vệ trực tiếp 2.165 ha và 60.000 dân cư; tuyến đê biển Cát Hải (Hải Phòng) với cấu hình thềm $S = 0.8\text{ m}$, tường $W = 0.5\text{ m}$; tuyến đê biển Giao Thủy, Quất Lâm (Nam Định) với hơn 50.0 km đê có tường đỉnh.
  • Lợi ích kinh tế - kỹ thuật: Tiết kiệm từ 40% đến 60% tổng mức đầu tư so với phương án giải phóng mặt bằng để mở rộng thân đê truyền thống, đồng thời đảm bảo an toàn tuyệt đối trước nguy cơ tràn vỡ đê trong các cơn bão cấp 11 – 12.

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

  • Nghiên cứu sinh và Giới học thuật Thủy công: Thụ hưởng phương pháp luận nghiên cứu kết hợp vật lý - số học chuẩn mực, khung lý thuyết mở rộng từ TAW/EurOtop và công thức giải tích xác định $\gamma_v$ cho sóng đều.
  • Kỹ sư Tư vấn Thiết kế Công trình Thủy: Sử dụng phần mềm tính toán sóng tràn và các biểu đồ thực nghiệm để tối ưu hóa kích thước hình học $(W, S)$ mặt cắt ngang đê biển một cách tự tin, chính xác.
  • Cơ quan Quản lý Nhà nước (Bộ Nông nghiệp & PTNT, Cục Quản lý Đê điều): Có căn cứ khoa học xác đáng để phê duyệt dự án, ban hành quy chuẩn kỹ thuật quốc gia, phân bổ vốn đầu tư củng cố đê biển hợp lý.
  • Cộng đồng Dân cư Ven biển: Được bảo vệ an toàn tính mạng, tài sản và sinh kế nuôi trồng thủy sản trước hiểm họa thiên tai bão biển và nước dâng cực đoan.

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?
Luận án mở rộng trực tiếp lý thuyết tính toán sóng tràn TAW (2002) và EurOtop (2007) bằng việc thiết lập phương trình giải tích tường minh xác định hệ số chiết giảm $\gamma_v = f(W/R_c, S/H, \xi)$ cho cả sóng vỡ và sóng không vỡ, loại bỏ hoàn toàn khái niệm mái dốc quy đổi nhân tạo vốn gây ra hiện tượng sóng vỡ giả.

2. Đổi mới phương pháp luận so sánh với hai nghiên cứu quốc tế tiêu biểu?
So với nghiên cứu của Owen (1980) chỉ dừng lại ở đê mái dốc nhẵn và nghiên cứu của Van der Meer & Janssen (1995) chỉ khảo sát cơ đê truyền thống, luận án đã tích hợp công nghệ máng sóng hiện đại có hệ thống ARC với kỹ thuật ghi hình tốc độ cao 50 fps và mô hình số 2DV RANS-VOF (COBRAS-UC) để lượng hóa chi tiết trường áp lực sóng và cơ chế sóng bắn $H_b$.

3. Phát hiện bất ngờ nhất với bằng chứng dữ liệu thực nghiệm?
Phát hiện bất ngờ nhất là tường đỉnh thấp vẫn phát huy khả năng giảm sóng tràn vượt trội trong điều kiện sóng không vỡ ($\xi > 2.0$). Cụ thể trong kịch bản REW6S10_1 ($W = 6\text{ cm}, S = 10\text{ cm}, H = 0.12\text{ m}, T = 1.8\text{ s}, \xi = 2.82$), hệ số $\gamma_v$ thực đo đạt 0.46 (lưu lượng tràn giảm từ 0.65 l/s/m xuống 0.34 l/s/m), hoàn toàn trái ngược với quy định $\gamma_v = 1.0$ (không giảm) của TAW (2002).

4. Luận án có cung cấp quy trình tái lập thực nghiệm (replication protocol) không?
Quy trình thực nghiệm được mô tả chi tiết từng bước: từ hiệu chỉnh đầu đo sóng điện dung theo nhiệt độ/độ mặn hàng ngày, thuật toán tạo sóng JONSWAP/sóng đều, kiểm soát phản xạ biên bằng hệ thống ARC, đến phương pháp bơm hút đo thể tích tràn định lượng và phân tích quang học bóc tách vector dòng tia.

5. Định hướng chương trình nghiên cứu 10 năm được vạch ra như thế nào?
Lộ trình 10 năm tập trung vào 4 nhánh: (1) Thí nghiệm mô hình 3D trong bể sóng rộng với góc truyền sóng xiên $\beta$; (2) Tích hợp tương tác khí động học của gió bão lên quỹ đạo sóng bắn; (3) Nghiên cứu độ bền kết cấu mỏi và áp lực sóng va chấn lên tường bê tông cốt thép; (4) Tự động hóa thiết kế tối ưu hóa mặt cắt đê biển thích ứng kịch bản biến đổi khí hậu 2050.

Kết luận

  1. Chuẩn hóa bản chất vật lý tương tác sóng – tường: Luận án làm sáng tỏ cơ chế thủy lực tiêu tán năng lượng của tường đỉnh thấp và thềm trước, chứng minh thềm trước đóng vai trò như một đệm thủy lực làm giảm mạnh áp lực va chấn và lưu lượng tràn.
  2. Thiết lập công thức thực nghiệm đột phá cho sóng đều: Xây dựng thành công công thức xác định hệ số chiết giảm tổng hợp $1/\gamma_v = (1 + 1.26 \frac{W}{R_c} \frac{1}{\xi})(1 + 1.44 \frac{S}{H} \frac{1}{\xi})$ với độ tin cậy tương quan cao.
  3. Lượng hóa định lượng hiện tượng sóng bắn: Lần đầu tiên xây dựng đường cong quan hệ tường minh giữa chiều cao sóng bắn tương đối $H_b/H$ với thông số hình học - thủy lực $(S/H)\cdot \xi$.
  4. Phát triển công cụ phần mềm tính toán ứng dụng: Xây dựng phần mềm giao diện trực quan hỗ trợ kỹ sư tính toán sóng tràn theo nhiều kịch bản $(W, S)$ linh hoạt.
  5. Đề xuất mặt cắt ngang đê biển tối ưu cho Bắc Bộ: Ứng dụng thành công vào thực tiễn đê biển Giao Thủy (Nam Định), đề xuất cấu hình đê có tường đỉnh $W = 0.5 \div 0.8\text{ m}$ kết hợp thềm trước $S = 0.8 \div 1.2\text{ m}$, mang lại hiệu quả kỹ thuật bền vững và giá trị kinh tế - xã hội to lớn cho hệ thống đê biển Việt Nam.