Tổng quan nghiên cứu

Việt Nam sở hữu nguồn tài nguyên năng lượng tái tạo dồi dào với lượng bức xạ mặt trời trung bình hàng năm dao động từ 4,3 đến 5,7 triệu kWh/m2, cùng số giờ nắng tại các khu vực Nam Trung Bộ và Tây Nguyên đạt từ 2.000 đến 2.600 giờ/năm. Theo Chiến lược phát triển năng lượng tái tạo của Việt Nam, tỷ lệ thâm nhập của nguồn năng lượng sạch này được định hướng đạt 0,5% vào năm 2020, nâng lên 6% vào năm 2030 và hướng tới mốc 20% vào năm 2050. Tuy nhiên, rào cản kỹ thuật lớn nhất hiện nay nằm ở chi phí công nghệ và độ tin cậy của các giải pháp khai thác quy mô công nghiệp.

Công nghệ tháp đối lưu sử dụng năng lượng mặt trời (Solar Chimney Power Plant) là giải pháp đột phá kết hợp hiệu ứng tích nhiệt nhà kính, đối lưu thẳng đứng và turbine phát điện. Thách thức cốt lõi của hệ thống nằm ở kết cấu tháp siêu cao chịu tác động đồng thời của dòng lưu chất nhiệt bên trong và áp lực khí động học bên ngoài. Luận văn tập trung nghiên cứu phương pháp số nhằm phân tích trường cặp đôi lưu chất - nhiệt - kết cấu (Fluid-Structure Interaction - FSI) trên phần mềm ANSYS.

Nghiên cứu lấy mô hình tháp Manzanares tại Tây Ban Nha làm chuẩn đối sánh với chiều cao thân tháp 194,6 m và bán kính mái thu nhiệt 122 m. Mục tiêu cụ thể là xác định trường phân bố vận tốc, biến thiên nhiệt độ và đánh giá vị trí tối ưu để bố trí turbine phát điện 50 kW. Kết quả mô phỏng số đạt độ chính xác cao với sai số vận tốc chỉ 2,26% so với dữ liệu thực nghiệm, tạo tiền đề khoa học vững chắc cho việc thiết kế và nội địa hóa các công trình điện mặt trời dạng tháp tại Việt Nam.

Cơ sở lý thuyết và phương pháp nghiên cứu

Khung lý thuyết áp dụng

Nghiên cứu được xây dựng trên nền tảng của thuyết nhiệt động lực học đối lưu tự nhiên, cơ học lưu chất tính toán (Computational Fluid Dynamics - CFD) và lý thuyết tương tác lưu chất - kết cấu. Hệ thống tháp vận hành theo nguyên lý hiệu ứng ống khói: không khí dưới mái kính hấp thụ bức xạ nhiệt mặt trời bị đốt nóng, giảm khối lượng riêng và chuyển động hướng tâm vào cổ tháp tạo thành luồng khí đối lưu liên tục làm quay turbine.

Mô hình toán học của dòng khí không nén được (số Mach dưới 0,3) được mô tả qua hệ phương trình vi phân đạo hàm riêng phi tuyến:

  1. Phương trình liên tục bảo toàn khối lượng.
  2. Hệ phương trình Navier-Stokes biểu diễn chuyển động của lưu chất thực có ma sát, kết hợp toán tử nhớt động học và gradient áp suất.
  3. Phương trình vi phân năng lượng và phương trình biến thiên động lượng tích phân trên thể tích kiểm soát.

Các khái niệm cơ bản được chuẩn hóa trong nghiên cứu bao gồm: mật độ lưu chất, hệ số nhớt động lực học, ứng suất ma sát lớp biên và suất đàn hồi thể tích của chất khí.

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

Nghiên cứu áp dụng quy trình tính toán số hiện đại thông qua sự kết hợp giữa Phương pháp Thể tích Hữu hạn (Finite Volume Method - FVM) trên mô-đun ANSYS CFX và Phương pháp Phần tử Hữu hạn (Finite Element Method - FEM). Dữ liệu đầu vào sử dụng các thông số kỹ thuật thực tế từ mô hình tháp Manzanares với chiều cao trung bình mái thu 1,85 m và công suất thiết kế 50 kW.

Để tối ưu hóa tài nguyên phần cứng máy tính, mô hình tính toán được thiết lập đối xứng với tỷ lệ 1/64 toàn bộ hình học tháp. Quy mô phân chia lưới (cỡ mẫu tính toán) đạt 26.802 nút (nodes) và 54.283 phần tử (elements) dạng tứ diện (tetrahedrons). Kỹ thuật chia lưới lớp biên Inflation được áp dụng sát bề mặt thành tháp để thu thập chính xác gradient vận tốc và ứng suất tiếp tuyến.

Điều kiện biên được thiết lập nghiêm ngặt: áp suất đầu vào 1 atm với nhiệt độ môi trường 23 °C (296 K); áp suất đầu ra tại đỉnh tháp 194,6 m đạt 0,9771 atm với nhiệt độ 21,6 °C; bề mặt hấp thụ nhận thông lượng nhiệt mặt trời 1.000 W/m2 với độ nhám 15 mm. Lý do lựa chọn FVM kết hợp FEM là khả năng giải quyết chính xác các tương tác phi tuyến phức tạp giữa dòng khí chuyển động nhiệt và biến dạng đàn hồi của vỏ tháp mà các công thức giải tích cổ điển không thể thực hiện được.

Kết quả nghiên cứu và thảo luận

Những phát hiện chính

Phân tích số liệu từ mô hình mô phỏng ANSYS CFX ghi nhận các kết quả định lượng nổi bật:

Thứ nhất, trường vận tốc dòng khí gia tăng đột biến khi chuyển tiếp từ khu vực thu nhiệt rộng lớn vào cổ tháp. Vận tốc dòng khí tại vị trí đặt turbine ở cao độ 9 m đạt giá trị 42,6199 m/s. Khi đối sánh với số liệu thực nghiệm đo đạc thực tế tại trạm Manzanares là 41,677 m/s, mức sai số tương đối chỉ đạt khoảng 2,26%, khẳng định độ tin cậy xuất sắc của mô hình số.

Thứ hai, trường nhiệt độ không khí bên dưới mái thu nhiệt tăng trưởng tịnh tiến từ mép ngoài vào trung tâm. Không khí được gia nhiệt từ 296 K (23 °C) lên 311,1 K (38,1 °C), tạo ra mức chênh lệch nhiệt độ $\Delta T = 15,1\text{ K}$. Kết quả này hoàn toàn tương thích với dải gia nhiệt tham chiếu từ 15 K đến 20 K của hệ thống thu nhiệt thực tế.

Thứ ba, sự thay đổi tiết diện đột ngột tại chân tháp tạo ra gradient áp suất âm lớn, đóng vai trò là động lực chính duy trì lưu lượng khí ổn định qua hệ thống 4 cánh turbine biên dạng FX W-151-A, bảo đảm hiệu suất phát điện liên tục.

Thảo luận kết quả

Cơ chế gia tăng vận tốc đột ngột tại cao độ 9 m được giải thích bởi định luật bảo toàn dòng liên tục kết hợp với lực đẩy nổi Archimedes. Khi dòng khí nóng bị nén từ diện tích thu nhiệt bán kính 122 m vào đường kính tháp chỉ 10,16 m (bán kính 5,08 m), động năng dòng khí tăng mạnh tương ứng với sự sụt giảm thế năng áp suất.

Dữ liệu mô phỏng có thể được trực quan hóa sinh động qua biểu đồ đường dòng vận tốc (Streamlines) theo trục chiều cao từ 0 m đến 194,6 m và bảng đối chiếu các giá trị nhiệt động lực học. Biểu đồ phân bố vector vận tốc trên mặt cắt 1/2 mô hình chứng minh dòng chảy đạt trạng thái ổn định, không xuất hiện các xoáy khí cục bộ làm tiêu hao năng lượng tại vùng chân tháp.

So với các công trình nghiên cứu quốc tế trước đây, việc áp dụng mô hình đối xứng 1/64 trong nghiên cứu này giúp giảm hơn 80% thời gian xử lý dữ liệu tính toán mà vẫn bảo toàn độ chính xác của các đường đặc tính khí động học. Đây là cơ sở khoa học quan trọng để mở rộng mô phỏng cho các siêu dự án tháp năng lượng cao tới 1.000 m.

Đề xuất và khuyến nghị

Nhằm hiện thực hóa công nghệ tháp năng lượng mặt trời tại các khu vực giàu bức xạ, nghiên cứu đưa ra 4 nhóm giải pháp kỹ thuật và chính sách cụ thể:

  1. Tối ưu hóa vị trí lắp đặt turbine phát điện: Đơn vị thiết kế cần bố trí cụm turbine gió tại dải cao độ từ 8,5 m đến 9,5 m tính từ mặt đất. Đây là vùng hội tụ dòng khí có vận tốc đạt đỉnh trên 42,6 m/s, giúp tăng hiệu suất thu hồi động năng thêm 5% đến 8% trước năm 2027. Đơn vị thực hiện: Kỹ sư thiết kế khí động học và Viện nghiên cứu năng lượng.
  2. Tích hợp hệ thống ống trữ nhiệt tuần hoàn: Doanh nghiệp dự án cần lắp đặt mạng lưới ống chứa nước kín dưới 100% diện tích mái thu nhiệt 122 m. Giải pháp này giúp duy trì độ chênh nhiệt $\Delta T \ge 12\text{ K}$ vào ban đêm, bảo đảm hệ thống vận hành phát điện ổn định 24/7 từ quý 3/2026. Đơn vị thực hiện: Đơn vị tổng thầu EPC và kỹ sư vận hành nhà máy.
  3. Chuẩn hóa quy trình mô phỏng FSI đa môi trường: Nhóm nghiên cứu kết cấu cần ban hành bộ quy chuẩn thiết lập lưới Inflation và phân tích cặp đôi FSI với quy mô tối thiểu 50.000 phần tử trên ANSYS cho các kết cấu tháp thép từ 200 m đến 1.000 m, hoàn thành trước tháng 12/2026. Đơn vị thực hiện: Nhóm tính toán R&D kết cấu công trình.
  4. Xây dựng khung cơ chế ưu đãi đầu tư điện mặt trời dạng tháp: Cơ quan quản lý nhà nước cần ban hành biểu giá điện hỗ trợ FIT chuyên biệt và chính sách vay vốn ưu đãi cho công nghệ tháp đối lưu, hướng tới mục tiêu hạ giá thành điện năng xuống dưới 7 US cents/kWh giai đoạn 2028 - 2030. Đơn vị thực hiện: Bộ Công Thương và Tập đoàn Điện lực Việt Nam.

Đối tượng nên tham khảo luận văn

Nội dung và kết quả nghiên cứu mang lại giá trị thực tiễn sâu sắc cho 4 nhóm đối tượng trọng tâm:

  1. Kỹ sư Cơ học lưu chất tính toán (CFD) và Kết cấu: Nắm vững quy trình thiết lập mô hình thể tích hữu hạn 26.802 nút, kỹ thuật xử lý điều kiện biên áp suất 0,9771 atm và phương pháp giải bài toán trường cặp đôi FSI trên phần mềm ANSYS CFX để áp dụng vào các kết cấu dạng ống siêu cao tầng.
  2. Nhà phát triển dự án năng lượng tái tạo: Tiếp cận dữ liệu kỹ thuật chi tiết của tháp 194,6 m công suất 50 kW để xây dựng báo cáo nghiên cứu tiền khả thi, tính toán diện tích mặt bằng mái thu và phân tích hiệu quả tài chính tại các tỉnh duyên hải miền Trung có bức xạ từ 4,3 triệu kWh/m2.
  3. Giảng viên, Học viên cao học và Sinh viên kỹ thuật: Sử dụng hệ thống lý thuyết về phương trình Navier-Stokes, phương trình Euler và phương pháp phần tử hữu hạn FVM/FEM làm tài liệu tham khảo giảng dạy, học tập và phát triển đề tài nghiên cứu chuyên sâu.
  4. Chuyên gia hoạch định chiến lược năng lượng: Tham khảo mô hình công nghệ phát điện không tiêu thụ nguồn nước và không phát thải khí nhà kính để xây dựng các kịch bản chuyển dịch năng lượng đạt mốc 20% điện mặt trời vào năm 2050.

Câu hỏi thường gặp

Hệ thống tháp năng lượng mặt trời có thể phát điện vào ban đêm không? Hệ thống hoàn toàn có khả năng phát điện 24/7 nhờ mạng lưới ống nước ngầm bố trí dưới mái thu nhiệt bán kính 122 m. Nước hấp thụ nhiệt lượng vào ban ngày và giải phóng nhiệt chậm vào ban đêm, duy trì độ chênh nhiệt từ 8 K đến 12 K để tạo luồng đối lưu làm quay turbine 50 kW liên tục.

Tại sao luận văn chỉ mô phỏng phân đoạn 1/64 của tháp Manzanares? Do tháp có kết cấu hình trụ tròn cao 194,6 m mang tính đối xứng hình học và thủy nhiệt hoàn hảo qua trục. Việc mô phỏng 1/64 giúp giảm số lượng phần tử tính toán xuống còn 54.283 phần tử, tiết kiệm hơn 80% thời gian xử lý máy tính mà vẫn giữ sai số dưới 2,3% so với thực nghiệm.

Vận tốc dòng khí 42,6199 m/s đạt được ở vị trí nào và có vai trò gì? Vận tốc đạt cực đại tại vị trí cao độ 9 m tính từ chân tháp, nơi tiết diện thu hẹp từ mái kính sang thân ống bán kính 5,08 m. Kết quả 42,6199 m/s khớp với lý thuyết 41,677 m/s, chứng minh đây là vị trí tối ưu nhất để đặt turbine gió 4 cánh thu hồi động năng.

Phương pháp Thể tích Hữu hạn (FVM) mang lại ưu thế gì trong bài toán này? FVM bảo toàn chính xác các đại lượng vật lý như khối lượng, động lượng và năng lượng trên từng thể tích kiểm soát rời rạc. Phương pháp này vượt trội khi xử lý phương trình Navier-Stokes cho dòng khí đối lưu tự nhiên, trong khi FEM được kết hợp để giải bài toán biến dạng ứng suất vỏ tháp.

Công nghệ tháp năng lượng mặt trời có phù hợp với khí hậu Việt Nam không? Việt Nam có bức xạ mặt trời cao từ 4,3 đến 5,7 triệu kWh/m2/năm và nhiều vùng đất khô cằn. Công nghệ này không sử dụng nước làm mát, hấp thụ được cả bức xạ trực tiếp lẫn khuếch tán trong ngày mây, rất thích hợp triển khai tại các tỉnh Nam Trung Bộ và Tây Nguyên.

Kết luận

  • Ứng dụng thành công phương pháp số trên phần mềm ANSYS CFX để mô phỏng chính xác trường cặp đôi lưu chất - nhiệt trong tháp năng lượng mặt trời.
  • Kiểm chứng vận tốc cực đại tại vị trí đặt turbine 9 m đạt 42,6199 m/s với sai số rất nhỏ chỉ 2,26% so với dữ liệu thực nghiệm Manzanares.
  • Xác định mức gia nhiệt dòng khí qua mái thu nhiệt đạt $\Delta T = 15,1\text{ K}$, bảo đảm động lực đối lưu nhiệt liên tục cho hệ thống.
  • Khẳng định tính ưu việt của giải pháp chia lưới 1/64 với 26.802 nút và 54.283 phần tử giúp tối ưu hóa tài nguyên tính toán chuyên sâu.
  • Đóng góp cơ sở phương pháp luận vững chắc để mở rộng phân tích tương tác kết cấu FSI cho các dự án tháp quy mô lớn 1.000 m trong giai đoạn 2026 - 2030.

Công trình nghiên cứu mở ra hướng tiếp cận công nghệ xanh đột phá và bền vững cho ngành năng lượng tái tạo. Hãy áp dụng ngay các giải pháp mô phỏng số tiên tiến này để nâng cao hiệu quả thiết kế và hiện thực hóa các dự án tháp năng lượng mặt trời tại Việt Nam!