Giới thiệu dự án

Chất lượng môi trường trong nhà (Indoor Environmental Quality - IEQ) và tiện nghi nhiệt (Thermal Comfort - TC) là hai yếu tố then chốt ảnh hưởng trực tiếp đến sức khỏe, trạng thái tâm sinh lý và năng suất làm việc của con người. Các nghiên cứu từ Hiệp hội Kỹ sư Nhiệt lạnh và Điều hòa Không khí Hoa Kỳ (ASHRAE) chỉ ra rằng nhân viên văn phòng dành trung bình từ 8 đến 10 giờ mỗi ngày trong không gian kín; việc tối ưu hóa vi khí hậu có thể gia tăng hiệu suất làm việc từ 10% đến 15% và giảm thiểu đáng kể các hội chứng bệnh nhà kín (Sick Building Syndrome - SBS). Tại Việt Nam, theo báo cáo chất lượng không khí năm 2022, nồng độ bụi mịn $\text{PM}_{2.5}$ và chỉ số ô nhiễm đang ở mức đáng báo động (xếp hạng 36/117 quốc gia), đặt ra yêu cầu cấp thiết về việc duy trì hệ thống cấp gió tươi kết hợp kiểm soát vi khí hậu trong nhà xưởng và khối nhà văn phòng điều hành.

Vấn đề kỹ thuật cốt lõi tại văn phòng điều hành nhà máy Dalat Milk (tỉnh Lâm Đồng) là tình trạng phân bố nhiệt độ và trường vận tốc không đồng đều khi vận hành cụm dàn lạnh âm trần cassette 4 hướng thổi thuộc hệ thống VRV. Hiện tượng gió lùa cục bộ (Draft Rate - DR vượt ngưỡng cho phép) tác động trực tiếp vào vùng làm việc của nhân viên, đi kèm với độ chênh lệch nhiệt độ không khí theo chiều thẳng đứng giữa đầu và mắt cá chân (Vertical Air Temperature Difference - VATD), gây cảm giác khó chịu và suy giảm tiện nghi nhiệt.

Mục tiêu cụ thể của dự án nghiên cứu bao gồm:

  1. Tính toán chính xác cân bằng nhiệt và phụ tải lạnh thừa của văn phòng theo phương pháp Carrier tiêu chuẩn.
  2. Xây dựng mô hình 3D hình học chi tiết và thiết lập mô phỏng động lực học dòng chảy (Computational Fluid Dynamics - CFD) trên phần mềm ANSYS Workbench 2022 R1 (Fluent).
  3. Ứng dụng phương pháp quy hoạch thực nghiệm Taguchi với ma trận trực giao $L_9$ trên phần mềm Minitab để khảo sát ảnh hưởng đồng thời của 3 thông số: góc mở cánh đảo gió cassette ($30^\circ, 40^\circ, 50^\circ$), lưu lượng gió cấp ($810 - 1230\text{ m}^3/\text{h}$) và lưu lượng gió tươi ($200, 400, 600\text{ m}^3/\text{h}$).
  4. Đánh giá định lượng mức độ tiện nghi nhiệt theo các chỉ số tiêu chuẩn quốc tế ISO 7730 và ASHRAE Standard 55 (PMV, PPD, DR, VATD) cùng chất lượng không khí thông qua nồng độ $\text{CO}_2$.
  5. Đề xuất giải pháp cấu hình vận hành tối ưu và phương án cải tiến hệ thống phân phối không khí cho công trình.

Phương pháp tiếp cận dựa trên mô phỏng số 3D CFD kết hợp quy hoạch thực nghiệm Taguchi giúp rút ngắn thời gian thử nghiệm thực tế, giảm chi phí vận hành thử nghiệm và dự đoán chính xác trường phân bố vận tốc, nhiệt độ trong không gian làm việc.

Dự án đặt ra các chỉ số đầu ra kỳ vọng có thể lượng hóa:

  • Chỉ số dự đoán trung bình theo phiếu đánh giá (Predicted Mean Vote - PMV) nằm trong vùng tiện nghi: $-0.5 \le \text{PMV} \le +0.5$.
  • Tỷ lệ phần trăm không hài lòng dự đoán (Predicted Percentage of Dissatisfied - PPD) đạt mức tối ưu: $\text{PPD} < 10%$.
  • Tỷ lệ khó chịu do gió lùa (Draft Rate - DR) tại vùng hoạt động của con người: $\text{DR} \le 15%$ (thấp hơn giới hạn tiêu chuẩn $20%$).
  • Độ chênh lệch nhiệt độ theo chiều đứng (VATD) từ độ cao mắt cá chân ($0.1\text{ m}$) đến đầu người ngồi ($1.1\text{ m}$): $\text{VATD} < 3.0^\circ\text{C}$.
  • Nồng độ $\text{CO}_2$ trung bình trong phòng được duy trì dưới ngưỡng $1000\text{ ppm}$.

Phạm vi nghiên cứu tập trung vào không gian văn phòng tầng 2 thuộc nhà máy Dalat Milk tại Lâm Đồng với diện tích sàn $36\text{ m}^2$ ($8.0\text{ m} \times 4.5\text{ m} \times 3.0\text{ m}$), điều kiện biên khí hậu ngoài trời tại vĩ độ $12^\circ\text{B}$ ($t_N = 30^\circ\text{C}$, độ ẩm $\varphi_N = 83%$) và thông số thiết kế trong nhà ($t_T = 23^\circ\text{C}$, $\varphi_T = 60%$). Giới hạn của đề tài là khảo sát dòng chảy dừng (steady-state), giả định mật độ nhân sự và công suất thiết bị văn phòng hoạt động ổn định.

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

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

Trong thiết kế hệ thống HVAC truyền thống, các kỹ sư thường dựa trên phương pháp tính toán phụ tải tập trung (lumped thermal load calculation) và bố trí miệng gió theo kinh nghiệm. Cách tiếp cận này bộc lộ nhiều hạn chế khi không thể dự đoán được các vùng xoáy quẩn khí động, hiện tượng ngắn mạch luồng không khí (air short-circuiting) hoặc các điểm tập trung nhiệt cục bộ.

Phương pháp đánh giá Ưu điểm Nhược điểm Chi phí & Thời gian Độ chính xác không gian
Tính toán truyền thống (Carrier / TCVN 5687) Nhanh chóng, đơn giản, xác định tổng công suất lạnh tổng thể. Không hiển thị trường nhiệt độ, vận tốc 3D; không đánh giá được DR, VATD cục bộ. Chi phí rất thấp; Thời gian thực hiện: 1 - 2 ngày. Rất thấp (chỉ cho giá trị trung bình toàn phòng).
Đo đạc thực nghiệm hiện trường Đo lường dữ liệu thực tế tại các điểm cảm biến vật lý. Chi phí cảm biến cao; can thiệp vào vận hành; khó thay đổi nhiều kịch bản vận hành. Chi phí cao; Thời gian thực hiện: 2 - 4 tuần. Trung bình (phụ thuộc vào mật độ đặt đầu đo sensor).
Mô phỏng số 3D CFD (ANSYS Fluent) Trực quan hóa toàn diện trường 3D (áp suất, vận tốc, nhiệt độ, $\text{CO}_2$); thử nghiệm linh hoạt đa kịch bản. Yêu cầu năng lực tính toán phần cứng mạnh; đòi hỏi kỹ thuật chia lưới và gán điều kiện biên chuẩn xác. Chi phí phần mềm/nhân lực chuyên sâu; Thời gian: 3 - 7 ngày. Rất cao (độ phân giải đến từng phần tử thể tích hữu hạn).

Yêu cầu kỹ thuật của bài toán được xác lập theo mô hình phân loại MoSCoW:

  • Must have: Đảm bảo công suất lạnh bù đắp toàn bộ nhiệt hiện ($Q_{ht}$) và nhiệt ẩn ($Q_{at}$); các chỉ số PMV, PPD, DR đạt chuẩn ISO 7730 loại B.
  • Should have: Tối ưu hóa góc nghiêng cánh đảo gió dàn lạnh cassette để loại bỏ hiện tượng gió thổi trực tiếp vào vùng đầu của nhân viên làm việc; duy trì gradient nhiệt độ thẳng đứng $\text{VATD} < 2^\circ\text{C}$.
  • Could have: Tích hợp hệ thống phân phối gió tươi qua dàn xử lý sơ bộ PAU (Primary Air Unit) đồng bộ với lưu lượng hút khí thải.
  • Won't have: Mô phỏng quá trình truyền nhiệt động học phi dừng (transient dynamic) có sự di chuyển liên tục của con người trong phòng.

Thách thức kỹ thuật nằm ở việc xử lý lớp biên nhiệt sát bề mặt kính tiếp xúc trực tiếp bức xạ mặt trời hướng Tây Bắc và Đông Bắc (diện tích kính lớn $33.75\text{ m}^2$) với kính Calorex màu xanh ($\alpha_k = 0.57$) và màn che sáng ($\alpha_m = 0.56$), tạo nên dòng nhiệt bức xạ tức thời $Q_{11}$ phức tạp tác động vào không gian.

Thiết kế hệ thống

Kiến trúc giải pháp mô phỏng số bao gồm chuỗi liên kết tính toán vật lý và xử lý dữ liệu vi khí hậu:

[Mô hình Hình học CAD 3D (SpaceClaim/DesignModeler)]
                         │
                         ▼
[Chia lưới thể tích Poly-Hexcore (Fluent Meshing 2022 R1)]
                         │
                         ▼
[Thiết lập mô hình Vật lý & Điều kiện biên (ANSYS Fluent)]
 ├── Phương trình Năng lượng (Energy Equation)
 ├── Mô hình Rối k-epsilon Realizable / Enhanced Wall Treatment
 └── Phương trình Truyền chất Species Transport (Nồng độ CO2)
                         │
                         ▼
[Quy hoạch thực nghiệm Taguchi L9 (Minitab 21)]
 ├── Yếu tố A: Góc mở cánh đảo gió (30°, 40°, 50°)
 ├── Yếu tố B: Lưu lượng gió cấp cassette (810, 960, 1230 m³/h)
 └── Yếu tố C: Lưu lượng gió tươi PAU (200, 400, 600 m³/h)
                         │
                         ▼
[Trích xuất & Hậu xử lý dữ liệu vi khí hậu (CFD Post)]
 └── Tính toán các chỉ số: PMV, PPD, DR, VATD, CO2

Hệ thống công cụ và công nghệ sử dụng trong nghiên cứu:

  • Phần mềm mô phỏng CFD: ANSYS Mechanical Workbench 2022 R1 / ANSYS Fluent 2022 R1.
  • Phần mềm phân tích thống kê: Minitab v21.1.0 (sử dụng module Taguchi Design of Experiments và ANOVA Regression).
  • Tiêu chuẩn áp dụng: ASHRAE Standard 55-2020, ISO 7730:2005, TCVN 5687:2010.

Công thức xác định phụ tải lạnh tổng thể theo phương pháp Carrier:

$$Q_0 = Q_t = Q_{ht} + Q_{at}$$

Dòng nhiệt bức xạ qua kính có kể đến hệ số tác dụng tức thời $\eta_t = 0.81$:

$$Q_{11} = \eta_t \cdot F \cdot R_K \cdot c \cdot d_s \cdot m_m \cdot k_h \cdot \alpha_m \cdot \alpha_r$$

Trong đó dòng nhiệt bức xạ vào phòng qua kính Calorex $R_K = 143\text{ W/m}^2$, diện tích kính hướng Tây Bắc $F_{TB} = 21.6\text{ m}^2$ cho ra $Q_{11,TB} \approx 2502.2\text{ W}$.

Công thức xác định tỷ lệ người khó chịu do gió lùa (Draft Rate - DR):

$$\text{DR} = \left(34 - t_a\right) \left(v - 0.05\right)^{0.62} \left(0.37 \cdot v \cdot Tu + 3.14\right) \quad [%]$$

Trong đó:

  • $t_a$: Nhiệt độ không khí tại điểm khảo sát ($^\circ\text{C}$).
  • $v$: Vận tốc không khí trung bình cục bộ ($\text{m/s}$), với điều kiện $v < 0.05\text{ m/s}$ thì gán $\text{DR} = 0%$.
  • $Tu$: Cường độ dòng chuyển động quẩn (Turbulence Intensity, tính bằng $%$).

Chỉ số PMV được tính toán theo mô hình cân bằng nhiệt Fanger:

$$\text{PMV} = \left[0.303 \cdot e^{-0.036 \cdot M} + 0.028\right] \cdot L$$

$$\text{PPD} = 100 - 95 \cdot \exp\left(-\left(0.03353 \cdot \text{PMV}^4 + 0.2179 \cdot \text{PMV}^2\right)\right)$$

Methodology

Dự án áp dụng phương pháp nghiên cứu kết hợp giữa tính toán nhiệt kỹ thuật truyền thống, mô phỏng số thể tích hữu hạn (Finite Volume Method - FVM) và quy hoạch thực nghiệm Taguchi nhằm tối ưu hóa các tham số vận hành.

Giai đoạn 1: Khảo sát & Tính toán nhiệt (Tuần 1 - 3)
 ├── Thu thập bản vẽ kiến trúc, thông số kết cấu tường 230mm, kính Calorex
 └── Tính toán phụ tải lạnh Carrier (Q0 = Qht + Qat)
        │
Giai đoạn 2: Xây dựng Mô hình Hình học & Lưới 3D (Tuần 4 - 6)
 ├── Dựng CAD 3D phòng 8m x 4.5m x 3m với cassette và miệng gió tươi/thải
 └── Khảo sát tính độc lập của lưới Poly-hexcore (1.0M, 1.4M, 1.9M phần tử)
        │
Giai đoạn 3: Thực hiện Mô phỏng CFD & Ma trận Taguchi L9 (Tuần 7 - 10)
 ├── Thiết lập điều kiện biên (Velocity-inlet, Pressure-outlet, Heat flux)
 └── Chạy 9 ca mô phỏng hội tụ trên ANSYS Fluent 2022 R1
        │
Giai đoạn 4: Phân tích ANOVA & Đánh giá Tiện nghi (Tuần 11 - 12)
 ├── Nhập dữ liệu DR, VATD, PMV/PPD vào Minitab
 └── Phân tích hồi quy, xác định mức ảnh hưởng của từng thông số
        │
Giai đoạn 5: Đề xuất Giải pháp Cải tiến & Kết luận (Tuần 13 - 14)
 └── Đánh giá phương án thay thế cassette bằng dàn lạnh giấu trần nối ống gió

Bảng quản trị rủi ro kỹ thuật:

  • Rủi ro sai lệch hội tụ số: Hiện tượng Residuals dao động không đạt tiêu chí $10^{-4}$ cho phương trình xung lượng và $10^{-6}$ cho phương trình năng lượng $\rightarrow$ Giải pháp: Ứng dụng thuật toán ghép áp suất - vận tốc COUPLED kết hợp thuật toán chia lưới Poly-Hexcore để nâng cao chất lượng Orthogonal Quality $> 0.2$ và giảm Skewness $< 0.8$.
  • Rủi ro sai lệch điều kiện biên nhiệt: Tường tiếp xúc gián tiếp và trực tiếp có hệ số dẫn nhiệt khác nhau $\rightarrow$ Giải pháp: Tính toán chi tiết từng lớp vật liệu (vữa trát ngoài $15\text{ mm}$, gạch $200\text{ mm}$, vữa trong $15\text{ mm}$, $k_{tt} = 2.81\text{ W/m}^2\text{K}$).

Implementation và kết quả

Development process

Quá trình chia lưới thể tích đóng vai trò quyết định đến độ chính xác của nghiệm CFD. Đề tài tiến hành so sánh 3 cấu trúc lưới: Tetrahedral/Hexcore, Polyhedral và Poly-Hexcore. Cấu trúc Poly-Hexcore được lựa chọn nhờ khả năng giảm 40% số lượng phần tử so với lưới tứ diện truyền thống nhưng vẫn bảo toàn độ chính xác cao tại các vùng biên gradient vận tốc lớn.

Ba bộ lưới với mật độ khác nhau được thiết lập để kiểm tra sự hội tụ và tính độc lập của lưới:

  • Lưới thô (Mesh 1): Xấp xỉ $1.0 \times 10^6$ phần tử.
  • Lưới trung bình (Mesh 2): Xấp xỉ $1.4 \times 10^6$ phần tử.
  • Lưới mịn (Mesh 3): Xấp xỉ $1.9 \times 10^6$ phần tử.

Chỉ số Skewness tối đa được kiểm soát ở mức $0.78$ (nằm trong ngưỡng an toàn $< 0.85$) và Aspect Ratio trung bình duy trì ở mức $1.85$. Kết quả sai lệch vận tốc và nhiệt độ giữa Lưới 2 ($1.4\text{M}$) và Lưới 3 ($1.9\text{M}$) nhỏ hơn $1.2%$, do đó Lưới 2 ($1.4$ triệu phần tử) được chọn nhằm tối ưu hóa chi phí tính toán.

Thiết lập biểu thức Custom Field Function trong ANSYS Fluent để tính toán trường Draft Rate:

; Fluent Scheme / Expression setup for Draft Rate (DR) Calculation
(cx-macro-define 'calculate-draft-rate
  '(lambda ()
     (ti-menu-load-string 
       "define/custom-field-functions/define \"draft-rate\" 
        MAX(0.0, (34.0 - temp_c) * pow(MAX(0.0, mag(v) - 0.05), 0.62) * (0.37 * mag(v) * turb_intensity + 3.14))"
     )
   )
)

Điều kiện biên gán cho mô hình:

  • Miệng thổi cassette (Cassette Inlets - 4 hướng): Velocity-inlet với vận tốc gió tương ứng theo các cấp lưu lượng $810, 960, 1230\text{ m}^3/\text{h}$, nhiệt độ gió cấp $t_c = 16^\circ\text{C}$, góc lệch cánh hướng dòng so với mặt phẳng ngang lần lượt là $30^\circ, 40^\circ, 50^\circ$.
  • Miệng gió tươi (PAU Inlet): Velocity-inlet với lưu lượng $200, 400, 600\text{ m}^3/\text{h}$, nhiệt độ $t_{fresh} = 20^\circ\text{C}$, nồng độ $\text{CO}_2 = 400\text{ ppm}$.
  • Miệng hút hồi cassette & Hút thải: Pressure-outlet với áp suất tĩnh $P_{gauge} = 0\text{ Pa}$.
  • Nguồn nhiệt nội vi: Nhân viên văn phòng phát tán nhiệt hiện $Q_{4h} = 75\text{ W/người}$, nồng độ phát thải $\text{CO}_2 = 0.005\text{ l/s/người}$; thiết bị máy tính phát nhiệt $Q_3 = 150\text{ W/bộ}$.

Testing và validation

Ma trận trực giao Taguchi $L_9(3^3)$ được thiết lập nhằm khảo sát 9 tổ hợp tham số:

Trường hợp Yếu tố A: Góc mở cánh đảo ($^\circ$) Yếu tố B: Lưu lượng gió cấp ($\text{m}^3/\text{h}$) Yếu tố C: Lưu lượng gió tươi ($\text{m}^3/\text{h}$) $\text{DR}_{\text{max}}$ Mặt cắt 1 (%) $\text{VATD}_{\text{max}}$ Mặt cắt 1 ($^\circ\text{C}$) $\text{PMV}_{\text{tb}}$ Vùng làm việc
Case 1 30 810 200 12.4 1.12 -0.15
Case 2 30 960 400 14.8 1.25 -0.22
Case 3 30 1230 600 18.2 1.41 -0.38
Case 4 40 810 400 16.5 1.38 -0.10
Case 5 40 960 600 21.3 1.56 -0.26
Case 6 40 1230 200 25.7 1.82 -0.42
Case 7 50 810 600 22.1 1.64 +0.05
Case 8 50 960 200 26.8 1.95 -0.18
Case 9 50 1230 400 31.4 2.24 -0.35

Kết quả phân tích hồi quy phương sai (ANOVA) trên phần mềm Minitab đối với chỉ số gió lùa DR và gradient nhiệt độ VATD tại hai mặt cắt khảo sát (Mặt cắt 1: ngang vị trí ngồi làm việc; Mặt cắt 2: dọc trục dàn lạnh) cho thấy:

  • Góc mở cánh đảo gió là nhân tố chi phối mạnh nhất đến chỉ số DR với tỷ lệ đóng góp ảnh hưởng $\text{Contribution} = 64.82%$ ($p\text{-value} = 0.003 < 0.05$).
  • Lưu lượng gió cấp đóng góp $27.15%$ vào sự biến thiên của DR ($p\text{-value} = 0.018$).
  • Lưu lượng gió tươi có mức độ ảnh hưởng không đáng kể đến DR ($8.03%$) nhưng đóng vai trò quyết định đến nồng độ $\text{CO}_2$.

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

Đồ thị phân bố vi khí hậu chỉ ra rằng:

  1. Chỉ số DR: Ở góc mở cánh đảo $30^\circ$, luồng không khí lạnh bám sát trần (hiệu ứng Coanda) trước khi hòa trộn và rơi êm dịu vào vùng làm việc, duy trì $\text{DR} < 15%$ tại đại đa số vị trí ngồi. Ngược lại, ở góc mở $50^\circ$, luồng gió lạnh bị thổi dốc trực tiếp xuống sàn, tạo nên vùng xoáy cục bộ với vận tốc $> 0.35\text{ m/s}$, đẩy chỉ số DR vượt $25% - 31.4%$, vi phạm tiêu chuẩn tiện nghi nhiệt ASHRAE 55.
  2. Gradient nhiệt độ thẳng đứng (VATD): Dao động từ $1.12^\circ\text{C}$ đến $2.24^\circ\text{C}$ trên toàn bộ các ca mô phỏng, hoàn toàn thỏa mãn yêu cầu tiêu chuẩn ($\text{VATD} < 3^\circ\text{C}$).
  3. Chỉ số PMV và PPD: Tại 9 điểm đo phân bố đều trong văn phòng, chỉ số PMV dao động từ $-0.38$ đến $+0.05$, tương ứng với tỷ lệ PPD dưới $8.5%$, chứng minh môi trường nhiệt đạt cấp độ tiện nghi cao nhất (Class A theo ISO 7730).
  4. Nồng độ $\text{CO}_2$: Với lưu lượng gió tươi $400\text{ m}^3/\text{h}$ và $600\text{ m}^3/\text{h}$, nồng độ $\text{CO}_2$ trong phòng duy trì ổn định ở mức $550 - 680\text{ ppm}$, đáp ứng hoàn hảo yêu cầu thông gió bảo vệ sức khỏe.

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

Nghiên cứu mang lại những đóng góp học thuật và kỹ thuật thực tiễn:

  1. Kết hợp đồng bộ CFD - Taguchi - ANOVA: Khác với các phương pháp thử - sai (trial-and-error) tốn kém, việc tích hợp ma trận Taguchi $L_9$ giúp cắt giảm $66.7%$ số lượng ca tính toán cần thiết (từ $3^3 = 27$ ca mô phỏng xuống còn 9 ca) mà vẫn rút ra được mô hình hồi quy đa biến có độ tin cậy $R^2 > 96.5%$.
  2. Định lượng hóa quy luật khí động của cánh đảo gió cassette: Làm sáng tỏ cơ chế tác động của góc mở cánh gió đến sự hình thành vi xoáy và hiện tượng bẫy nhiệt cục bộ trong phòng làm việc kín có tải bức xạ mặt trời lớn.
  3. Đề xuất chuyển đổi công nghệ phân phối khí: Đồ án chứng minh rằng mặc dù máy lạnh cassette có thể cải thiện tiện nghi bằng cách đặt góc mở $30^\circ$, nhưng giải pháp tối ưu bền vững nhất cho các văn phòng hiện đại là sử dụng dàn lạnh âm trần nối ống gió (Ducted Concealed Unit) kết hợp miệng thổi khuếch tán dạng khe (Linear Slot Diffuser) hoặc miệng khuếch tán vuông (Square Diffuser). Giải pháp này giúp phân phối trường nhiệt đều hơn $45%$, triệt tiêu hoàn toàn hiện tượng gió thổi dốc và nâng cao tính thẩm mỹ kiến trúc.
Tiêu chí kỹ thuật Thiết kế truyền thống (Kinh nghiệm) Mô phỏng thực nghiệm cục bộ Đề xuất tối ưu (CFD + Dàn giấu trần)
Kiểm soát gió lùa ($\text{DR}$) Không kiểm soát ($\text{DR} > 25%$) Tối ưu góc mở cassette ($\text{DR} \approx 14.8%$) Phân phối qua ống gió ($\text{DR} < 10%$)
Độ đồng đều nhiệt độ ($\Delta T$) Lệch cục bộ $2.5 - 3.5^\circ\text{C}$ Chênh lệch $1.2 - 1.8^\circ\text{C}$ Chênh lệch $< 0.8^\circ\text{C}$ toàn phòng
Chất lượng không khí ($\text{CO}_2$) Cấp tươi không kiểm soát ($> 1200\text{ ppm}$) Kiểm soát qua PAU ($550 - 680\text{ ppm}$) Tích hợp PAU và bộ lọc $\text{PM}_{2.5}$ ($< 600\text{ ppm}$)
Độ thỏa mãn tiện nghi ($\text{PPD}$) $18% - 25%$ $7% - 10%$ $< 6%$ (Tiện nghi tối ưu)

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

Kịch bản ứng dụng thực tế

Giải pháp từ nghiên cứu có khả năng áp dụng trực tiếp cho các khối nhà văn phòng nhà máy chế biến thực phẩm, dược phẩm, sữa (tương tự quy chuẩn của Dalat Milk, Vinamilk, TH True Milk) - nơi đòi hỏi kiểm soát nghiêm ngặt cả về vi khí hậu phòng làm việc lẫn phòng kiểm nghiệm chất lượng (QA/QC Laboratory).

Chiến lược triển khai và vận hành

  1. Cấu hình vận hành tức thời cho hệ thống hiện hữu:
    • Khóa cố định góc nghiêng cánh đảo gió dàn cassette ở mức $30^\circ$ đối với chế độ làm lạnh.
    • Cài đặt vận tốc quạt dàn lạnh ở mức trung bình (tương ứng lưu lượng $960\text{ m}^3/\text{h}$ cho phòng $36\text{ m}^2$).
    • Duy trì hệ thống quạt cấp gió tươi PAU ở mức lưu lượng $400\text{ m}^3/\text{h}$ trong suốt giờ làm việc.
  2. Kế hoạch nâng cấp dài hạn:
    • Khi cải tạo công trình, chuyển đổi sang dàn lạnh giấu trần nối ống gió công suất $5.6\text{ kW}$ ($18,000\text{ BTU/h}$) tích hợp 4 miệng thổi khuếch tán bố trí đối xứng.

Phân tích hiệu quả kinh tế (Cost-Benefit Analysis)

  • Chi phí đầu tư bổ sung: Việc tối ưu hóa góc đảo gió và lưu lượng trên hệ thống hiện hữu không phát sinh chi phí phần cứng (0 VNĐ).
  • Tiết kiệm năng lượng: Thiết lập góc thổi $30^\circ$ tận dụng hiệu quả đối lưu tự nhiên giúp nâng nhiệt độ cài đặt từ $22^\circ\text{C}$ lên $24^\circ\text{C}$ mà vẫn đảm bảo độ tiện nghi tương đương, giúp tiết kiệm $12% - 16%$ điện năng tiêu thụ của hệ thống nén lạnh VRV.
  • Thời gian hoàn vốn (ROI): Khi đầu tư hệ thống giấu trần nối ống gió mới, chi phí tăng thêm khoảng $15 - 20$ triệu VNĐ/phòng sẽ được hoàn vốn trong vòng 1.8 năm nhờ cắt giảm điện năng và gia tăng năng suất lao động.

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

Nghiên cứu đã giải quyết trọn vẹn mục tiêu đề ra nhưng vẫn tồn tại một số hạn chế kỹ thuật:

  • Giả định dòng chảy dừng (Steady-state): Chưa xét đến tác động tức thời của việc đóng/mở cửa ra vào hoặc sự di chuyển qua lại của nhân viên.
  • Mô hình bức xạ đơn giản hóa: Chưa tích hợp đầy đủ mô hình bức xạ bề mặt Discrete Ordinates (DO) đối với toàn bộ các chi tiết nội thất phức tạp.

Các hướng nghiên cứu tiếp nối trong tương lai:

  1. Phát triển mô hình CFD động học (Transient CFD): Khảo sát quá trình khuếch tán sol khí (aerosol) và vi khuẩn trong không gian kín nhằm phục vụ công tác phòng chống dịch bệnh đường hô hấp.
  2. Tích hợp thuật toán điều khiển thông minh AI/IoT: Kết nối dữ liệu mô phỏng với hệ thống quản lý tòa nhà (BMS) để tự động điều chỉnh cánh đảo gió và van gió tươi VAV theo thời gian thực dựa trên cảm biến $\text{CO}_2$ và nhiệt độ IoT.

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

  • Sinh viên & Học viên chuyên ngành Kỹ thuật Nhiệt - HVAC: Tiếp cận quy trình chuẩn về việc kết hợp phương pháp tính toán tải nhiệt Carrier với mô phỏng CFD 3D và xử lý thống kê Minitab Taguchi.
  • Kỹ sư tư vấn thiết kế MEP/HVAC: Sở hữu bảng dữ liệu thực nghiệm chuẩn xác về góc mở cánh đảo gió cassette và lưu lượng khí để áp dụng vào các dự án văn phòng thực tế.
  • Chủ đầu tư & Doanh nghiệp: Nâng cao tiện nghi môi trường làm việc, bảo vệ sức khỏe nhân viên, giảm tỷ lệ vắng mặt do ốm đau và tối ưu hóa chi phí điện năng vận hành hàng tháng từ $12% - 16%$.
  • Nhà nghiên cứu khoa học: Cung cấp cơ sở dữ liệu vi khí hậu kiểm chứng cho các công trình nghiên cứu về tiện nghi nhiệt trong điều kiện khí hậu nhiệt đới gió mùa tại Việt Nam.

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

1. Cấu hình phần cứng và phần mềm yêu cầu để thực hiện mô phỏng CFD tương tự là gì?

Hệ thống cần tối thiểu CPU 8 nhân / 16 luồng (Intel Core i7/i9 hoặc AMD Ryzen 7/9), 32 GB RAM (khuyến nghị 64 GB cho lưới $> 2$ triệu phần tử) và card đồ họa chuyên dụng hỗ trợ tăng tốc tính toán CUDA. Phần mềm sử dụng là ANSYS Workbench 2022 R1 (module Fluent Meshing & Fluent Solver) và Minitab phiên bản 19 trở lên.

2. Tại sao góc mở cánh đảo gió $30^\circ$ lại mang lại chỉ số DR tốt hơn góc $50^\circ$?

Ở góc nghiêng $30^\circ$ so với phương ngang, luồng gió lạnh tốc độ cao bám dính vào trần nhờ hiệu ứng Coanda, giúp luồng khí có đủ thời gian hòa trộn với khối không khí ấm trong phòng và giảm dần vận tốc trước khi rơi vào vùng hoạt động của con người. Ở góc $50^\circ$, luồng gió bị dội trực tiếp xuống sàn với vận tốc lớn ($> 0.35\text{ m/s}$), gây hiện tượng gió lùa lạnh buốt (Draft Rate cao).

3. Phương pháp Taguchi $L_9$ có đảm bảo độ tin cậy so với mô phỏng toàn phần không?

Có. Ma trận trực giao $L_9$ phân bố đều các mức thử nghiệm trong không gian tham số 3 yếu tố x 3 mức. Khi kết hợp với phân tích hồi quy ANOVA, phương pháp cho phép đánh giá mức độ tương quan chính xác với hệ số xác định $R^2 > 96%$, giúp tiết kiệm 18 ca mô phỏng mà không làm mất tính tổng quát của kết quả.

4. Chi phí bảo trì khi vận hành hệ thống cấp gió tươi PAU là bao nhiêu?

Hệ thống PAU chỉ yêu cầu vệ sinh hoặc thay thế phin lọc bụi định kỳ (phin lọc sơ cấp G4 vệ sinh 3 tháng/lần, phin lọc tinh F7/F8 thay thế 6 - 12 tháng/lần) với chi phí trung bình khoảng $1.5 - 3$ triệu VNĐ/năm cho một không gian văn phòng tiêu chuẩn.

5. Làm thế nào để áp dụng kết quả này cho các phòng có diện tích khác $36\text{ m}^2$?

Kỹ sư có thể sử dụng nguyên lý đồng dạng khí động và giữ nguyên mật độ thông gió: duy trì bội số tuần hoàn gió làm lạnh từ $8 - 12$ lần/giờ, lưu lượng gió tươi tối thiểu $25 - 30\text{ m}^3/\text{h/người}$, và cài đặt góc thổi cassette ở $30^\circ$ nhằm đảm bảo vận tốc gió trong vùng làm việc luôn nằm trong dải tối ưu $0.15 - 0.25\text{ m/s}$.

Kết luận

Đồ án tốt nghiệp "Nghiên cứu các điều kiện tiện nghi nhiệt bằng phương pháp mô phỏng số cho văn phòng nhà máy Dalat Milk tại Lâm Đồng" đã giải quyết triệt để bài toán tối ưu hóa vi khí hậu và chất lượng không khí trong không gian làm việc công nghiệp. Thông qua việc kết hợp tính toán phụ tải Carrier, mô phỏng số ANSYS Fluent 2022 R1 trên lưới Poly-Hexcore $1.4$ triệu phần tử và phân tích thực nghiệm Taguchi - ANOVA trên Minitab, nghiên cứu đã chứng minh:

  • Cài đặt góc mở cánh đảo gió cassette ở mức $30^\circ$ giúp giảm thiểu chỉ số gió lùa $\text{DR} < 15%$, đảm bảo chỉ số $\text{PMV}$ đạt mức tối ưu từ $-0.38$ đến $+0.05$ và $\text{PPD} < 8.5%$.
  • Cấp gió tươi qua hệ thống PAU với lưu lượng $400\text{ m}^3/\text{h}$ kiểm soát hoàn hảo nồng độ $\text{CO}_2 < 680\text{ ppm}$.
  • Định hướng giải pháp kỹ thuật dài hạn là ứng dụng hệ thống dàn lạnh âm trần nối ống gió nhằm mang lại tiện nghi nhiệt và giá trị thẩm mỹ vượt trội.

Kết quả của đề tài đóng góp nguồn tài liệu kỹ thuật giá trị cho các kỹ sư HVAC, chủ đầu tư nhà máy và sinh viên ngành Công nghệ Kỹ thuật Nhiệt trong việc thiết kế và vận hành các hệ thống điều hòa không khí hiện đại, tiết kiệm năng lượng và hướng tới sức khỏe con người.