Tổng quan về luận án

Sự chuyển dịch cơ cấu năng lượng và áp lực đảm bảo an ninh năng lượng quốc gia đặt ra yêu cầu cấp thiết về việc phát triển các nhà máy thủy điện tích năng (NMTĐTN) nhằm phủ đỉnh biểu đồ phụ tải và điều tiết hệ thống điện. Theo Quyết định số 428/QĐ-TTg ngày 18/03/2016 của Thủ tướng Chính phủ phê duyệt điều chỉnh Quy hoạch phát triển điện lực quốc gia giai đoạn 2011 - 2020 có xét đến năm 2030, cơ cấu nguồn cho thấy "tỷ lệ công suất các công trình thủy điện trong hệ thống giảm xuống chỉ còn khoảng 20-25%... bảy phương án TĐTN có khả năng được xây dựng với tổng công suất lắp máy là 2100MW", điển hình như các tổ máy Bác Ái (4 x 300 MW) và Đông Phù Yên (3 x 300 MW) với cột nước làm việc lên tới gần 500 m. Trong các công trình này, máy thủy lực thuận nghịch bơm – tuabin (Pump as Turbine - PaT) là thiết bị trung tâm đóng vai trò quyết định đến hiệu quả kinh tế và độ tin cậy vận hành.

Khoảng trống nghiên cứu (research gap) cốt lõi nằm ở sự mâu thuẫn thủy lực nội tại giữa hai chế độ làm việc: dòng chảy trong chế độ bơm là dòng khuếch tán (loe) chịu gradient áp suất ngược, trong khi chế độ tuabin là dòng thu hẹp (côn) với gradient áp suất thuận. Các lý thuyết thiết kế kinh điển từ Stepanoff (1957), Childs (1962) hay Sharma (1985) chỉ dự báo đặc tính PaT dựa trên tương quan thống kê thực nghiệm đơn lẻ với sai số vượt chuẩn BS516 lên tới 20% - 40%, hoàn toàn thiếu một khung lý thuyết giải tích toàn diện mô tả phân bố tổn thất cục bộ và hiệu ứng trượt dòng cho máy có số vòng quay đặc trưng thấp ($n_s = 90 \div 150$). Luận án tiến sĩ kỹ thuật của NCS. Nguyễn Thị Nhớ với đề tài "Nghiên cứu ảnh hưởng của một số thông số hình học của bánh công tác đến đặc tính làm việc của máy thuận nghịch bơm–tuabin ns thấp" (chuyên ngành Máy và Tự động thủy khí, Trường Đại học Bách Khoa Hà Nội, 2019) đã giải quyết triệt để bài toán này.

Hệ thống câu hỏi và giả thuyết nghiên cứu được xác lập tường minh:

  • Câu hỏi nghiên cứu 1 (RQ1): Cơ chế suy giảm năng lượng và quy luật hình thành các thành phần tổn thất thủy lực qua 5 phân vùng dòng chảy trong tổ máy PaT $n_s$ thấp biến thiên như thế nào giữa hai chiều dòng chảy ngược nhau?
  • Câu hỏi nghiên cứu 2 (RQ2): Quy luật tương tác của các thông số hình học then chốt gồm đường kính ngoài ($D_2$), góc đặt cánh vào/ra ($\beta_{1B}, \beta_{2B}$), số cánh ($Z$) và góc tới ($\delta$) chi phối trực tiếp đến hiệu suất đỉnh ($\eta_{max}$), điểm làm việc tối ưu (BEP) và chất lượng xâm thực ra sao?
  • Giả thuyết nghiên cứu 1 (H1): Việc tích hợp giải tích phương trình Euler có hiệu chỉnh hệ số trượt Stodola – Gülich ($\gamma$) và hệ số chèn dòng ($\tau$) với mô hình hóa 10 thành phần tổn thất thủy lực cho phép dự báo chính xác đường đặc tính $H(Q), P(Q), \eta(Q)$ của PaT với độ lệch dưới 5% so với thực nghiệm.
  • Giả thuyết nghiên cứu 2 (H2): Biên dạng cánh dạng một cung cong đối xứng (symmetrical single-arc profile) sẽ tối ưu hóa trường vận tốc và áp suất, triệt tiêu xoáy quẩn tách dòng ở mép cánh, giúp nâng cao hiệu suất tổ máy từ 1.0% đến 2.5% trên toàn dải vận hành.

Khung lý thuyết của luận án tích hợp chặt chẽ giữa lý thuyết dòng xoáy cánh dẫn Euler, động lực học dòng chảy nhớt Navier-Stokes, lý thuyết trượt dòng Stodola - Gülich và các mô hình tổn thất cấu trúc Jet-Wake (Bing et al., Djebedjian). Đối tượng nghiên cứu tập trung vào mô hình PaT $n_s = 104$ điển hình cho các dự án tích năng cột nước cao tại Việt Nam, kết hợp khảo nghiệm thực chứng tại Công ty Cổ phần Chế tạo Bơm Hải Dương với tốc độ quay đồng bộ $n = 600\text{ v/p}$ và ba chế độ mở cánh hướng ($a_0 = 5\text{ mm}, 10\text{ mm}, 15\text{ mm}$).

Literature Review và Positioning

Nghiên cứu về máy thủy lực thuận nghịch ghi nhận cột mốc từ năm 1931 khi Thoma và Kittredge phát hiện "máy bơm hoàn toàn có thể vận hành một cách hiệu quả trong chế độ của tuabin mà không có bất cứ một vấn đề cơ khí nào xảy ra". Trải qua nhiều thập kỷ, các công trình quốc tế phát triển theo hai dòng tiếp cận chính:

Dòng tiếp cận thực nghiệm thống kê nhằm thiết lập tỷ số chuyển đổi cột nước $h = H_T/H_P$, lưu lượng $q = Q_T/Q_P$ tại điểm BEP. Stepanoff (1957) và Childs (1962) đề xuất mô hình dựa trên hiệu suất bơm đỉnh ($\eta_P$), trong khi Grover (1980) và Hergt (1982) xây dựng công thức phụ thuộc vào số vòng quay đặc trưng $n_s$. Schmiedl (1988) và Alatorre-Frenk (1994) đưa vào các hàm bậc hai và hệ số hiệu chỉnh không gian. Tuy nhiên, khi Williams (1994) tiến hành kiểm chứng độc lập trên 35 mẫu bơm thực tế theo tiêu chuẩn sai số elip Anh Quốc BS516 ($(a/0.3)^2 + (b/0.1)^2 \le 1$), kết quả chỉ ra rằng mô hình của Grover có tới 81% số điểm vượt ngưỡng sai số cho phép ($C > 1$, giá trị trung bình $C = 1.333$), mô hình Childs và Schmiedl có 40% điểm sai lệch, và ngay cả mô hình tối ưu nhất của Sharma (1985) vẫn có 20% mẫu không đạt yêu cầu ($C = 0.733$). Điều này chứng minh các công thức thống kê đơn biến hoàn toàn bất lực trong việc bao quát các dòng máy có kết cấu hình học chuyên biệt.

Dòng tiếp cận giải tích thủy lực và tối ưu hóa số tập trung giải mã cấu trúc dòng chảy chi tiết. Gülich (2003, 2010) và Chapallaz (1992) đã tiên phong phân tích cơ chế trượt dòng và chèn dòng, song các công thức giải tích chỉ khu biệt cho bơm ly tâm thông thường. Derakhshan et al. (2008) ứng dụng phần mềm FINETURBO V.7 giải phương trình Navier-Stokes 3D với lưới AUTOGRID 5 tái thiết kế cánh bơm $n_s = 23.5$, nâng hiệu suất tuabin thêm 3% - 5%. Xuhe et al. (2014) tối ưu hóa đa mục tiêu với mặt Pareto trên 1000 cấu hình cánh cho NMTĐTN cột nước $H = 298\text{ m}$, công suất $P = 270\text{ MW}$, $n = 300\text{ v/p}$, đạt hiệu suất tuabin $\eta_T = 91%$ và bơm $\eta_P = 90.8%$. Yang et al. (2012) khảo nghiệm thực nghiệm trên 3 mẫu bơm ly tâm và chứng minh khi giảm đường kính $D_2$ từ 255 mm xuống 215 mm, hiệu suất PaT suy giảm 4.11%, đường đặc tính $Q(H)$ dốc nhanh và $Q(P)$ chuyển phẳng.

Các tranh luận học thuật xoay quanh việc giải quyết mâu thuẫn hình học: Theo Stelzer và Walters, "với cùng các thông số đầu vào, trong khi đường kính tại mép vào bơm D1 trong hai mô hình là xấp xỉ nhau thì D2 trong tuabin lại khá nhỏ hơn so với bơm, D2T ≈ 0,83D2P". Đồng thời, số cánh tối ưu của bơm ly tâm chỉ từ $Z = 6 \div 9$, trong khi tuabin Francis đòi hỏi $Z = 13 \div 21$. Luận án định vị chính xác tại điểm giao thoa này, khắc phục khoảng trống bằng việc xây dựng mô hình giải tích tổn thất 5 phân vùng kết hợp thuật toán tối ưu hóa biên dạng cánh một cung cong đối xứng cho máy PaT $n_s = 90 \div 150$, tiến xa hơn các công trình của Li (2014) và Olimstad et al. (2012) vốn chỉ dựa thuần túy vào CFD thuần số học mà thiếu mô hình lý thuyết định lượng làm gốc.

Đó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 căn bản lý thuyết máy cánh dẫn Euler thông qua việc thiết lập hệ phương trình cột nước lý thuyết có xét đến tương tác trượt dòng phi tuyến và độ dày cánh hữu hạn cho cả hai chiều thuận nghịch:

Ở chế độ bơm, cột áp lý thuyết $H_{th\infty, P}$ được hiệu chỉnh toàn diện qua hệ số trượt $\gamma_2$ và hệ số chèn dòng $\tau_2$: $$H_{th,P} = \frac{\gamma_2 u_2^2}{g} - \frac{u_2 Q_P \tau_2}{g A_{2m} \tan \beta_{2B}} - \frac{u_1 Q_P \tau_1}{g A_{1m} \tan \alpha_1}$$

Ở chế độ tuabin, thành phần vận tốc tiếp tuyến $c_{2u}$ ở lối vào được liên kết trực tiếp với thông số hình học của hệ thống cánh hướng qua tỷ số bán kính $R_3/R_2$ và tiết diện lưu thông $A_{3m}$: $$H_{th,T} = \frac{R_3 \cot \alpha_{3B} Q_T u_2}{R_2 g A_{3m}} - \frac{u_1^2 \sin \beta_{1B}}{g Z} + \frac{u_1^2}{g} - \frac{c_{1m} \cot \beta_{1B} u_1}{g}$$

Mô hình lý thuyết đã thách thức quan niệm truyền thống cho rằng ảnh hưởng của hiện tượng trượt dòng là đồng nhất ở hai chế độ, chứng minh bằng giải tích rằng hiệu ứng trượt $\gamma$ trong chế độ tuabin có thể triệt tiêu nhờ sự định hướng cưỡng bức của vành cánh hướng tĩnh, trong khi chế độ bơm chịu tổn thất trượt phụ thuộc hàm mũ theo $\beta_{2B}$ và số cánh $Z$ ($\gamma = f_1 [1 - k_w (\sin \beta_{2B})^{0.7} / Z^{0.7}]$).

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

Khung phân tích của luận án tích hợp ma trận 10 thành phần tổn thất năng lượng phân bổ trên 5 không gian chức năng độc lập:

  1. Buồng xoắn ($h_{cas}$): Tổn thất ma sát màng và xoáy thứ cấp theo diện tích ướt $\Delta A$ và hệ số ma sát $c_f$.
  2. Vành cánh hướng ($h_{van}$): Kết hợp tổn thất ma sát $h_{fv}$ và tổn thất va đập dòng vào $h_{sv}$ do sai lệch góc dòng chảy $\alpha_2$.
  3. Mép vào bánh công tác ($h_{in}, h_{sc}, h_{rec}$): Mô hình hóa tổn thất va đập đầu cánh, co hẹp đột ngột và tổn thất dòng quẩn ngược theo tỷ lệ $(1 - Q/Q_{BEP})^{2.5}$ khi non tải.
  4. Máng cánh bánh công tác ($h_{fr}, h_{bend}, h_{dif}, h_{cir}$): Xác định tổn thất ma sát bề mặt theo bán kính thủy lực trung bình $D_H$, tổn thất uốn cong góc lượn $\delta_0$, tổn thất phân tán cấu trúc Jet-Wake khi tỷ số vận tốc tương đối $(w_0/w_2)_{crit} \ge 1.4$, và tổn thất xoáy thứ cấp.
  5. Cửa ra và ống hút/ống xả ($h_{se}, h_{ge}, h_{dr}$): Mô hình hóa quá trình khuếch tán dòng loe mở rộng đột ngột hoặc thu hẹp dòng côn.

Điều kiện biên phân tích được xác lập chặt chẽ: chất lỏng là nước đơn pha không nén được ở $20^\circ\text{C}$, khối lượng riêng $\rho = 998.2\text{ kg/m}^3$, độ nhớt động học $\nu = 1.004 \times 10^{-6}\text{ m}^2/\text{s}$, vận hành trong vùng số vòng quay đặc trưng $n_s = 90 \div 150$.

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ệ hình thực chứng (positivism) kết hợp chặt chẽ giữa tam giác phương pháp: Mô hình hóa giải tích 1D $\rightarrow$ Mô phỏng số động lực học dòng chảy 2D/3D CFD $\rightarrow$ Khảo nghiệm thực nghiệm thủy lực vòng kín (closed-loop test rig). Thiết kế nghiên cứu đa tầng cho phép kiểm chứng chéo độc lập từng biến số hình học trước khi tích hợp vào tổ máy hoàn chỉnh.

Quy trình nghiên cứu rigorous

Quy trình nghiên cứu trải qua các giai đoạn nghiêm ngặt:

  • Giai đoạn 1 - Thiết kế mặt cắt kinh tuyến và khảo sát đơn biến: Thiết lập kích thước cơ bản ($D_1, D_2, b_1, b_2$) từ các thông số năng lượng đầu vào; phân tích độ nhạy của hiệu suất $\eta$ theo các biến số hình học đơn lẻ ($D_2, Z, \beta_{1B}, \beta_{2B}, \delta$).
  • Giai đoạn 2 - Thiết kế và sàng lọc 2D biên dạng cánh: Xây dựng 3 phương án biên dạng cánh (Mẫu 1: cánh cong truyền thống; Mẫu 2: cánh biến thiên độ cong bậc cao; Mẫu 3: cánh cung cong đối xứng). Tiến hành mô phỏng 2D khảo sát phân bố trường áp suất và động năng rối ($k$).
  • Giai đoạn 3 - Mô phỏng số 3D Navier-Stokes toàn tổ máy: Thiết lập bài toán mô phỏng 3D ở chế độ dừng và quay (MRF) cho toàn bộ hệ thống gồm buồng xoắn, vành cột trụ, cánh hướng, bánh công tác và ống hút. Đánh giá chi tiết nguy cơ xâm thực thông qua trường phân bố áp suất tĩnh $P_s$ và áp suất bốc hơi bão hòa $P_v$.
  • Giai đoạn 4 - Khảo nghiệm thực nghiệm chế độ kép: Chế tạo mẫu thử thực tế tại Công ty CP Chế tạo Bơm Hải Dương. Thực hiện quy trình đo đạc theo tiêu chuẩn quốc tế cho cả chế độ bơm và tuabin.

Data và phân tích

Phân tích số 3D sử dụng hệ phương trình Reynolds-Averaged Navier-Stokes (RANS) đóng kín bằng mô hình rối hai phương trình $k-\epsilon$ tiêu chuẩn và $k-\omega\text{ SST}$ trên nền tảng phần mềm chuyên dụng. Lưới tính toán cấu trúc đa khối (structured multi-block mesh) và lưới phi cấu trúc lăng trụ tại lớp biên với giá trị $y^+ < 5$ đảm bảo phân giải chính xác lớp biên dính nhớt.

Hệ thống thực nghiệm tại Hải Dương bao gồm hệ thống đường ống vòng kín đường kính DN200/DN300, trang bị cảm biến áp suất chênh áp dải đo $0 \div 1.6\text{ MPa}$ (cấp chính xác $\pm 0.1%$), lưu lượng kế điện từ (độ chính xác $\pm 0.2%$), thiết bị đo mô-men xoắn và tốc độ quay động cơ/máy phát (sai số $\pm 0.05%$). Phân tích sai số tổng hợp hệ thống tuân thủ nghiêm ngặt tiêu chuẩn ISO 9906, xác định sai số ngẫu nhiên dưới $\pm 0.85%$ và sai số toàn phần của phép đo hiệu suất không vượt quá $\pm 1.25%$.

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

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

  1. Quy luật chi phối của đường kính $D_2$ và góc đặt cánh $\beta_{2B}$: Khi tăng đường kính $D_2$, cột áp bơm $H_P$ và công suất tuabin $P_T$ tăng tuyến tính, song tổn thất ma sát đĩa $h_{disk}$ tăng theo lũy thừa bậc 5 của bán kính. Góc đặt cánh $\beta_{2B}$ tối ưu cho PaT $n_s = 104$ nằm trong dải hẹp $26^\circ \div 32^\circ$, giúp dung hòa giữa yêu cầu tăng cột áp bơm và hạn chế tách dòng ở lối vào chế độ tuabin.
  2. Tác động phi tuyến của số cánh $Z$: Kết quả phân tích xác định số cánh tối ưu cho máy PaT $n_s$ thấp là $Z = 7 \div 8$ cánh. Khi $Z < 6$, hệ số trượt $\gamma$ giảm mạnh làm suy giảm cột nước lý thuyết; khi $Z > 9$, tổn thất chèn dòng $\tau$ và tổn thất ma sát máng cánh $h_{fr}$ tăng vọt $14.2%$, triệt tiêu hoàn toàn lợi thế dẫn dòng.
  3. Ưu thế vượt trội của mẫu biên dạng cánh cung cong đối xứng: So với mẫu cánh cong thông thường, mẫu cánh cung cong đối xứng làm giảm gradient áp suất đột ngột tại mép vào ở chế độ bơm, loại bỏ vùng xoáy quẩn tách dòng tại mặt hút, qua đó tăng hiệu suất trung bình toàn hệ thống thêm $2.15%$ và mở rộng dải vận hành không xâm thực thêm $18%$.
  4. Ảnh hưởng điều tiết của độ mở cánh hướng ($a_0$): Khảo nghiệm thực tế với 3 độ mở $a_0 = 5\text{ mm}, 10\text{ mm}, 15\text{ mm}$ chỉ ra rằng việc tăng độ mở cánh hướng dịch chuyển điểm BEP về phía vùng lưu lượng lớn, đồng thời thay đổi tỷ số cột nước $h = H_T/H_P$ từ $1.22$ xuống $1.08$, cung cấp công cụ linh hoạt để đồng bộ hóa vùng làm việc tối ưu giữa hai chế độ.
  5. Độ tương hợp cao giữa mô hình lý thuyết, CFD và thực nghiệm: Sai số tương đối về cột áp và hiệu suất giữa tính toán giải tích 1D, mô phỏng 3D CFD và kết quả đo thực nghiệm tại điểm BEP luôn duy trì dưới $3.8%$ ở chế độ bơm và dưới $4.5%$ ở chế độ tuabin.

Implications đa chiều

  • Ý nghĩa lý thuyết: Xác lập mô hình toán học giải tích tường minh đầu tiên tại Việt Nam cho máy PaT $n_s$ thấp, đặt nền móng khoa học vững chắc thay thế các phương pháp nội suy kinh nghiệm thiếu tin cậy.
  • Đổi mới phương pháp luận: Chuẩn hóa quy trình thiết kế tích hợp 4 bước (Giải tích $\rightarrow$ Tối ưu hình học $\rightarrow$ CFD 3D $\rightarrow$ Thực nghiệm kiểm chứng) có khả năng chuyển giao cho các dòng máy thủy lực có $n_s$ khác nhau.
  • Ứng dụng thực tiễn công nghiệp: Cung cấp bộ thông số thiết kế hoàn chỉnh cho các nhà máy chế tạo cơ khí thủy khí trong nước (tiêu biểu như Công ty CP Chế tạo Bơm Hải Dương), nâng cao tỷ lệ nội địa hóa thiết bị cơ điện thủy điện tích năng.
  • Tác động chính sách: Cung cấp cơ sở khoa học định lượng phục vụ Quy hoạch điện VIII và các chiến lược phát triển nguồn điện phủ đỉnh quốc gia, giúp tối ưu hóa chi phí đầu tư các dự án NMTĐTN hàng nghìn tỷ đồng.

Limitations và Future Research

Nghiên cứu ghi nhận một số giới hạn học thuật cần tiếp tục mở rộng:

  1. Giới hạn đối tượng khảo sát tập trung chủ yếu vào dòng máy PaT ly tâm một cấp có số vòng quay đặc trưng thấp ($n_s = 90 \div 150$), chưa bao quát toàn diện các dòng máy hỗn lưu ($n_s > 200$) hoặc đa cấp.
  2. Mô hình dòng chảy 3D CFD ở chế độ dừng chưa mô tả đầy đủ các hiện tượng dao động áp suất không dừng (unsteady pressure pulsations) do tương tác thủy lực giữa bánh công tác quay và vành cánh hướng tĩnh (RSI - Rotor-Stator Interaction).
  3. Khảo nghiệm thực nghiệm mới chỉ thực hiện trên mô hình thử nghiệm tại phòng thí nghiệm công nghiệp, chưa tiến hành đo đạc dài hạn trong điều kiện vận hành thực địa với bùn cát và biến động thủy lực lưới điện.

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

  • Phát triển thuật toán tối ưu hóa đa mục tiêu tự động kết hợp trí tuệ nhân tạo (AI/Machine Learning) và thuật toán di truyền (NSGA-II) trên toàn bộ không gian tham số 3D của lá cánh.
  • Ứng dụng mô phỏng dòng không dừng LES/DES (Large Eddy Simulation) để khảo sát hiện tượng đập áp lực và vi bọt khí thực trong điều kiện chuyển tiếp chế độ vận hành (transient mode).
  • Mở rộng nghiên cứu sang dòng máy PaT cánh chéo (Deriaz) và các giải pháp điều tốc biến tần nhằm nâng cao hiệu suất vận hành non tải.

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

Luận án tạo ra dấu ấn học thuật tiên phong tại Việt Nam với tư cách là công trình nghiên cứu tiến sĩ toàn diện đầu tiên về thiết kế máy thủy lực thuận nghịch cho thủy điện tích năng. Các kết quả của luận án mở đường cho việc đăng tải các công trình trên hệ thống tạp chí cơ học thủy khí chuyên ngành quốc tế (ISI/Scopus), ước tính thu hút sự quan tâm trích dẫn cao từ cộng đồng nghiên cứu PaT toàn cầu.

Về mặt kinh tế - xã hội, việc làm chủ công nghệ thiết kế và chế tạo nội địa máy PaT giúp cắt giảm $30% - 45%$ chi phí nhập khẩu thiết bị từ các tập đoàn nước ngoài (như Alstom, Voith, Andritz hay Toshiba), tiết kiệm hàng triệu USD cho mỗi dự án NMTĐTN. Đồng thời, công trình thúc đẩy ngành cơ khí chế tạo máy thủy lực trong nước tiếp cận trình độ công nghệ tiên tiến, hỗ trợ đắc lực cho mục tiêu phát thải ròng bằng không (Net Zero) của Việt Nam thông qua tích hợp hiệu quả các nguồn năng lượng tái tạo gió và mặt trời vào lưới điện quốc gia.

Đố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 khung lý thuyết giải tích tổn thất 5 phân vùng hoàn chỉnh và quy trình mô phỏng CFD chuẩn mực để tiếp tục phát triển các đề tài chuyên sâu.
  • Các nhà khoa học và giảng viên chuyên ngành: Sở hữu tài liệu tham khảo có giá trị học thuật cao phục vụ công tác giảng dạy máy thủy lực và nghiên cứu năng lượng tái tạo tại các trường đại học kỹ thuật.
  • Kỹ sư R&D tại các doanh nghiệp chế tạo bơm và tuabin: Ứng dụng trực tiếp bản vẽ hình học, quy tắc đặt góc cánh và quy trình thử nghiệm để sản xuất các tổ máy PaT thương mại đáp ứng yêu cầu kỹ thuật khắt khe.
  • Các nhà hoạch định chính sách và chủ đầu tư năng lượng (EVN, Bộ Công Thương): Có căn cứ khoa học chính xác để thẩm định giải pháp kỹ thuật, tính toán tổng mức đầu tư và xây dựng tiêu chuẩn kỹ thuật quốc gia cho thiết bị thủy điện tích năng.

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ì?
Đóng góp độc đáo nhất là việc xây dựng thành công mô hình giải tích 1D tích hợp hệ số trượt Stodola - Gülich và hệ số chèn dòng $\tau$ với hệ thống phương trình tính toán 10 thành phần tổn thất năng lượng thủy lực trên 5 phân vùng cấu trúc cho dòng máy PaT $n_s$ thấp ở cả hai chiều dòng chảy thuận nghịch, giải quyết hạn chế của các mô hình kinh nghiệm đơn biến quốc tế.

2. Đổi mới phương pháp luận so với các công trình quốc tế trước đây?
Khác biệt với Derakhshan et al. (2008) hay Xuhe et al. (2014) vốn chỉ sử dụng CFD thuần túy để tối ưu hình học ngẫu nhiên, luận án kết hợp chặt chẽ giữa giải tích thủy lực cơ bản để định hướng dải thông số ban đầu, sàng lọc 2D/3D CFD và kiểm chứng thực nghiệm trực tiếp trên hệ thống đo đạc vòng kín tiêu chuẩn cao tại Hải Dương.

3. Phát hiện bất ngờ nhất từ dữ liệu thực nghiệm là gì?
Phát hiện bất ngờ nhất là biên dạng cánh dạng cung cong đối xứng đơn giản lại mang lại hiệu quả thủy lực và chống xâm thực vượt trội hơn hẳn các mẫu cánh biến thiên độ cong phức tạp, do loại bỏ được sự đảo chiều áp suất đột ngột khi dòng chảy đổi hướng giữa chế độ bơm và tuabin.

4. Giao thức tái lập nghiên cứu (Replication Protocol) có được cung cấp đầy đủ không?
Luận án cung cấp đầy đủ và chi tiết toàn bộ hệ thống công thức giải tích, kích thước hình học mặt cắt kinh tuyến, thông số lưới và điều kiện biên mô phỏng CFD, cùng bảng thông số thiết bị đo và sai số khảo nghiệm thực nghiệm, cho phép cộng đồng khoa học hoàn toàn có thể tái lập và mở rộng nghiên cứu.

5. Lộ trình nghiên cứu 10 năm tiếp theo được định hình như thế nào?
Lộ trình 10 năm hướng tới việc mở rộng dải $n_s$ cho dòng PaT hỗn lưu, ứng dụng thuật toán tối ưu hóa thông minh AI/CFD, nghiên cứu động lực học không dừng chống rung chấn RSI, và chế tạo thử nghiệm tổ máy PaT công suất lớn (vài chục MW) lắp đặt thực tế tại các công trình tích năng phân tán của Việt Nam.

Kết luận

  1. Luận án đã giải quyết trọn vẹn bài toán thiết kế máy thủy lực thuận nghịch bơm – tuabin có số vòng quay đặc trưng thấp ($n_s = 90 \div 150$), phục vụ trực tiếp chiến lược phát triển thủy điện tích năng quốc gia giai đoạn 2011 - 2030.
  2. Xây dựng thành công cơ sở lý thuyết giải tích dự báo đặc tính năng lượng $H(Q), P(Q), \eta(Q)$ trên cơ sở mô hình hóa tổn thất 5 phân vùng và tương tác trượt dòng - chèn dòng chính xác.
  3. Đề xuất và kiểm chứng thành công mẫu biên dạng cánh mới dạng một cung cong đối xứng, nâng cao hiệu suất tổ máy thêm $2.15%$ và cải thiện vượt bậc chất lượng xâm thực so với cánh truyền thống.
  4. Xác lập bộ thông số hình học tối ưu ($Z = 7 \div 8$, $\beta_{2B} = 26^\circ \div 32^\circ$, $\delta = 2^\circ \div 3^\circ$) cho dòng máy PaT $n_s = 104$, thiết lập tương quan tin cậy giữa độ mở cánh hướng và điểm BEP.
  5. Mở ra 3 hướng nghiên cứu khoa học mới: tối ưu hóa đa mục tiêu thông minh cho máy thủy lực thuận nghịch, phân tích động lực học tương tác Rotor - Stator không dừng, và công nghệ chế tạo nội địa tổ máy tích năng công suất lớn.
  6. Đặt nền tảng vững chắc đưa công nghệ cơ khí thủy khí Việt Nam hội nhập trình độ khu vực và quốc tế, tạo tiền đề khoa học quan trọng cho việc tự chủ thiết kế và chế tạo các thiết bị thủy điện tích năng trong tương lai.