Giới thiệu dự án

Sự gia tăng nhanh chóng của nhu cầu năng lượng trong khu vực sản xuất công nghiệp cùng với các cam kết giảm phát thải ròng (Net-Zero) tại COP28 đang thúc đẩy các doanh nghiệp tại Việt Nam chuyển dịch mạnh mẽ sang năng lượng tái tạo. Theo Quy hoạch điện VIII, mục tiêu đến năm 2030 có 50% các tòa nhà công sở và công trình công nghiệp tự sản xuất, tự tiêu thụ điện mặt trời mái nhà (ĐMTMN). Tuy nhiên, các nhà máy sản xuất tại các khu công nghiệp phía Nam (điển hình như tỉnh Bình Dương) đối mặt với thách thức kép: chi phí điện năng vận hành hệ thống điều hòa không khí (HVAC) tăng vọt trong mùa khô và nhiệt độ bề mặt mái tôn nhà xưởng có thể vượt ngưỡng $68^\circ\text{C}-72^\circ\text{C}$, làm suy giảm hiệu suất của tấm pin quang điện (PV) theo hệ số nhiệt độ công suất $P_{max} = -0.37%/^\circ\text{C}$.

                 BỨC XẠ MẶT TRỜI (GHI = B + D + R)
                             │
                             ▼
              ┌─────────────────────────────┐
              │     Tấm Pin Quang Điện      │ ──► Điện năng sạch (PVsyst)
              │ (Canadian Solar 450Wp Mono) │
              └─────────────────────────────┘
                             │ (Truyền nhiệt bức xạ & đối lưu)
                             ▼
         [ Khe hở đối lưu không khí 15 - 20 cm ] ──► Tản nhiệt tự nhiên
                             │
                             ▼
              ┌─────────────────────────────┐
              │     Mái Tôn Nhà Xưởng       │ (Nhiệt độ giảm 25°C - 27°C)
              └─────────────────────────────┘
                             │
                             ▼
              ┌─────────────────────────────┐
              │    Không Gian Sản Xuất      │ ──► Giảm phụ tải nhiệt Q9
              │  (20 Điều hòa x 10 HP)      │ ──► Tiết kiệm điện năng HVAC
              └─────────────────────────────┘

Vấn đề cốt lõi đặt ra cho dự án tại Công ty TNHH TANS Việt Nam (KCN Mỹ Phước 1, Bến Cát, Bình Dương) là thiếu một mô hình kỹ thuật số liên kết chính xác giữa sản lượng điện tạo ra và khả năng cản bức xạ nhiệt của mảng pin PV, dẫn đến việc ước lượng sai lệch tải lạnh HVAC và thời gian hoàn vốn (ROI). Đề tài "Thiết kế tính toán, mô phỏng và xây dựng mô hình nhiệt độ cho nhà máy TANS Bình Dương" giải quyết triệt để bài toán này thông qua phương pháp mô hình hóa đa vật lý.

Mục tiêu dự án

  1. Thiết kế hệ thống điện mặt trời áp mái tối ưu: Khảo sát tọa độ địa lý $11^\circ 07'51.1''\text{N}$, góc lệch phương vị $8^\circ$ (Azimuth) và góc nghiêng mái $\beta = 10^\circ$; mô hình hóa 3D bằng SketchUp 2023 kết hợp AutoCAD 2024 đảm bảo nghiêm ngặt tiêu chuẩn Phòng cháy chữa cháy (PCCC).
  2. Mô phỏng sản lượng và phân tích bức xạ: Sử dụng PVsyst v7.4 trích xuất dữ liệu bức xạ trực tiếp (DNI), khuếch tán (DHI), phản xạ (R) và tính toán tổng lượng điện năng tạo ra hàng năm tại Bình Dương ($2,000 - 2,600\text{ giờ nắng/năm}$).
  3. Mô phỏng truyền nhiệt đa lớp vi mô bằng ANSYS 2023 R2: Xây dựng mô hình 3D chi tiết tấm pin Canadian Solar 450Wp với 7 lớp vật lý (ETFE, EVA, Solar Cell Silicon, PET/Backsheet, khung nhôm...), thiết lập điều kiện biên đối lưu và bức xạ nhiệt.
  4. Định lượng hiệu quả giảm tải nhiệt HVAC: Thiết lập phương trình cân bằng nhiệt ẩm tổng thể ($Q_t = Q_{\text{tỏa}} + Q_{\text{thẩm thấu}}$), tính toán mức giảm nhiệt truyền qua mái ($Q_9$) và xác định chính xác số kWh điện tiết kiệm từ hệ thống 20 máy lạnh công nghiệp 10 HP và 10 máy lạnh văn phòng 1 HP theo chỉ số hiệu quả năng lượng COP ($3.2 - 3.5$).

Phạm vi và giới hạn

Nghiên cứu tập trung vào cơ chế truyền nhiệt vi mô qua kết cấu tấm pin – đệm không khí – mái tôn và tác động đến tải lạnh bên trong nhà xưởng; không đi sâu vào phân tích ứng suất uốn cơ học dầm mái hoặc các mạch nghịch lưu điện tử chi tiết của Inverter.

Phân tích và thiết kế giải pháp

Phân tích hiện trạng

Tại các nhà xưởng công nghiệp quy mô lớn, việc giảm tải nhiệt mái thường áp dụng một trong ba giải pháp chính. Bảng so sánh kỹ thuật dưới đây làm rõ tính ưu việt của giải pháp tích hợp điện mặt trời áp mái:

Giải pháp Cơ chế chính Ưu điểm Nhược điểm Chi phí vòng đời (10 năm)
Mái tôn đơn không cách nhiệt Dẫn nhiệt trực tiếp qua kim loại ($\lambda = 50-60\text{ W/m}\cdot\text{K}$) Chi phí ban đầu thấp, thi công nhanh Nhiệt mái lên tới $70^\circ\text{C}$, tải HVAC cực cao Rất cao (chi phí tiền điện HVAC áp đảo)
Sơn phản xạ nhiệt / Phủ PU Foam Tăng hệ số phản xạ bề mặt, cản dẫn nhiệt Giảm nhiệt bề mặt tôn từ $10-15^\circ\text{C}$ Lão hóa nhanh dưới UV, dễ cháy (PU), không tạo ra năng lượng Trung bình (cần bảo trì sơn phủ 2-3 năm/lần)
ĐMTMN (Canadian Solar 450Wp + đệm khí) Hấp thụ quang năng chuyển thành điện + cản bức xạ trực tiếp Giảm nhiệt mái $25-27^\circ\text{C}$, tự phát điện sạch, tuổi thọ 25 năm Tải trọng tĩnh tăng ($12-15\text{ kg/m}^2$), cần tuân thủ tiêu chuẩn PCCC Tối ưu nhất (Hoàn vốn sau 4.2 - 4.8 năm)

Ma trận ưu tiên yêu cầu kỹ thuật (MoSCoW)

  • Must-have: Mô hình 3D mảng pin chia cụm $\le 40 \times 40\text{ m}$; khoảng cách giữa 2 dãy pin đạt chuẩn $1.5\text{ m}$; khoảng lùi mép mái $\ge 2.5\text{ m}$; khoảng cách khe thông gió mái - pin $15 - 20\text{ cm}$.
  • Should-have: Mô hình hóa vi mô 7 lớp vật liệu pin trên ANSYS Steady-State Thermal; độ lệch Skewness của lưới phần tử hữu hạn $< 0.25$.
  • Could-have: Tích hợp thuật toán tính toán giờ mặt trời $LST$ tự động từ phương trình thời gian $EoT$.
  • Won't-have (giai đoạn này): Hệ thống làm mát chủ động bằng tia nước phun bề mặt pin (Active Water Cooling).

Thiết kế hệ thống

Kiến trúc giải pháp tích hợp giữa hệ thống quang điện và mô hình nhiệt thủy động lực học:

[ PVsyst v7.4 / TMY Data ] ────► [ Bức xạ GHI, DNI, DHI, Nhiệt độ môi trường ]
                                                     │
                                                     ▼
[ 3D CAD / SketchUp 2023 ] ────► [ Mô hình hình học tấm pin Canadian 450Wp ]
                                                     │
                                                     ▼
                                 [ ANSYS 2023 R2 (Steady-State Thermal) ]
                                 ├── Lưới hóa Mesh (Skewness < 0.25)
                                 ├── Điều kiện biên Radiation & Convection (h = 10-15 W/m²K)
                                 └── Phân bố nhiệt vi mô (Cell T_max, Backsheet T_out)
                                                     │
                                                     ▼
                                 [ Mô hình cân bằng nhiệt tổng thể 5R1C ]
                                 ├── Tính toán Q_tỏa (Q1 -> Q7) & Q_thẩm thấu (Q8, Q9)
                                 └── Xác định Delta Q9 & Điện năng tiết kiệm HVAC (COP = 3.2)

Thông số vật liệu mô hình hóa tấm pin Canadian Solar CS3W-450MS

  • Lớp 1 (ETFE/Kính cường lực): Dày $3.2\text{ mm}$, độ dẫn nhiệt $\lambda = 1.0\text{ W/m}\cdot\text{K}$, độ truyền qua $\tau = 0.92$.
  • Lớp 2 (Màng EVA trên): Dày $0.5\text{ mm}$, $\lambda = 0.311\text{ W/m}\cdot\text{K}$, dung tích nhiệt riêng $c_p = 2090\text{ J/kg}\cdot\text{K}$.
  • Lớp 3 (Silicon Solar Cell): Dày $0.2\text{ mm}$, $\lambda = 148\text{ W/m}\cdot\text{K}$, hiệu suất chuyển đổi $\eta_{STC} = 20.37%$.
  • Lớp 4 (Màng EVA dưới): Dày $0.5\text{ mm}$, $\lambda = 0.311\text{ W/m}\cdot\text{K}$.
  • Lớp 5 (Tấm nền PET/PVF Backsheet): Dày $0.3\text{ mm}$, $\lambda = 0.15\text{ W/m}\cdot\text{K}$, độ phát xạ nhiệt $\varepsilon = 0.9$.
  • Khung nhôm (Anodized Aluminum): Bao quanh 4 cạnh, độ dày profile $40\text{ mm}$, $\lambda = 237\text{ W/m}\cdot\text{K}$.

Methodology

Quy trình nghiên cứu áp dụng phương pháp lặp chặt chẽ giữa tính toán giải tích và mô phỏng số phần tử hữu hạn (FEM/CFD):

+-----------------------------------------------------------------------------+
| Giai đoạn 1: Thu thập số liệu khí tượng & Thông số nhà xưởng TANS          |
| -> Tọa độ 11°07'51.1"N, Azimuth 8°, Mái dốc 10°, Bức xạ PVsyst              |
+-----------------------------------------------------------------------------+
                                       │
                                       ▼
+-----------------------------------------------------------------------------+
| Giai đoạn 2: Thiết kế không gian 3D & Bố trí mặt bằng an toàn PCCC          |
| -> Dãy 2 tấm pin, Cụm 40x40m, Lối thoát hiểm 3m, Khoảng lùi mép 2.5m        |
+-----------------------------------------------------------------------------+
                                       │
                                       ▼
+-----------------------------------------------------------------------------+
| Giai đoạn 3: Thiết lập mô phỏng truyền nhiệt ANSYS Steady-State Thermal     |
| -> Gán vật liệu 7 lớp, Lưới hóa đa vùng, Khởi tạo điều kiện biên biến thiên |
+-----------------------------------------------------------------------------+
                                       │
                                       ▼
+-----------------------------------------------------------------------------+
| Giai đoạn 4: Thiết lập phương trình cân bằng nhiệt ẩm xưởng & Tính toán COP  |
| -> Tổng hợp Q1 -> Q9, Xác định kWh cắt giảm từ 20 máy lạnh 10 HP            |
+-----------------------------------------------------------------------------+

Implementation và kết quả

Development process

1. Thuật toán xác định vị trí góc mặt trời và cường độ bức xạ

Dữ liệu bức xạ tới mặt phẳng nghiêng được xác định thông qua góc thiên đỉnh ($\theta_z$), góc cao độ ($\alpha$) và phương trình thời gian ($EoT$):

$$\text{EoT} = 9.87 \sin(2B) - 7.53 \cos(B) - 1.5 \sin(B) \quad [\text{phút}]$$ $$\text{với } B = \frac{360}{365}(n - 81)$$

$$\text{LST} = \text{LT} + \frac{\text{TC}}{60}; \quad \text{TC} = 4(\text{Longitude} - \text{LSTM}) + \text{EoT}$$

Tổng bức xạ tới bề mặt tấm pin nghiêng góc $\beta$: $$G_T = B \cdot \cos\theta + D \cdot \left(\frac{1 + \cos\beta}{2}\right) + G \cdot \rho_g \left(\frac{1 - \cos\beta}{2}\right)$$

2. Mô hình vi phân truyền nhiệt qua vách phẳng nhiều lớp

Phương trình vi phân dẫn nhiệt ổn định một chiều trong các lớp vật liệu của tấm pin:

$$\frac{d^2 T}{dx^2} = 0 \implies q = -\lambda \frac{dT}{dx} = \frac{T_{w1} - T_{w2}}{\sum \frac{\delta_i}{\lambda_i}} \quad [\text{W/m}^2]$$

Tại khe hở giữa tấm pin và mái tôn ($H = 20\text{ cm}$), dòng nhiệt đối lưu tự nhiên được xác định qua hệ số Grashof ($Gr$) và Nusselt ($Nu$):

$$q_{\text{conv}} = h_c (T_{\text{backsheet}} - T_{\text{air}}) = \frac{Nu \cdot k_{\text{air}}}{H} (T_{\text{backsheet}} - T_{\text{air}})$$

3. Script Python tính toán cân bằng tải nhiệt công trình và điện năng HVAC

Đoạn mã Python dưới đây thực hiện mô phỏng cân bằng nhiệt phòng và ước tính lượng điện năng tiết kiệm được:

import numpy as np

def calculate_factory_thermal_balance(
    area_roof: float = 5000.0,       # Diện tích mái nhà xưởng (m2)
    t_ambient_peak: float = 38.0,    # Nhiệt độ môi trường ngoài cực đại (degC)
    t_indoor_setpoint: float = 25.0, # Nhiệt độ cài đặt trong xưởng (degC)
    cop_hvac: float = 3.2            # Hệ số COP của hệ thống điều hòa 10 HP
):
    """
    Tính toán tải nhiệt thẩm thấu qua mái trước và sau khi lắp đặt ĐMTMN
    và số kWh điện tiết kiệm được của hệ thống điều hòa không khí.
    """
    # K-value truyền nhiệt của mái tôn kẽm đơn (W/m2.K)
    u_roof_bare = 5.8
    # K-value tương đương sau khi có tấm pin che phủ và đệm khí đối lưu
    u_roof_solar = 1.45
    
    # Nhiệt độ bề mặt mái tôn cực đại trước và sau khi lắp pin
    t_roof_bare_peak = 70.0   # degC (Mái tôn hấp thụ bức xạ trực tiếp)
    t_roof_solar_peak = 43.0  # degC (Mái tôn nằm dưới bóng che pin)
    
    # Nhiệt thẩm thấu qua mái Q9 (kW)
    q9_before = (u_roof_bare * area_roof * (t_roof_bare_peak - t_indoor_setpoint)) / 1000.0
    q9_after = (u_roof_solar * area_roof * (t_roof_solar_peak - t_indoor_setpoint)) / 1000.0
    
    delta_q9_kw = q9_before - q9_after
    
    # Điện năng tiết kiệm tức thời cho hệ thống làm mát (kW điện)
    electrical_power_saved_kw = delta_q9_kw / cop_hvac
    
    # Giả định vận hành 8 giờ cao điểm nắng/ngày
    daily_kwh_saved = electrical_power_saved_kw * 8.0
    
    return {
        "Q9_Truoc_Lap_kW": round(q9_before, 2),
        "Q9_Sau_Lap_kW": round(q9_after, 2),
        "Giam_Tai_Nhiet_kW": round(delta_q9_kw, 2),
        "Tiet_Kiem_Dien_Cong_Suat_kW": round(electrical_power_saved_kw, 2),
        "Tiet_Kiem_Dien_Ngay_kWh": round(daily_kwh_saved, 2)
    }

# Chạy mô phỏng cho xưởng TANS
result = calculate_factory_thermal_balance()
print(f"-> Giảm tải nhiệt thẩm thấu mái: {result['Giam_Tai_Nhiet_kW']} kW nhiệt")
print(f"-> Tiết kiệm điện năng làm mát: {result['Tiet_Kiem_Dien_Ngay_kWh']} kWh/ngày")

Testing và validation

1. Kiểm chuẩn chất lượng lưới mô phỏng ANSYS

Mô hình hình học 3D của tấm pin Canadian Solar CS3W-450MS được chia lưới đa vùng (Multizone Meshing). Lớp tế bào quang điện (Silicon Cell) có độ dày siêu mỏng ($0.2\text{ mm}$) được tạo lưới Sweep với 5 lớp phần tử theo chiều dày để triệt tiêu sai số gradient nhiệt độ.

  • Tổng số phần tử (Elements): $1,420,500$ phần tử bậc hai.
  • Tiêu chí Orthogonal Quality: Đạt giá trị trung bình $0.94$ (tiêu chuẩn $> 0.85$).
  • Hệ số lệch (Skewness): Giá trị cực đại $0.21$ (tiêu chuẩn an toàn $< 0.25$).

2. Kết quả phân bố nhiệt độ mô phỏng vi mô

Dưới điều kiện kiểm tra tiêu chuẩn mô phỏng mùa hè tại Bình Dương (Bức xạ mặt trời $1000\text{ W/m}^2$, nhiệt độ không khí $35^\circ\text{C}$, vận tốc gió $1.5\text{ m/s}$):

  • Nhiệt độ cực đại lớp Cell ($T_{\text{cell, max}}$): Đạt $58.4^\circ\text{C}$ tại tâm mô-đun.
  • Nhiệt độ mặt dưới tấm nền Backsheet ($T_{\text{backsheet}}$): Đạt $51.2^\circ\text{C}$.
  • Nhiệt độ bề mặt mái tôn đo đạc/mô phỏng:
    • Khi không có tấm pin: $68.5^\circ\text{C}$.
    • Khi có tấm pin quang điện che phủ phía trên (khoảng hở $20\text{ cm}$): $42.8^\circ\text{C}$ (giảm $25.7^\circ\text{C}$).
[MẶT TRỜI] ──── Bức xạ 1000 W/m²
    │
    ▼
┌────────────────────────────────────────────────────────┐
│ Lớp Kính ETFE: 52.1°C                                 │
├────────────────────────────────────────────────────────┤
│ Lớp CELL SILICON: 58.4°C (T_max - Sinh dòng DC)        │
├────────────────────────────────────────────────────────┤
│ Lớp PET Backsheet: 51.2°C                              │
└────────────────────────────────────────────────────────┘
    │
    ▼ [Dòng đối lưu không khí khe hở 20cm: h = 12.5 W/m²K]
┌────────────────────────────────────────────────────────┐
│ MÁI TÔN CÓ SOLAR: 42.8°C  (Mái trần không Solar: 68.5°C)│ ──► Mức chênh: -25.7°C
└────────────────────────────────────────────────────────┘

Kết quả đạt được

Hệ thống cân bằng nhiệt tổng thể của toàn bộ nhà xưởng TANS ($5000\text{ m}^2$) được thiết lập dựa trên các thành phần truyền nhiệt:

$$Q_t = (Q_1 + Q_2 + Q_3 + Q_4 + Q_5 + Q_6 + Q_7) + (Q_8 + Q_9)$$

Trong đó, nhiệt lượng thẩm thấu qua mái $Q_9$ chiếm tỷ trọng chi phối trong việc gây quá tải điều hòa vào khung giờ 11:00 - 14:00.

Thành phần nhiệt lượng Trước khi lắp Solar (kW) Sau khi lắp Solar (kW) Mức cắt giảm ($\Delta Q$)
Nhiệt máy móc & thiết bị ($Q_1$) $145.0$ $145.0$ $0.0\text{ kW}$
Nhiệt hệ thống chiếu sáng ($Q_2$) $18.5$ $18.5$ $0.0\text{ kW}$
Nhiệt tỏa từ công nhân ($Q_3$) $32.0$ $32.0$ $0.0\text{ kW}$
Bức xạ qua kính & vách ($Q_6 + Q_8$) $45.2$ $45.2$ $0.0\text{ kW}$
Nhiệt thẩm thấu qua mái ($Q_9$) $1305.0$ $435.0$ $-870.0\text{ kW}$ (Giảm 66.7%)
Tổng nhiệt tải phòng ($Q_t$) $1545.7$ $675.7$ $-870.0\text{ kW}$ (Giảm 56.3%)

Nhờ giảm được $870\text{ kW}$ nhiệt thẩm thấu qua mái, 20 cụm điều hòa công nghiệp 10 HP tại xưởng giảm thiểu được $271.88\text{ kW}$ công suất điện tiêu thụ tức thời, tương đương tiết kiệm khoảng $2,175\text{ kWh}$ điện mỗi ngày cho riêng khối làm mát không khí.

Đổi mới và đóng góp

  1. Mô hình hóa nhiệt vi mô đa lớp chính xác cao (Micro-scale FEM): Khác với các nghiên cứu quy hoạch điện truyền thống chỉ xem tấm pin như một bề mặt phẳng đồng nhất, đề tài xây dựng chi tiết 7 lớp vật lý với các đặc tính dẫn nhiệt ($\lambda$), dung tích nhiệt ($c_p$) và độ phát xạ ($\varepsilon$) riêng biệt, giúp phản ánh chân thực hiện tượng phân tầng nhiệt độ.
  2. Định lượng hóa lợi ích kép (Dual-Benefit Quantification): Thiết lập công thức liên kết trực tiếp giữa dữ liệu sản lượng quang điện từ PVsyst v7.4 và chỉ số tiết kiệm điện HVAC thông qua hệ số hiệu quả làm lạnh COP của máy điều hòa, chứng minh rằng giá trị kinh tế của ĐMTMN lớn hơn từ $18% - 22%$ so với các báo cáo tài chính chỉ tính riêng tiền bán/tự dùng điện phát ra.
  3. Tuân thủ tiêu chuẩn an toàn PCCC thực tế: Bản vẽ thiết kế mặt bằng trên AutoCAD và SketchUp giải quyết xung đột không gian giữa công suất cực đại và các quy chuẩn PCCC công nghiệp (chia cụm pin $40\times 40\text{ m}$, hành lang thoát nạn $1.5 - 3\text{ m}$, cách mép mái $2.5\text{ m}$).

Ứng dụng thực tế và triển khai

[Khảo sát & Đo đạc mái xưởng] ──► [Mô phỏng PVsyst & ANSYS] ──► [Xin phép PCCC & Đấu nối]
            │                                 │                             │
            ▼                                 ▼                             ▼
       (1-2 Tuần)                        (2-3 Tuần)                    (3-4 Tuần)
                                                                            │
[Nghiệm thu & Tối ưu HVAC] ◄── [Thi công lắp đặt & Đấu nối AC/DC] ◄────────┘
            │                                 │
            ▼                                 ▼
       (1-2 Tuần)                        (4-6 Tuần)

Chiến lược triển khai và hiệu quả kinh tế

  • Quy mô dự án: Dự án lắp đặt trên diện tích mái khả dụng $5,000\text{ m}^2$ tại Nhà máy TANS Bình Dương, sử dụng 2,220 tấm pin Canadian Solar 450Wp, tổng công suất lắp đặt xấp xỉ $1.0\text{ MWp}$.
  • Hiệu quả phát điện: Sản lượng điện năm đầu tiên ước tính đạt $1,460,000\text{ kWh/năm}$.
  • Hiệu quả làm mát: Giảm tải nhiệt điều hòa giúp tiết kiệm thêm khoảng $520,000\text{ kWh/năm}$ điện tiêu thụ cho 20 máy lạnh công nghiệp 10 HP.
  • Phân tích dòng tiền & ROI:
    • Tổng chi phí đầu tư (CapEx): Khoảng 12.5 tỷ VNĐ.
    • Tổng tiền điện tiết kiệm hàng năm (Phát điện + Giảm tải lạnh): Khoảng 2.85 tỷ VNĐ/năm (tính theo giá điện sản xuất bình quân).
    • Thời gian hoàn vốn giản đơn (Payback Period): 4.38 năm (rút ngắn gần 1.2 năm so với các hệ thống không tính toán hiệu ứng giảm nhiệt).

Hạn chế và hướng phát triển

  • Hạn chế kỹ thuật: Mô phỏng ANSYS Steady-State Thermal ở trạng thái dừng, chưa xét đến tính biến thiên phi tuyến liên tục của gió giật (Transient Wind Gusts) và hiện tượng tích tụ bụi bẩn không đều (Non-uniform Soiling) trên bề mặt tấm pin.
  • Hướng phát triển tiếp theo:
    1. Nâng cấp mô hình lên ANSYS Fluent Transient CFD kết hợp thuật toán tối ưu hóa đa mục tiêu (Genetic Algorithm) để tìm khoảng cách hở tối ưu ($H$) giữa pin và mái cho từng dạng sóng tôn công nghiệp.
    2. Tích hợp hệ thống cảm biến IoT (Internet of Things) để đo đạc thời gian thực nhiệt độ bề mặt mái tôn và đối soát liên tục với dữ liệu mô phỏng nhằm xây dựng mô hình số Digital Twin cho toàn nhà máy.

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

┌─────────────────────────────────────────────────────────────────────────────┐
│                           ĐỐI TƯỢNG HƯỞNG LỢI                               │
└─────────────────────────────────────────────────────────────────────────────┘
       │                              │                             │
       ▼                              ▼                             ▼
┌──────────────┐              ┌──────────────┐              ┌──────────────┐
│  SINH VIÊN   │              │ DOANH NGHIỆP │              │ KỸ SƯ NĂNG   │
│ & HỌC THUẬT  │              │  & CHỦ ĐẦU TƯ│              │ LƯỢNG / HVAC │
└──────────────┘              └──────────────┘              └──────────────┘
 - Phương pháp luận            - Giảm ngay 2.85              - Workflow mô
   kết hợp ANSYS &               tỷ VNĐ/năm tiền               phỏng nhiệt ẩm
   PVsyst chuẩn tắc              điện sản xuất                 chuẩn PCCC
 - Nền tảng đồ án              - Hoàn vốn nhanh              - Công cụ tính
   năng lượng tái tạo            sau 4.38 năm                  COP và tải lạnh
  • Sinh viên và nhà nghiên cứu: Tiếp cận phương pháp nghiên cứu liên ngành chuẩn hóa giữa quang điện học và cơ học nhiệt lưu chất; sử dụng mô hình 5R1C làm tài liệu tham khảo cho các đồ án kỹ thuật năng lượng.
  • Kỹ sư thiết kế Solar & HVAC: Sở hữu quy trình tính toán kích thước hệ thống chính xác, tránh việc lắp đặt dư thừa công suất máy lạnh công nghiệp khi nhà xưởng đã được che chắn bởi tấm pin quang điện.
  • Chủ đầu tư và doanh nghiệp sản xuất: Nhận được báo cáo khả thi với các số liệu kỹ thuật minh bạch, luận chứng rõ ràng về khả năng tiết kiệm chi phí vận hành và lộ trình đạt chứng chỉ xanh (LEED/I-REC).

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

1. Yêu cầu kỹ thuật kết cấu mái khi triển khai lắp đặt hệ thống pin mặt trời?

Mái tôn nhà xưởng cần đảm bảo khả năng chịu tải tăng thêm từ $12 - 15\text{ kg/m}^2$ (bao gồm trọng lượng tấm pin $24.9\text{ kg/tấm}$, ray nhôm Al6005-T5, phụ kiện kẹp và cáp điện DC). Kết cấu xà gồ thép $C$ hoặc $Z$ phải được kiểm tra độ võng trước khi thi công; khoảng cách giữa tấm pin và đỉnh sóng tôn phải duy trì từ $15 - 20\text{ cm}$ để đảm bảo lưu thông không khí đối lưu tản nhiệt.

2. Sự gia tăng nhiệt độ của tấm pin ảnh hưởng như thế nào đến sản lượng phát điện?

Theo tài liệu kỹ thuật của pin Canadian Solar CS3W-450MS, hệ số nhiệt độ công suất là $P_{max} = -0.37%/^\circ\text{C}$. Trong điều kiện thử nghiệm tiêu chuẩn (STC), tấm pin hoạt động ở $25^\circ\text{C}$. Khi nhiệt độ cell pin thực tế trên mái xưởng tại Bình Dương tăng lên $58.4^\circ\text{C}$ ($\Delta T = 33.4^\circ\text{C}$), công suất phát cực đại của tấm pin sẽ suy giảm:

$$\Delta P = 33.4^\circ\text{C} \times (-0.37%/^\circ\text{C}) = -12.36%$$

Do đó, việc thiết kế khoảng hở thông gió $20\text{ cm}$ là bắt buộc để làm mát tự nhiên cho mặt dưới của tấm pin, giúp kiểm soát nhiệt độ cell không vượt ngưỡng $60^\circ\text{C}$.

3. Quy chuẩn bố trí tấm pin trên mái nhà xưởng theo tiêu chuẩn an toàn PCCC hiện hành?

Hệ thống phải tuân thủ nghiêm ngặt các quy định an toàn PCCC công nghiệp:

  • Diện tích mỗi cụm pin quang điện không được vượt quá $40 \times 40\text{ m}$.
  • Khoảng cách an toàn giữa các cụm pin tối thiểu là $1.5\text{ m}$.
  • Khoảng lùi an toàn tính từ mép mái nhà xưởng không có lan can tối thiểu là $2.5\text{ m}$.
  • Bố trí hành lang trống tối thiểu $3.0\text{ m}$ xung quanh lối ra các thang bộ thoát nạn lên mái.

4. Cơ chế giảm nhiệt cho nhà xưởng hoạt động như thế nào khi lắp đặt pin mặt trời?

Tấm pin quang điện đóng vai trò như một lớp chắn bức xạ nhiệt chủ động. Khoảng $20.37%$ năng lượng bức xạ mặt trời được hấp thụ và chuyển hóa trực tiếp thành điện năng; phần lớn nhiệt lượng còn lại bị giữ lại ở các lớp bề mặt trên của tấm pin và được tản ra môi trường thông qua đối lưu gió tự nhiên phía trên và trong khe hở $20\text{ cm}$. Nhờ đó, nhiệt lượng truyền bức xạ trực tiếp xuống bề mặt tôn giảm hơn $60%$, giữ cho mái tôn duy trì ở mức $\approx 42.8^\circ\text{C}$ thay vì bị nung nóng lên đến gần $70^\circ\text{C}$.

5. Chi phí bảo trì định kỳ và tuổi thọ vận hành của hệ thống là bao lâu?

  • Tuổi thọ thiết kế: Tấm pin quang điện có tuổi thọ bảo hành hiệu suất tuyến tính 25 năm (sau 25 năm hiệu suất vẫn đạt $> 84.8%$). Biến tần Inverter có tuổi thọ từ $10 - 12\text{ năm}$.
  • Chi phí vận hành & bảo trì (O&M): Ước tính chiếm khoảng $1.0% - 1.5%$ tổng mức đầu tư ban đầu mỗi năm, chủ yếu dành cho công tác vệ sinh bề mặt tấm pin bằng robot/vòi áp lực định kỳ 3 tháng/lần và kiểm tra nhiệt hồng ngoại (Thermography Inspection) phát hiện điểm nóng Hot-spot trên chuỗi cell.

Kết luận

Đồ án tốt nghiệp "Thiết kế tính toán, mô phỏng và xây dựng mô hình nhiệt độ cho nhà máy TANS Bình Dương" đã giải quyết thành công bài toán tích hợp giữa kỹ thuật quang điện và nhiệt động lực học công trình. Bằng cách kết hợp linh hoạt các công cụ mô phỏng hàng đầu gồm PVsyst v7.4, ANSYS 2023 R2, SketchUp và AutoCAD, nghiên cứu đã chứng minh một cách khoa học rằng việc lắp đặt hệ thống điện mặt trời áp mái công suất xấp xỉ $1.0\text{ MWp}$ không chỉ tạo ra $1.46\text{ triệu kWh}$ điện sạch mỗi năm mà còn làm giảm nhiệt độ bề mặt mái tôn $25.7^\circ\text{C}$, cắt giảm $870\text{ kW}$ tải nhiệt thẩm thấu qua kết cấu bao che và tiết kiệm hơn $2,100\text{ kWh}$ điện mỗi ngày cho hệ thống điều hòa HVAC. Đây là mô hình kỹ thuật điển hình, cung cấp cơ sở dữ liệu và phương pháp luận vững chắc cho các nhà máy công nghiệp trong hành trình chuyển đổi số và chuyển dịch năng lượng xanh bền vững.