Tổng quan về luận án
Trong kỹ thuật xây dựng công trình thủy lợi – thủy điện, công trình tháo và hệ thống tiêu năng hạ lưu đóng vai trò quyết định đến sự an toàn công trình và ổn định lòng dẫn. Hiện tượng chuyển tiếp thủy lực từ trạng thái chảy xiết sang trạng thái chảy êm ($Fr > 1 \rightarrow Fr < 1$) mang mật độ động năng cục bộ cực lớn, đe dọa trực tiếp đến độ bền kết cấu và gây xói lở nghiêm trọng lòng sông. Luận án tiến sĩ kỹ thuật chuyên ngành Kỹ thuật Xây dựng Công trình Thủy (Mã số: 62580202) của NCS. Nguyễn Quốc Huy dưới sự hướng dẫn khoa học của PGS.TS Lê Văn Nghị tại Viện Khoa học Thủy lợi Việt Nam (2017) với đề tài: "Nghiên cứu một số đặc trưng thủy động lực học của dòng nối tiếp hỗn hợp mặt – đáy – ngập 3 xoáy sau bậc thụt" đã giải quyết một bài toán thủy lực biên động phức tạp chưa từng được mô hình hóa toàn diện.
Khoảng trống nghiên cứu (research gap) cốt lõi xuất phát từ thực tế các nghiên cứu kinh điển về nối tiếp sau bậc thụt trước đây chỉ khu trú ở ba nhóm hình thức:
- Bậc thụt phẳng hoặc góc hất nhỏ $\theta \le 15^\circ$ với $a/P = 0.25 \div 0.35$ (Cumin, 1948; Astaficheva, 1952; Ivanov, 2004; Kaverin, 2013).
- Bồn tiêu năng đa xoáy có góc hất lớn nhưng chiều cao bậc thụt gần như triệt tiêu $a = 0.05R$ (Peterka, 1958; Rajan & Shivashanakara Rao, 1980; Hager, 1992).
- Mũi hất phóng xa ở cao độ lớn với dòng phun tự do không ngập.
Dạng nối tiếp hỗn hợp mặt – đáy – ngập 3 xoáy (gọi tắt là dòng chảy phễu) hình thành từ bậc thụt có mũi hất cong ngược, góc hất lớn $\theta = 25^\circ \div 51^\circ$ và chiều cao bậc đáng kể ($a/P = 0.14 \div 0.46$) mới chỉ được Viện Nghiên cứu Thủy lực Nam Kinh (1985) khảo sát bán thực nghiệm đơn lẻ trên một vài công trình cụ thể, thiếu vắng các quy luật định lượng tổng quát.
Luận án thiết lập hệ thống câu hỏi nghiên cứu và giả thuyết khoa học chặt chẽ:
- Câu hỏi nghiên cứu 1 (RQ1): Điều kiện thủy lực và hình học biên nào (tỷ lệ $a/P$, góc hất $\theta$, số Froude $Fr$) xác lập ranh giới xuất hiện và duy trì ổn định chế độ dòng chảy phễu?
- Câu hỏi nghiên cứu 2 (RQ2): Cấu trúc không gian 3 vùng xoáy, chiều cao nước vồng ($h_v$) và chiều dài cuộn xoáy biến đổi theo quy luật toán học nào dưới tác động của mực nước hạ lưu $h_h$?
- Câu hỏi nghiên cứu 3 (RQ3): Phân bố lưu tốc trung bình, lưu tốc đáy lớn nhất ($U_m$) và hiệu suất tiêu tán năng lượng của dòng chảy phễu đạt mức định lượng bao nhiêu so với tiêu năng đáy truyền thống?
- Giả thuyết khoa học 1 (H1): Tồn tại một dải chiều sâu ngập hạ lưu xác định $\Delta h_{gh} = h_{max} - h_{min}$ phụ thuộc phi tuyến vào số Froude co hẹp $Fr_c$ và thông số hình học bậc hất ($a/P, \theta, R/h_c$), mà trong khoảng đó dòng chảy tự tổ chức thành hệ động lực 3 xoáy triệt tiêu động năng nội tại.
- Giả thuyết khoa học 2 (H2): Cấu trúc dòng chảy phễu nâng dòng chủ lưu mang lưu tốc cao lên bề mặt tự do, làm giảm triệt để lưu tốc đáy sát lòng dẫn hạ lưu so với nước nhảy đáy hoàn chỉnh, giúp thu hẹp đáng kể chiều dài gia cố sân sau.
Khung lý thuyết của công trình tích hợp giữa Thuyết tương tự thủy lực Froude, Định lý số $\pi$ Buckingham trong phân tích thứ nguyên và Lý thuyết dòng tia rối mở rộng (Turbulent Jet Theory). Công trình đạt đóng góp đột phá khi định lượng hóa toàn bộ miền làm việc của dòng chảy phễu qua 70+ kịch bản thực nghiệm, kiểm soát dải thông số Froude $Fr = 1.35 \div 4.5$, góc hất $\theta = 25^\circ \div 51^\circ$, và tỷ lệ bậc $a/P = 0.14 \div 0.46$. Nghiên cứu mở ra giải pháp công nghệ kinh tế – kỹ thuật vượt trội cho các công trình thủy lợi quy mô lớn tại Việt Nam như đập tràn Bản Mồng, Khe Bố, Hồi Xuân.
Literature Review và Positioning
Lịch sử nghiên cứu hiện tượng nối tiếp và tiêu năng sau công trình xả lũ ghi nhận nhiều trường phái học thuật quốc tế. Dòng chảy mặt gắn liền với bậc thụt được phát hiện lần đầu tiên bởi A. Sabaneyev trong thập niên 1920 tại đập thủy điện Vol-khov (Liên Xô cũ) (Ivanov, 2004). Tiếp đó, D. Cumin (1948) đã xây dựng bảng phân loại kinh điển gồm 9 trạng thái nối tiếp sau bậc thụt phẳng ($\theta = 0^\circ$), từ chảy đáy, chảy đáy ngập, chảy mặt, chảy mặt đáy hỗn hợp cho đến chảy tràn. M. Astaficheva (1952) phát triển hệ phương trình động lượng bổ sung với các biểu thức xác định cột nước áp suất dưới dòng cong $h_a$, thiết lập các công thức độ sâu giới hạn $h_{gh1}/h_k$ và $h_{gh2}/h_k$ cho bậc thụt có $\theta \le 15^\circ$ và $a/P \ge 0.2$.
Trường phái phương Tây, đại diện bởi Peterka (1958) tại Cục Khai hoang Hoa Kỳ (USBR), tập trung vào bồn tiêu năng đa xoáy (Roller Bucket - Bể tiêu năng số VII). Peterka khảo sát các mũi hất có $\theta = 8^\circ \div 16^\circ$, bán kính cong $R = (15.24 \div 45.72)\text{ cm}$, nhưng khống chế độ cao bậc thụt ở mức rất nhỏ ($a = 0.05R$). Hager (1992), Rajan & Shivashanakara Rao (1980), và Rajan et al. (1982) hoàn thiện phương pháp tính toán bồn tiêu năng dạng khối, thiết lập bán kính tối thiểu $R_{min}/h_c \cong 2.2 Fr_c^{0.5}$ cho phạm vi $Fr = 3 \div 10$.
Các mâu thuẫn học thuật nổi bật trong y văn bao gồm:
- Tranh luận về cơ chế ổn định của dòng mặt: Sabaneyev và Cumin nhận định dòng chảy mặt sau bậc thẳng đứng rất nhạy cảm với mực nước hạ lưu, dễ chuyển hóa đột ngột thành dòng chảy đáy gây xói sâu nguy hiểm. Ngược lại, Viện Nghiên cứu Thủy lực Nam Kinh (1985) và Thomas (1957) chứng minh rằng việc bố trí góc hất lớn ($\theta \ge 25^\circ$) kết hợp độ cong cung tròn sẽ ghìm giữ dòng xoáy cuộn bề mặt ổn định hơn nhiều nhờ lực ly tâm.
- Xung đột trong xác định chiều cao bậc nhỏ nhất $a_{min}$: Công thức của D. Cumin ($a_{min}/h_1 = \sqrt[3]{Fr_1^2 - 1}$), M. Skladnev ($a_{min}/h_k = 2.7(\sqrt[3]{Fr_1^2} - 1.6)$), và Kaverin (2013) cho các giá trị sai lệch đáng kể khi mở rộng góc hất $\theta > 15^\circ$.
| Công trình / Tác giả | Quốc gia | Loại kết cấu | Phạm vi thông số | Hạn chế cốt lõi |
|---|---|---|---|---|
| Peterka (1958) | Hoa Kỳ (USBR) | Bồn tiêu năng đa xoáy | $a = 0.05R$, $\theta = 8^\circ \div 16^\circ$ | Không áp dụng được cho bậc thụt cao ($a/P > 0.1$) |
| Astaficheva (1952) | Liên Xô | Bậc thụt đập tràn | $\theta \le 15^\circ$, $a/P \ge 0.2$ | Chỉ đúng với góc hất nhỏ, dạng phẳng |
| Kaverin (2012, 2013) | Nga | Bậc thụt thấp máng kính | $a/h_1 = 0.67 \div 4$, $\theta = 0^\circ$ | Bỏ qua phân bố áp suất phi thủy tĩnh trong luồng cong |
| Nam Kinh (1985) | Trung Quốc | Mũi hất phễu liên tục | $Fr > 3$, nghiên cứu theo dự án | Thiếu hệ số tương quan cho dải $a/P$ biến thiên rộng |
| Luận án (2017) | Việt Nam | Mũi hất cong bậc cao | $\theta = 25^\circ \div 51^\circ$, $a/P = 0.14 \div 0.46$ | Khái quát hóa toàn diện bằng thực nghiệm đa biến |
So sánh với các công trình quốc tế quy chuẩn, nghiên cứu của Nguyễn Quốc Huy đã định vị chính xác khoảng trống: kết hợp ưu điểm tạo xoáy bề mặt của bồn tiêu năng Peterka với năng lực phân tán thế năng của bậc thụt cao kiểu Nga, thiết lập mô hình thủy lực định lượng cho dòng chảy phễu trong dải $a/P = 0.14 \div 0.46$ và $\theta = 25^\circ \div 51^\circ$.
Đó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 Lý thuyết nước nhảy mặt kinh điển của A. Sabaneyev và M. Astaficheva sang miền thủy lực có góc lệch lớn và dòng gia tốc uốn cong:
- Chuẩn hóa phương trình biến thiên động lượng phi tuyến: Bổ sung lực gia tốc ly tâm của dòng chảy cong tại mũi bậc vào thành phần phản lực đáy $ha$, làm biến đổi phân bố áp suất từ dạng thủy tĩnh tuyến tính sang đường cong hyperbol thực tế.
- Khái quát hóa mô hình 3 miền xoáy tương hỗ: Thay vì chỉ xét một xoáy mặt đơn lẻ như nước nhảy mặt truyền thống, mô hình lý thuyết định nghĩa cấu trúc động lực học "3 xoáy 1 sóng":
- Xoáy 1 ($L_1$): Xoáy cuộn đáy dưới mũi phễu do phản lực tách dòng.
- Xoáy 2 ($L_2$): Xoáy cuộn bề mặt uốn ngược ngược chiều kim đồng hồ mang năng lượng phân tán chính.
- Xoáy 3 ($L_3$): Xoáy ngập sau đỉnh nước vồng liên hiệp với dòng chảy hạ lưu.
- Sóng dâng nước vồng ($h_v$): Đỉnh sóng nâng luồng chủ lưu khỏi đáy kênh.
Mặt nước thượng lưu
\~~~~~~\
\ (Mặt tràn)
\
\_________ Xoáy mặt (L2) <--- Sóng vồng (hv)
/ ~~~~~~~~~~~~~~~~~~~~~ /\~~~~~~~ Mặt nước hạ lưu (hh)
/ / \ / \~~~~~~
(Mũi hất) / \/
[θ, R] / Chủ lưu mặt Xoáy ngập (L3)
| / =================> ------------>
|____/ Xoáy đáy (L1)
| a (Dưới mũi hất)
|
==================+=================================================== Đáy kênh hạ lưu
- Thiết lập các mệnh đề toán học về ranh giới trạng thái:
- Mệnh đề 1: Độ sâu giới hạn dưới $h_{min}$ là hàm đơn điệu tăng theo số Froude co hẹp $Fr_c$ và tỷ số chiều cao bậc $a/P$.
- Mệnh đề 2: Độ sâu giới hạn trên $h_{max}$ tỷ lệ thuận với bán kính cong $R/h_c$ và góc hất $\theta$, xác định ranh giới chuyển hóa từ dòng chảy phễu (HT3) sang dòng chảy đáy ngập (HT5).
Khung phân tích độc đáo
Khung phân tích của luận án tích hợp đa chiều 3 nền tảng lý thuyết: Lý thuyết dòng tia rối biên (Schlichting Boundary Layer Theory), Phương trình tích phân động lượng Karman Navier-Stokes, và Định lý hàm số $\pi$ Buckingham:
$$\Pi_1 = f(\Pi_2, \Pi_3, \Pi_4, \Pi_5) \Longleftrightarrow \frac{h_h}{h_k} = f\left(\frac{a}{P}, \frac{a}{H}, \theta, Fr_c, \frac{R}{h_c}\right)$$
Điều kiện biên xác định (Boundary Conditions):
- Bài toán phẳng hai chiều (2D), dòng chảy không đều biến đổi dần.
- Dòng tự do qua đỉnh tràn không bị cản trở bởi cửa van điều tiết.
- Số Froude luồng tới: $Fr = 1.35 \div 4.5$.
- Tỷ lệ bậc thụt: $a/P = 0.14 \div 0.46$.
- Góc hất mũi bậc cong liên tục: $\theta = 25^\circ \div 51^\circ$, đỉnh mũi hất ngập hoàn toàn dưới mực nước hạ lưu.
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ủ triết lý khoa học thực chứng (Positivism), kết hợp chặt chẽ giữa mô hình hóa giải tích thứ nguyên và quy hoạch thực nghiệm vật lý quy mô lớn. Nghiên cứu sử dụng mô hình vật lý thủy lực bài toán phẳng tỷ lệ hình học xác định theo Tiêu chuẩn Tương tự Froude ($\lambda_L = \lambda_v^2 = \lambda_a$), bảo toàn đồng dạng trường trọng lực và quán tính, đồng thời kiểm soát nghiêm ngặt số Reynolds mô hình ($Re_m > Re_{gh} = 10^4$) để loại trừ hoàn toàn ảnh hưởng của hiệu ứng nhớt (scale effect).
Quy trình nghiên cứu rigorous
Quy trình thu thập và xử lý số liệu được triển khai tại Phòng Thí nghiệm Trọng điểm Quốc gia về Động lực học Sông Biển:
- Hệ thống thiết bị: Máng kính chuyên dụng dài $13.26\text{ m}$, rộng $0.60\text{ m}$, cao $1.0\text{ m}$ với đáy máng và thành bên bằng kính hữu cơ trong suốt, cho phép quan sát trực tiếp chuyển động laser bọt khí và trường xoáy.
- Dụng cụ đo đạc chuẩn xác cao:
- Ống đo áp lực vi sai điện tử và cảm biến biến áp đa kênh đo mạch động áp suất đáy và mặt đập với tần số lấy mẫu cao.
- Vận tốc dòng chảy được đo bằng máy đo lưu tốc cánh quạt vi cơ điện tử kết hợp ống Pitot-Prandtl chuẩn hóa, xác định trường vector vận tốc tại 15+ mặt cắt dọc từ mũi bậc đến hạ lưu.
- Kim đo mực nước điện tử có độ chính xác $\pm 0.1\text{ mm}$.
- Quy trình kiểm định chéo (Triangulation): Kết hợp kiểm chứng tam giác giữa số liệu thực nghiệm mô hình vật lý, các công thức giải tích kinh điển của Nga - Trung Quốc, và phương pháp thống kê phân tích tương quan đa biến.
Data và phân tích
Hệ thống dữ liệu gồm hơn 70 chuỗi thí nghiệm đầy đủ biến thiên các thông số $a/P$, $\theta$, $q$, $H_0$, $h_h$. Phương pháp Holdout được áp dụng để đánh giá mô hình hồi quy: 80% tập số liệu dùng để thiết lập hàm toán học và 20% tập số liệu độc lập dùng để kiểm định độ tin cậy.
Các chỉ tiêu thống kê kiểm định:
- Hệ số tương quan bội: $R^2 \ge 0.92 \div 0.98$.
- Sai số chuẩn (Standard Error of Estimate): $SE < 0.05$.
- Sai số tương đối giữa tính toán và thực đo: $\delta < 8.5%$.
- Hệ số tin cậy kiểm định Fisher $F > F_{critical}$ ở mức ý nghĩa $p < 0.001$.
Phát hiện đột phá và implications
Những phát hiện then chốt
- Xác lập ranh giới định lượng hình thành dòng chảy phễu: Tìm ra quy luật tương quan của chiều sâu hạ lưu nhỏ nhất ($h_{min}$) và lớn nhất ($h_{max}$) xuất hiện dòng chảy phễu qua các biểu thức kinh nghiệm không thứ nguyên: $$\frac{h_{min}}{h_k} = f\left(\frac{a}{P}, \frac{a}{H}, \theta\right) = 0.84\left(\frac{a}{h_k}\right)^{0.32} \left(\frac{a}{P}\right)^{-0.45} (\sin\theta)^{0.65}$$ $$\frac{h_{max}}{h_k} = 1.16\left(\frac{a}{h_k}\right)^{0.28} \left(\frac{a}{P}\right)^{-0.52} (\sin\theta)^{0.48} + 2.38$$
- Hiện tượng nước vồng bề mặt ($h_v$): Luận án chứng minh rằng chiều cao cột nước vồng $h_v$ đóng vai trò là độ sâu liên hiệp thứ hai của xoáy mặt ($h_{22}$). Chiều cao tương đối $h_v/h_k$ đạt cực đại khi góc hất $\theta = 40^\circ \div 45^\circ$, giải phóng toàn bộ năng lượng dư thừa lên không gian mặt nước thay vì tác động vào đáy kênh.
- Quy luật suy giảm đột biến của lưu tốc đáy: Tại vị trí $x \approx (3 \div 5)a$ sau bậc thụt, lưu tốc đáy $U_m$ trong chế độ chảy phễu giảm từ $60%$ đến $75%$ so với lưu tốc đáy của nước nhảy đáy thông thường tại cùng mặt cắt: $$U_{max,day} \le (0.25 \div 0.35) V_c$$
- Hiệu suất tiêu tán năng lượng vượt trội: Hiệu quả tiêu hao năng lượng tương đối của dòng chảy phễu đạt từ $40%$ đến $65%$ ngay trong phạm vi khu xoáy ($L \le 6a$), làm triệt tiêu xung lực đáy và ngăn ngừa hiện tượng xâm thực mài mòn bê tông.
Implications đa chiều
- Về mặt lý thuyết: Hoàn thiện lý thuyết thủy lực nối tiếp công trình tháo, làm cầu nối dung hòa giữa lý thuyết bồn tiêu năng Peterka và lý thuyết nước nhảy mặt Sabaneyev - Astaficheva.
- Về mặt phương pháp luận: Chuẩn hóa quy trình ứng dụng quy hoạch thực nghiệm kết hợp lý thuyết Buckingham trong việc giải quyết các bài toán thủy động lực học dòng xoáy đa pha phức tạp.
- Về mặt thực tiễn kỹ thuật: Cung cấp bộ công cụ tính toán thủy lực cho phép các kỹ sư thiết kế thay thế kết cấu bể tiêu năng đáy sâu tốn kém bằng kết cấu bậc thụt mũi hất cong treo cao, tiết kiệm chi phí xây dựng.
- Về mặt chính sách và quản lý: Định hình cơ sở khoa học để Bộ Nông nghiệp & PTNT và Bộ Công Thương ban hành các chỉ dẫn kỹ thuật cập nhật tiêu chuẩn thiết kế công trình thủy lợi - thủy điện vùng địa chất đá cứng và bán phong hóa.
Limitations và Future Research
- Giới hạn bài toán phẳng hai chiều (2D): Nghiên cứu được thực hiện trong máng hẹp với biên lòng dẫn đối xứng, chưa xét đến các hiệu ứng dòng chảy không gian 3D phức tạp phát sinh từ trụ pin đập tràn, dòng chảy xiên hoặc co hẹp dòng chảy từ các cửa xả cục bộ.
- Giới hạn chế độ dòng chảy không điều tiết: Toàn bộ công thức thực nghiệm được thiết lập cho dòng chảy tràn tự do qua đỉnh; chưa bao quát các kịch bản đóng mở cửa van cung lệch tâm gây dòng chảy ngầm áp lực cao.
- Giới hạn dải thông số hình học: Các quan hệ toán học đạt độ tin cậy tuyệt đối trong biên: $Fr_c = 1.35 \div 4.5$, $a/P = 0.14 \div 0.46$, $\theta = 25^\circ \div 51^\circ$. Khi công trình vượt ra ngoài các cận này, cần có các hiệu chỉnh bổ sung.
Chương trình nghiên cứu tiếp nối (Future Research Agenda):
- Phát triển mô phỏng số trị 3D sử dụng phần mềm Flow-3D và các mô hình rối tiên tiến ($k-\omega\text{ SST}$, Large Eddy Simulation - LES) để kiểm chứng trường mạch động áp suất 3 chiều.
- Mở rộng nghiên cứu cho tràn có cửa van phân phối lưu lượng không đều theo phương ngang đập.
- Đánh giá cơ chế tương tác giữa dòng chảy phễu và quá trình bồi xói bùn cát đáy vùng hạ lưu lòng dẫn mềm.
Tác động và ảnh hưởng
Nghiên cứu của tác giả Nguyễn Quốc Huy tạo tiếng vang lớn trong cộng đồng khoa học thủy lực và kỹ thuật công trình thủy:
- Tác động học thuật: Cung cấp nguồn trích dẫn nền tảng cho các đề tài nghiên cứu cao học, nghiên cứu sinh và các bài báo chuyên ngành trên Tạp chí Khoa học và Công nghệ Thủy lợi, Tạp chí Nông nghiệp & PTNT.
- Đột phá kinh tế - kỹ thuật ngành xây dựng thủy lợi: Ứng dụng thành công kết cấu bậc thụt dòng chảy phễu giúp cắt giảm từ $15%$ đến $30%$ khối lượng bê tông cốt thép so với phương án đào sâu làm bể tiêu năng đáy truyền thống, rút ngắn thời gian thi công hố móng mùa kiệt.
- Bảo vệ môi trường sinh thái hạ du: Nhờ phân tán dòng chảy lên bề mặt và giảm vận tốc sát đáy, chế độ chảy phễu hạn chế tối đa hiện tượng sạt lở bờ sông và phá hủy thảm thực vật thủy sinh hạ lưu đập.
Đối tượng hưởng lợi
- Nghiên cứu sinh và Giảng viên ngành Công trình Thủy: Tiếp cận hệ thống dữ liệu thực nghiệm chuẩn tắc, phương pháp phân tích thứ nguyên và các mô hình hồi quy đa biến có độ chính xác cao.
- Các Công ty Tư vấn Thiết kế Thủy lợi – Thủy điện (HEC, PECC): Sở hữu quy trình 7 bước tính toán định lượng kết cấu bậc thụt mũi cong để lựa chọn phương án tiêu năng tối ưu cho các dự án đập tràn có cột nước trung bình và lưu lượng lớn.
- Chủ đầu tư các công trình hồ chứa: Giảm thiểu rủi ro sự cố xói lở chân đập, tiết kiệm hàng chục tỷ đồng chi phí bảo trì, gia cố lòng dẫn sau nhiều năm vận hành xả lũ.
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 lý thuyết độc đáo nhất là việc mở rộng Lý thuyết Nước nhảy mặt Sabaneyev - Astaficheva sang miền thủy lực biên uốn cong góc lớn ($\theta = 25^\circ \div 51^\circ$). Tác giả đã mô hình hóa thành công cấu trúc "3 xoáy 1 sóng", xác định bản chất của chiều cao nước vồng $h_v$ chính là chiều sâu liên hiệp thứ hai của xoáy mặt ($h_{22}$), làm sáng tỏ cơ chế tự triệt tiêu động năng trong lòng cuộn xoáy phễu.
2. Điểm mới về phương pháp luận so với các nghiên cứu tiền nhiệm quốc tế?
So với nghiên cứu của Peterka (1958) chỉ giới hạn $a = 0.05R$ và Kaverin (2013) chỉ áp dụng bậc thẳng đứng $a/h_1 < 4$, luận án đã:
- Kết hợp hoàn hảo giữa Lý thuyết thứ nguyên Buckingham và Quy hoạch thực nghiệm đa biến DoE.
- Áp dụng phương pháp kiểm định Holdout (tách 80% train / 20% test), bảo đảm tính khách quan và triệt tiêu hiện tượng quá khớp (overfitting) trong xây dựng công thức thực nghiệm thủy lực.
3. Phát hiện bất ngờ và phản trực giác nhất từ số liệu thực nghiệm?
Trong dòng chảy phễu, khi tăng góc hất từ $25^\circ$ lên $45^\circ$, chiều dài khu vực tiêu năng và vận tốc đáy không tăng theo quán tính văng xa mà ngược lại suy giảm nhanh hơn. Lực ly tâm của mũi cong lớn đã ép luồng nước uốn ngược cực mạnh tạo phễu nước sâu, hút ngược nước hạ lưu vào vùng đáy chân đập và triệt tiêu vận tốc đáy ($U_m$ giảm tới $75%$).
4. Luận án có cung cấp quy trình nhân bản/ứng dụng (Replication Protocol) không?
Có. Luận án xây dựng quy trình tính toán thiết kế chuẩn mực gồm 7 bước tuần tự tại Chương 4:
- Xác định thông số công trình ($P, Q, q, H_0, h_h$).
- Tính toán thông số động học dòng chảy ($h_c, V_c, Fr_c, h_k$).
- Xác định bán kính cong mũi hất ($R$).
- Tính toán cao độ đặt đáy mũi hất ($a_0$) và chiều cao bậc ($a$).
- Kiểm tra điều kiện xuất hiện dòng chảy phễu ($h_{min} \le h_h \le h_{max}$).
- Xác định kích thước hình học ($h_v, L_1, L_2, L_3$).
- Tính toán phân bố lưu tốc đáy và kiểm tra ổn định xói lòng dẫn.
5. Kiểm chứng thực tế trên công trình cụ thể nào?
Quy trình đã được áp dụng kiểm chứng trực tiếp trên công trình Đập tràn xả lũ hồ chứa Bản Mồng (Nghệ An). Với các thông số $q = 18\text{ m}^3/\text{s/m}$, góc hất tối ưu $\theta = 40^\circ$, tỷ lệ $a/P = 0.39$, kết quả tính toán theo công thức của luận án cho thấy dòng chảy hoàn toàn duy trì ổn định ở chế độ chảy phễu, giúp công trình tiết giảm đáng kể chiều dài sân tiêu năng so với thiết kế đáy ngập truyền thống.
Kết luận
Luận án tiến sĩ của NCS. Nguyễn Quốc Huy là một công trình nghiên cứu khoa học công phu, nghiêm túc và có giá trị học thuật cũng như thực tiễn đặc biệt xuất sắc:
- Khái quát hóa toàn diện điều kiện thủy động lực học của dòng nối tiếp hỗn hợp mặt – đáy – ngập 3 xoáy sau bậc thụt mũi cong, lấp đầy khoảng trống tồn tại nhiều thập kỷ giữa trường phái bồn tiêu năng USBR (Hoa Kỳ) và trường phái nước nhảy mặt Liên Xô.
- Xây dựng hệ thống công thức thực nghiệm có độ tin cậy cao ($R^2 > 0.92, \delta < 8.5%$) xác định chính xác các đại lượng $h_{min}, h_{max}, h_v, L_r, U_m$ phục vụ thiết kế công trình thủy lợi.
- Phát hiện quy luật tiêu tán năng lượng nội tại đạt hiệu suất $40% \div 65%$ và cơ chế nâng dòng chủ lưu giúp triệt tiêu lưu tốc đáy sát lòng dẫn hạ lưu.
- Chuẩn hóa quy trình thiết kế thủy lực 7 bước rõ ràng, logic, có khả năng chuyển giao và áp dụng trực tiếp vào thực tiễn kỹ thuật công trình thủy lợi - thủy điện tại Việt Nam.
- Mở ra các hướng nghiên cứu liên ngành mới về mô hình hóa số trị 3D dòng xoáy đa pha và động lực học bồi xói lòng dẫn phức tạp sau công trình xả nước.
- Mang lại hiệu quả kinh tế - xã hội định lượng rõ rệt, đóng góp thiết thực vào sự nghiệp bảo đảm an toàn hồ đập và phát triển bền vững cơ sở hạ tầng tài nguyên nước quốc gia.