Tổng quan về luận án

Trong bức tranh chuyển dịch năng lượng và tối ưu hóa hệ thống điện quốc gia theo Quyết định số 428/QĐ-TTg ngày 18/03/2016 của Thủ tướng Chính phủ về việc 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, việc xây dựng các nhà máy thủy điện tích năng (TĐTN) như Bác Ái (công suất $4 \times 300\text{ MW}$) và Đông Phù Yên (công suất $3 \times 300\text{ MW}$) giữ vai trò chiến lược sống còn nhằm phủ đỉnh biểu đồ phụ tải và duy trì tính ổn định của lưới điện. Thiết bị trung tâm của các công trình này là máy thủy lực thuận nghịch bơm – tuabin (Pump-as-Turbine - PaT). Tuy nhiên, rào cản lớn nhất đối với ngành cơ khí động lực trong nước là phần lớn các tổ máy PaT cột nước cao, số vòng quay đặc trưng thấp ($n_s = 90 - 150$) đều phải nhập khẩu toàn bộ từ nước ngoài với chi phí đắt đỏ do tính phức tạp trong thủy động lực học hai chiều.

Khoảng trống nghiên cứu (research gap) cốt lõi xuất phát từ thực tế: dòng chảy trong chế độ bơm là dòng khuếch tán tăng áp (decelerating flow) với nguy cơ cao về phân kỳ dòng, tách dòng và xâm thực, trong khi dòng chảy ở chế độ tuabin là dòng gia tốc hội tụ (accelerating flow). Việc tích hợp hai chế độ có cơ chế thủy lực đối nghịch vào một bánh công tác duy nhất đặt ra thách thức cực lớn trong việc lựa chọn các thông số hình học và biên dạng cánh. Các nghiên cứu kinh điển trước đây của Stepanoff (1957), Childs (1962), Sharma (1985), hay Gülich (2010) chủ yếu tập trung vào việc chuyển đổi các bơm ly tâm thông thường sang chạy ở chế độ tuabin mà chưa có một quy trình giải tích hoàn chỉnh kết hợp mô phỏng số 3D và thực nghiệm để thiết kế bánh công tác thuận nghịch tối ưu ngay từ đầu cho dải tỷ tốc thấp ($n_s \approx 104$).

Luận án tiến sĩ kỹ thuật cơ khí động lực của nghiên cứu sinh Nguyễn Thị Nhớ (2019) tại Trường Đại học Bách Khoa Hà Nội, dưới sự hướng dẫn của PGS. Trương Việt Anh và TS. Vũ Văn Trường, đã giải quyết trực diện bài toán này. Đề tài tập trung vào 3 câu hỏi nghiên cứu (RQ) và 2 giả thuyết khoa học (H):

  • RQ1: Làm thế nào để thiết lập mô hình toán học giải tích dự báo chính xác các đường đặc tính năng lượng ($H, P, \eta$) cho cả hai chế độ làm việc từ cấu trúc hình học bánh công tác?
  • RQ2: Mức độ ảnh hưởng định lượng của các thông số hình học then chốt ($D_2, \beta_{1B}, \beta_{2B}, Z$) đến sự thay đổi trường vận tốc, tổn thất thủy lực và hiệu suất của máy PaT tỷ tốc thấp là gì?
  • RQ3: Mẫu biên dạng cánh nào có khả năng dung hòa tối ưu phân bố áp suất và trường dòng chảy giữa chế độ bơm và tuabin?
  • H1: Tồn tại một dải thông số hình học biên dạng cánh tối ưu dạng một cung cong đối xứng giúp cực tiểu hóa tổn thất xoáy và va đập ở cả hai chiều quay.
  • H2: Việc tích hợp mô hình tính toán tổn thất phân đoạn lý thuyết với mô phỏng động lực học dòng chảy CFD 3D cho phép dự báo sai số dưới 10% so với thực nghiệm vận hành thực tế.

Khung lý thuyết của nghiên cứu được xây dựng trên phương trình cơ bản của Euler cho máy thủy lực cánh dẫn, lý thuyết trượt vận tốc dòng chảy (slip factor), kết hợp hệ thống giải tích phân rã tổn thất thủy lực nội tại ($h_{loss}$). Phạm vi nghiên cứu giới hạn trên tổ máy thí nghiệm PaT có $n_s = 104$, thử nghiệm vật lý trên hệ thống bệ thử hiện đại tại Công ty Cổ phần Chế tạo Bơm Hải Dương, tạo nền tảng đột phá cho công nghệ tự chủ chế tạo máy thủy điện tích năng tại Việt Nam.

Literature Review và Positioning

Lịch sử nghiên cứu máy thuận nghịch bắt đầu từ công bố nền tảng của Thoma và Kittredge (1931), khi phát hiện bơm ly tâm có thể vận hành ổn định như một tuabin phát điện. Trong nhiều thập kỷ sau đó, một dòng nghiên cứu lớn tập trung vào việc xây dựng các công thức bán thực nghiệm nhằm dự báo điểm hiệu suất cao nhất (Best Efficiency Point - BEP) của chế độ tuabin từ đặc tính bơm. Các học giả như Stepanoff (1957) đề xuất tỷ số $h = H_T/H_P = 1/\eta_P$ và $q = Q_T/Q_P = 1/\sqrt{\eta_P}$; Childs (196Childs, 1962) đưa ra quan hệ đồng nhất $h = 1/\eta_P$ và $q = 1/\eta_P$; trong khi Sharma (1985) đề xuất mô hình phi tuyến tinh chỉnh $h = 1/\eta_P^{1.2}$ và $q = 1/\eta_P^{0.8}$. Williams (1994) khi tiến hành thử nghiệm trên 35 mẫu bơm khác nhau theo tiêu chuẩn sai số elip BS516 đã chỉ ra rằng lý thuyết của Sharma tuy có độ tin cậy cao nhất nhưng vẫn có tới 20% số kết quả đo nằm ngoài vùng dự báo chuẩn.

Một nhánh nghiên cứu khác tập trung vào cơ cấu thủy động lực học cục bộ và thiết kế tối ưu hóa hình học. Gülich (2010) và Chapallaz (1992) khẳng định không tồn tại một quy luật phổ quát cho mọi loại máy PaT mà đặc tính thủy lực phụ thuộc chặt chẽ vào số vòng quay đặc trưng $n_s$ và cấu trúc tổn thất nội dòng. Derakhshan và cộng sự (2008) sử dụng phần mềm NUMECA FINETURBO giải phương trình 3D Navier-Stokes và AutoGrid 5 để thiết kế lại biên dạng cánh bơm ($n_s = 23.5$), giúp tăng hiệu suất chế độ tuabin thêm từ 3% đến 5%. Xuhe và cộng sự (2014) kết hợp thuật toán tối ưu đa mục tiêu Pareto trên 1.000 phương án biên dạng cánh PaT cho công trình có cột nước $298\text{ m}$, công suất $270\text{ MW}$, đưa hiệu suất tuabin đạt 91% và bơm đạt 90.8%.

Tranh luận học thuật quốc tế tồn tại rõ rệt giữa hai quan điểm: (1) Nhóm quan điểm truyền thống (Gülich, Stepanoff) ưu tiên tối ưu hóa bánh công tác theo chế độ bơm vì chế độ bơm có cơ chế dòng chảy nhạy cảm và dễ phát sinh bất ổn định thủy lực hơn; (2) Nhóm quan điểm hiện đại (Xuhe, Miao Sen-chun, Yang Sun-Sheng) đề xuất tối ưu hóa đồng thời đa mục tiêu dựa trên bề mặt Pareto nhằm tìm kiếm điểm thỏa hiệp tối ưu giữa hai chế độ.

So sánh với các công trình trên thế giới như nhà máy TĐTN Kazunogawa (Nhật Bản, cột áp $728\text{ m}$, công suất $412\text{ MW}$, $n = 500\text{ vòng/phút}$) hay nhà máy Helms (Hoa Kỳ, cột áp $531.6\text{ m}$, công suất $414\text{ MW}$), nghiên cứu của Nguyễn Thị Nhớ định vị chính xác khoảng trống công nghệ tại Việt Nam: phát triển một phương pháp thiết kế gốc cho máy PaT đơn cấp $n_s = 104$, giải quyết trọn vẹn sự mất cân xứng hình học giữa bánh công tác bơm (thường có $Z = 6 - 9$ cánh) và tuabin tâm trục Francis ($Z = 13 - 21$ cánh, đường kính mép vào tuabin nhỏ hơn bơm $D_{2T} \approx 0.83 D_{2P}$ theo Stelzer & Walters).

Đó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 khung lý thuyết truyền năng lượng cánh dẫn của Euler bằng việc thiết lập mô hình tính toán giải tích phân rã tổn thất thủy lực toàn diện $h_{loss}$ cho cả hai chiều dòng chảy: $$h_{loss} = h_{in} + h_{fr} + h_{dif} + h_{cir} + h_{disk} + h_{leg} + h_{cas} + h_{van} + h_{dr}$$

Trong đó, các thành phần tổn thất va đập mép cánh ($h_{in}$), ma sát bề mặt cánh ($h_{fr}$), tổn thất khuếch tán do độ loe máng dẫn ($h_{dif}$), xoáy quẩn cục bộ ($h_{cir}$), ma sát đĩa ($h_{disk}$), rò rỉ qua khe rãnh ($h_{leg}$), và tổn thất trong các bộ phận dẫn dòng (buồng xoắn $h_{cas}$, cánh hướng $h_{van}$, ống hút xả $h_{dr}$) được mô hình hóa bằng các hàm toán học tường minh gắn trực tiếp với góc đặt cánh $\beta_B$, góc tới $\Delta \beta$, độ mở cánh hướng $a_0$, và số cánh $Z$.

Đóng góp lý thuyết mang tính chuyển dịch tư duy (paradigm shift) của luận án là việc chứng minh sự khả thi của mẫu biên dạng cánh dạng "một cung cong đối xứng". Thay vì sử dụng các biên dạng cánh phức tạp uốn nhiều chiều vốn gây tổn thất phân ly dòng chảy khi đổi chiều quay, biên dạng cánh một cung cong đối xứng tạo ra sự cân bằng động lực học chất lưu, hạn chế sự hình thành các túi xoáy cục bộ và dải áp suất âm gây xâm thực tại mép cánh ở cả chế độ thuận và nghịch.

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: (1) Lý thuyết cơ học chất lỏng ứng dụng Navier-Stokes 3D; (2) Lý thuyết nhiệt động lực học máy thủy lực cánh dẫn; và (3) Lý thuyết tối ưu hóa hình học đa thông số. Luận án đã lượng hóa và thiết lập điều kiện biên hình học nghiêm ngặt cho dòng máy PaT $n_s$ thấp ($n_s \approx 104$):

  1. Tỷ số đường kính bánh công tác: $d_2^* = D_2/D_1 = 1.74 - 1.91$
  2. Góc đặt cánh tại mép trong: $\beta_{1B} = 24^\circ - 29^\circ$
  3. Góc đặt cánh tại mép ngoài: $\beta_{2B} = 30^\circ - 35^\circ$
  4. Số lượng cánh dẫn bánh công tác: $Z = 9 - 13$ cánh

Sự kết hợp này thiết lập biên giới hạn vận hành (boundary conditions) chặt chẽ, loại bỏ vùng làm việc bất ổn định dạng chữ "S" (S-shaped characteristic curve) thường gặp tại vùng khởi động hoặc quá độ của tuabin Francis thuận nghịch.

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ý thực chứng (positivism) với thiết kế nghiên cứu đa tầng kết hợp tam giác đạc phương pháp (methodological triangulation): (1) Mô hình hóa giải tích toán học; (2) Mô phỏng số động lực học dòng chảy 3D CFD (Computational Fluid Dynamics); và (3) Kiểm chứng thực nghiệm vật lý trên mô hình thực tế chế tạo tại nhà máy. Thiết kế đa tầng này đảm bảo tính nhất quán từ các giả thiết vi mô của hạt lưu chất đến các thông số năng lượng vĩ mô của toàn tổ máy.

Quy trình nghiên cứu rigorous

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

  • Thiết kế hình học và mô phỏng 2D/3D: Sử dụng bộ phần mềm chuyên dụng ANSYS TurboGrid để tạo lưới cấu trúc phân giải cao tại các vùng biên máng dẫn và mép cánh, kết hợp giải hệ phương trình Reynolds-Averaged Navier-Stokes (RANS) 3D. Khảo sát chi tiết phân bố áp suất tĩnh $P_s$, năng lượng động rối (turbulent kinetic energy), và các vector vận tốc tương đối $w$, vận tốc kinh tuyến $c_m$ tại các mặt cắt 10%, 50%, và 90% bề rộng máng dẫn cánh.
  • Quy trình thực nghiệm: Chế tạo mô hình thực tế với bánh công tác hợp kim đồng bộ theo thông số tối ưu tại Công ty Cổ phần Chế tạo Bơm Hải Dương. Hệ thống thử nghiệm khép kín gồm đầy đủ buồng xoắn, vành cánh hướng điều chỉnh được ($a_0 = 5\text{ mm}, 10\text{ mm}, 15\text{ mm}$), cột trụ, ống hút/xả, động cơ - máy phát điều tốc biến tần, van điện từ và hệ thống cảm biến đo đạc tự động.

Data và phân tích

Độ tin cậy và giá trị đo lường được kiểm soát theo Tiêu chuẩn Quốc gia TCVN 8639:2011 về khảo nghiệm máy thủy lực cánh dẫn. Phép phân tích sai số lan truyền toàn phần (propagation of uncertainty) được thực hiện chi tiết cho mọi đại lượng đo:

  • Sai số giới hạn đo áp suất: $\delta p = 0.27%$
  • Sai số giới hạn đo lưu lượng: $\delta Q = 0.40%$
  • Sai số giới hạn đo mô men xoắn trục: $\delta M = 0.18%$
  • Sai số giới hạn đo số vòng quay: $\delta n = 0.07%$

Tính toán sai số tổng hợp hệ thống: $$\delta H = \sqrt{0.77^2 + 0.47^2} = 0.90% < 1.0% \text{ (chuẩn TCVN 8639)}$$ $$\delta Q = \sqrt{0.40^2 + 1.40^2} = 1.45% < 1.5% \text{ (chuẩn TCVN 8639)}$$ $$\delta P = \sqrt{0.18^2 + 0.20^2} = 0.27% < 0.5% \text{ (chuẩn TCVN 8639)}$$ $$\delta \eta = \sqrt{0.90^2 + 1.45^2 + 0.27^2} = 1.72% < 4.0% \text{ (chuẩn TCVN 8639)}$$

Tất cả các sai số đo đạc thực nghiệm đều nằm sâu trong giới hạn cho phép của tiêu chuẩn quốc gia, khẳng định tính chuẩn xác và khả năng tái lập hoàn toàn của tập dữ liệu thực nghiệm.

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

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

  1. Hiệu suất năng lượng và tính bất đối xứng thủy lực: Tại điểm làm việc tốt nhất (BEP), chế độ bơm đạt hiệu suất đỉnh $\eta_{P,\max} = 77.07%$ tại lưu lượng $Q = 242.5\text{ m}^3/\text{h}$ và cột áp $H = 8.34\text{ m}$. Trong khi đó, chế độ tuabin đạt hiệu suất đỉnh $\eta_{T,\max} = 72.80%$ tại lưu lượng $Q = 315.0\text{ m}^3/\text{h}$ và cột áp $H = 9.49\text{ m}$. Kết quả chứng minh rõ nét: khi vận hành ở chế độ tuabin, máy đòi hỏi lưu lượng và cột áp cao hơn so với chế độ bơm, nhưng hiệu suất lại thấp hơn từ $4.27%$ đến $6.39%$.
  2. Quy luật tỷ số BEP tương thích dữ liệu quốc tế: Tỷ số cột áp thực nghiệm đạt $h = H_T/H_P = 1.07 - 1.30$ và tỷ số lưu lượng đạt $q = Q_T/Q_P = 1.11 - 1.16$. Kết quả này nằm hoàn hảo trong vùng dự báo của Sharma (1985) và Alatorre-Frenk (1994), khẳng định tính đúng đắn của phương pháp thiết kế gốc.
  3. Bản đồ phân bố tổn thất thủy lực nội dòng: Dữ liệu mô phỏng 3D chỉ ra tổn thất trong vùng bánh công tác luôn chiếm tỷ trọng áp đảo ($>50%$) tổng tổn thất toàn máy. Cụ thể:
    • Chế độ bơm: Tổn thất bánh công tác chiếm $56.20%$, vùng cánh hướng van điều tiết chiếm $18.60%$, vùng ống hút chiếm $12.90%$.
    • Chế độ tuabin: Tổn thất bánh công tác chiếm $59.13%$, vùng buồng xoắn chiếm $11.69%$, vùng cánh hướng chỉ chiếm $7.92%$. Tổn thất rò rỉ qua khe rãnh ở hai chế độ là tương đương nhau và chiếm tỷ lệ nhỏ nhất.
  4. Tác động phi tuyến của độ mở cánh hướng ($a_0$): Khi tăng độ mở cánh hướng từ $a_{01} = 5\text{ mm}$ lên $a_{02} = 10\text{ mm}$ và $a_{03} = 15\text{ mm}$, điểm BEP có xu hướng dịch chuyển sang phải:
    • Chế độ bơm: $Q_{BEP}$ tăng từ $218.22\text{ m}^3/\text{h}$ lên $300.0\text{ m}^3/\text{h}$, trong khi cột áp $H$ giảm từ $8.82\text{ m}$ xuống $7.71\text{ m}$.
    • Chế độ tuabin: $Q_{BEP}$ tăng từ $270.0\text{ m}^3/\text{h}$ lên $320.9\text{ m}^3/\text{h}$, trong khi cột áp $H$ giảm từ $9.78\text{ m}$ xuống $8.94\text{ m}$. Phương án $a_{02} = 10\text{ mm}$ cho hiệu suất đỉnh cao nhất, tuy nhiên phương án $a_{01} = 5\text{ mm}$ lại là phương án tối ưu thực tế vì bao trọn vùng vận hành an toàn theo cột áp tĩnh thực tế ($H_0 = 8 - 10\text{ m}$).
  5. Mức độ tương thích mô hình giải tích: Sai số giữa mô hình tính toán lý thuyết và thực nghiệm trong dải vận hành tiêu chuẩn ($90% - 110% Q_{design}$) được kiểm chứng luôn dưới $10%$.

Implications đa chiều

  • Về mặt học thuật: Cung cấp cơ sở lý thuyết định lượng về cơ chế tương tác thủy lực giữa bánh công tác và các bộ phận dẫn dòng hai chiều, lấp đầy khoảng trống dữ liệu thực nghiệm cho các dòng máy thủy lực PaT $n_s$ thấp tại Đông Nam Á.
  • Về mặt công nghệ và công nghiệp: Cung cấp quy trình thiết kế hoàn chỉnh giúp các doanh nghiệp cơ khí trong nước (như Công ty CP Chế tạo Bơm Hải Dương) hoàn toàn làm chủ quy trình tính toán, đúc, gia công CNC và kiểm thử các tổ máy thủy lực tích năng mà không cần phụ thuộc vào bản quyền sáng chế nước ngoài.
  • Về mặt chính sách: Cung cấp luận cứ khoa học thực chứng để Tập đoàn Điện lực Việt Nam (EVN) và Bộ Công Thương xây dựng tiêu chuẩn kỹ thuật cho các dự án thủy điện tích năng trong Quy hoạch điện quốc gia, tối ưu hóa nguồn vốn đầu tư công.

Limitations và Future Research

Luận luận án thẳng thắn thừa nhận 4 hạn chế khoa học nội tại:

  1. Nghiên cứu mới dừng lại ở việc thử nghiệm mô hình vật lý đơn cấp với $n_s = 104$ trên quy mô phòng thí nghiệm; chưa khảo sát các hiệu ứng tỷ lệ (scale effects) khi phóng lớn kích thước lên các tổ máy công nghiệp 300 MW.
  2. Chưa phân tích sâu về tương tác động rotor-stator (RSI - Rotor-Stator Interaction) và dao động áp suất tần số cao phát sinh do khe hở giữa cánh bánh công tác và vành cánh hướng.
  3. Nghiên cứu tập trung chủ yếu vào trạng thái xác lập (steady state), chưa mô phỏng toàn diện các chế độ quá độ động lực học (transient regimes) như quá trình lồng tốc, quá trình đóng mở van đột ngột, khởi động và ngắt tải.
  4. Chưa tích hợp thuật toán học máy (Machine Learning) để tự động hóa quá trình tối ưu hóa đa mục tiêu biên dạng cánh 3D không gian phức tạp.

Các hướng nghiên cứu tiếp theo bao gồm:

  • Mở rộng khảo sát ảnh hưởng của hình học buồng xoắn và góc nghiêng cột trụ đến trường dòng chảy xoáy.
  • Thiết lập mô hình mô phỏng chuyển tiếp pha xâm thực 3D đa pha (multiphase cavitation model) trong điều kiện quá độ.
  • Nghiên cứu cơ chế giảm thiểu rung động và mỏi vật liệu do tương tác thủy lực không dừng.
  • Phát triển phần mềm tự động hóa thiết kế tối ưu hóa bánh công tác PaT ứng dụng trí tuệ nhân tạo.

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

  • Tác động học thuật: Các kết quả của luận án đã được công bố trên 9 công trình khoa học uy tín, bao gồm 2 bài báo trên Tạp chí Cơ khí Việt Nam, 4 báo cáo tại Hội nghị Cơ học Thủy khí toàn quốc và Hội nghị Khoa học Thường niên Đại học Thủy Lợi, cùng 3 công bố quốc tế tại Seventh International Conference and Exhibition on Water Resources and Renewable Energy Development in Asia. Dự kiến công trình sẽ là tài liệu tham khảo nền tảng với tiềm năng trích dẫn cao cho các nghiên cứu sinh chuyên ngành máy thủy lực.
  • Tác động ngành công nghiệp: Trực tiếp chuyển giao công nghệ cho ngành chế tạo bơm và tuabin trong nước, nâng cao tỷ lệ nội địa hóa thiết bị cơ khí thủy công lên tới 40 - 60% cho các dự án TĐTN vừa và nhỏ.
  • Tác động xã hội và môi trường: Thúc đẩy hiệu quả tích trữ năng lượng tái tạo (điện gió, điện mặt trời), giảm phát thải khí nhà kính từ các nhà máy nhiệt điện đốt than khi vận hành phủ đỉnh.

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

  • Nghiên cứu sinh và Giảng viên ngành Cơ khí Thủy khí: Tiếp cận hệ thống phương pháp luận chặt chẽ, mô hình toán học giải tích tường minh và bộ dữ liệu thực nghiệm chuẩn hóa để phát triển các đề tài chuyên sâu.
  • Kỹ sư R&D tại các nhà máy cơ khí chế tạo: Sở hữu bộ thông số thiết kế hình học định lượng ($\beta_{1B}, \beta_{2B}, D_2/D_1, Z$) và quy trình kiểm thử sai số theo TCVN 8639:2011 để ứng dụng trực tiếp vào sản xuất.
  • Ban Quản lý Dự án Thủy điện và EVN: Có công cụ kỹ thuật độc lập để thẩm tra, đánh giá hiệu quả kinh tế - kỹ thuật của các gói thầu cung cấp thiết bị thủy lực thuận nghịch nhập khẩu.
  • Cơ quan Hoạch định Chính sách Năng lượng: Cơ sở thực chứng để xây dựng lộ trình nội địa hóa thiết bị năng lượng tái tạo theo Chiến lược Phát triển Khoa học và Công nghệ Quốc gia.

Câu hỏi chuyên sâu

  1. Đóng góp lý thuyết độc đáo nhất của luận án là gì và mở rộng lý thuyết nào? Luận án đã mở rộng phương trình năng lượng Euler và lý thuyết dòng tia trượt (slip factor theory) bằng việc xây dựng mô hình phân rã tổn thất thủy lực nội dòng 9 thành phần ($h_{loss}$) gắn trực tiếp với góc đặt cánh và thông số hình học cho cả hai chiều quay. Đóng góp đột phá là việc chứng minh tính ưu việt của mẫu cánh dạng một cung cong đối xứng giúp triệt tiêu các túi xoáy cục bộ và dung hòa tổn thất giữa chế độ bơm và tuabin.

  2. Cải tiến phương pháp luận so với các nghiên cứu quốc tế trước đây? So với nghiên cứu của Derakhshan (2008) chỉ dùng mô phỏng CFD cho một bơm cụ thể hay Williams (1994) chỉ thuần túy thống kê sai số thực nghiệm, luận án đã thiết lập quy trình khép kín ba bước: Mô hình toán giải tích phân đoạn $\rightarrow$ Mô phỏng số 3D RANS phân tích trường dòng $\rightarrow$ Thực nghiệm vật lý trên trạm thử chuẩn công nghiệp có sai số được kiểm soát ngặt nghèo theo TCVN 8639:2011 ($\delta \eta = 1.72%$).

  3. Phát hiện nào gây bất ngờ nhất từ dữ liệu thực nghiệm? Phương án độ mở cánh hướng $a_{02} = 10\text{ mm}$ mang lại hiệu suất đỉnh lớn nhất cho cả bơm ($77.07%$) và tuabin ($72.80%$), nhưng khi đối chiếu với đường đặc tính cột áp tĩnh thực tế ($H_0 = 8 - 10\text{ m}$), phương án này lại không thể vận hành được do đòi hỏi cột áp vượt quá khả năng cấp của trạm. Trái lại, phương án $a_{01} = 5\text{ mm}$ tuy có hiệu suất đỉnh thấp hơn đôi chút nhưng lại bao trọn vùng vận hành an toàn của toàn tổ máy.

  4. Luận án có cung cấp quy trình tái lập (replication protocol) hoàn chỉnh không? Có. Luận án công bố chi tiết toàn bộ các phương trình toán học giải tích, kích thước hình học chuẩn hóa của bánh công tác ($D_1, D_2, b_1, b_2, Z$), hệ số lưới ANSYS TurboGrid, sơ đồ bố trí cảm biến đo đạc, công thức lan truyền sai số và quy trình vận hành bệ thử nghiệm tại Công ty CP Chế tạo Bơm Hải Dương.

  5. Chương trình nghiên cứu 10 năm tiếp theo được định hình ra sao? Mở rộng từ mô hình PaT đơn cấp $n_s = 104$ sang dòng máy nhiều cấp cột nước siêu cao ($H > 500\text{ m}$); nghiên cứu động lực học xâm thực phi dừng trong giai đoạn quá độ; ứng dụng vật liệu composite chống mài mòn xâm thực; và tích hợp AI trong điều khiển thông minh góc mở cánh hướng theo biểu đồ giá điện thời gian thực.

Kết luận

  1. Thiết lập thành công mô hình toán học giải tích: Dự báo chính xác đặc tính năng lượng ($H, P, \eta$) của máy PaT hai chiều với sai số dưới $10%$ so với thực nghiệm trong dải lưu lượng làm việc danh định.
  2. Sáng tạo mẫu biên dạng cánh một cung cong đối xứng: Giải quyết xuất sắc xung đột thủy động lực học giữa chiều quay khuếch tán (bơm) và hội tụ (tuabin), hạn chế tối đa nguy cơ phân tách dòng và xâm thực.
  3. Xác lập bộ thông số hình học tối ưu cho máy PaT $n_s = 104$: $\beta_{1B} = 24^\circ - 29^\circ$, $\beta_{2B} = 30^\circ - 35^\circ$, $d_2^* = 1.74 - 1.91$, $Z = 9 - 13$ cánh, cung cấp cẩm nang thiết kế chuẩn mực cho các kỹ sư cơ khí thủy khí.
  4. Lượng hóa chi tiết bức tranh tổn thất thủy lực: Chỉ ra tổn thất trong bánh công tác chiếm trên $50%$ tổng tổn thất hệ thống ($56.20%$ ở chế độ bơm và $59.13%$ ở chế độ tuabin), xác định chính xác trọng tâm cần tối ưu hóa hình học.
  5. Chế tạo và thử nghiệm thành công mô hình PaT đầu tiên tại Việt Nam: Kiểm chứng thực nghiệm đạt hiệu suất $\eta_P = 77.07%$ và $\eta_T = 72.80%$, sai số đo đạc đạt chuẩn quốc gia TCVN 8639:2011.
  6. Mở ra 3 hướng nghiên cứu khoa học chuyên sâu mới: Động lực học dòng chảy quá độ trong máy thủy điện tích năng; Tối ưu hóa tương tác thủy lực đa bộ phận dẫn dòng (buồng xoắn - cánh van - bánh công tác); và Công nghệ chế tạo nội địa hóa tổ máy thủy điện tích năng quy mô công nghiệp phục vụ an ninh năng lượng quốc gia.