ISSN 1859 1531 TẠP CHÍ KHOA HỌC VÀ CÔNG NGHỆ ĐẠI HỌC ĐÀ NẴNG, SỐ 5(102) 2016 137 MÔ PHỎNG SỐ ĐẶC TÍNH CỦA LƯU LƯỢNG KHÍ QUA ỐNG KHÓI NHIỆT LOẠI NGHIÊNG NUMERICAL SIMULATION OF CHARACTERISTICS OF THE I[.]
ISSN 1859-1531 - TẠP CHÍ KHOA HỌC VÀ CƠNG NGHỆ ĐẠI HỌC ĐÀ NẴNG, SỐ 5(102).2016 137 MÔ PHỎNG SỐ ĐẶC TÍNH CỦA LƯU LƯỢNG KHÍ QUA ỐNG KHĨI NHIỆT LOẠI NGHIÊNG NUMERICAL SIMULATION OF CHARACTERISTICS OF THE INDUCED AIR FLOWRATE THROUGH INCLINED SOLAR CHIMNEYS Nguyễn Quốc Ý Trường Đại học Bách khoa Tp HCM; nguyenquocy@hcmut.edu.vn Tóm tắt - Ống khói nhiệt hấp thụ nhiệt xạ mặt trời để tạo hiệu ứng nhiệt dịng khí để thơng gió tự nhiên cho nhà nhà cao tầng Trong nghiên cứu này, khảo sát ảnh hưởng thơng số liên quan kích thước (chiều dài, chiều rộng bề dày kênh dẫn khí), góc nghiêng ống khói nhiệt cường độ nguồn nhiệt phát bề mặt ống khói nhiệt lên đường đặc tính lưu lượng khí theo góc nghiêng Phương pháp mô số CFD (Computational Fluid Dynamics) sử dụng Kết mô trước hết so sánh kiểm tra với kết thực nghiệm Chen cộng (2003) sau dùng để phân tích ảnh hưởng yếu tố Kết cho thấy đường đặc tính lưu lượng phụ thuộc chủ yếu vào tỉ số chiều dài bề dày kênh dẫn khí thơng số khác ảnh hưởng không đáng kể Abstract - Solar chimneys absorb solar radiation heat to induce thermal effects and airflow for natural ventilation of dwellings In this study, we investigate effects of the dimensions of the chimney (length, width and thickness of the air channel), inclined angle and heat flux on the surface of the chimney on the characteristics of the induced air flow rate through the chimney We use CFD (Computational Fluid Dynamics) methods The numerical model is first validated by the experimental results by Chen et al (2003) and then is used for investigating the induced air flow rate when the examined factors are changed The results have shown that the most important factor is the ratio between the length and the thickness of the air channel while other factors have negligible effects on the characteristics of the induced flow rate Từ khóa - thơng gió tự nhiên; ống khói nhiệt; hiệu ứng nhiệt; lưu lượng; CFD Key words - natural ventilation; solar chimney; thermal effect; flow rate; CFD Giới thiệu Ống khói nhiệt hay ống khói mặt trời (solar chimney) thiết bị thơng gió tự nhiên cho nhà hay nhà cao tầng dựa hiệu ứng nhiệt hay hiệu ứng ống khói Ống khói nhiệt thơng thường có cấu tạo gồm kênh dẫn khí với tất mặt bên hấp thụ nhiệt xạ mặt trời Kênh dẫn khí có đầu vào (thấp hơn) nối với khơng gian cần thơng gió đầu (cao hơn) hở khí Ống khói nhiệt nhận nhiệt xạ mặt trời làm tăng nhiệt độ khối khơng khí bên kênh dẫn Khối khơng khí nóng có xu hướng di chuyển lên tạo chênh lệch áp suất để tạo dịng khí cho việc thơng gió tự nhiên Do vậy, lưu lượng thơng gió dựa hiệu ứng nhiệt phụ thuộc hai yếu tố chính: chênh lệch nhiệt độ khối khí chênh lệch độ cao đầu vào đầu ống khói nhiệt Nói cách khác, cường độ xạ mặt trời góc nghiêng ống khói nhiệt ảnh hưởng lớn đến đặc tính làm việc thiết bị Việt Nam, phía Nam, gần xích đạo nên nhận số nắng năm cao với cường độ xạ lớn Các đo đạc thực nghiệm chúng tơi TP Hồ Chí Minh cho thấy cường độ nhiệt xạ mặt trời lên đến 1100W/m2 vào trưa tháng 11 [1] Điều tạo điều kiện thuận lợi cho việc sử dụng ống khói nhiệt để thơng gió tự nhiên cho nhà nhà cao tầng Việt Nam, góp phần tiết kiệm lượng tiêu thụ cho hệ thống thơng gió điều hồ khơng khí khí đảm bảo chất lượng khơng khí nhà Do cường độ xạ mặt trời nhận khác vị trí khác Trái đất, Mathur cộng [2] xây dựng mơ hình tốn cho ống khói nhiệt dựa vào quỹ đạo Mặt trời để tính tốn góc nghiêng tối ưu để ống khói nhiệt nhận lượng nhiệt xạ mặt trời nhiều Họ đề xuất góc nghiêng tối ưu từ 45o đến 60o so với phương ngang cho vĩ độ từ đến 60o Harris Helwig [3] dùng mơ số để dự đốn lưu lượng qua ống khói nhiệt Edinburg (Anh) với vĩ độ 52o Kết họ cho thấy, vào thời điểm 12:00 góc nghiêng 67,5o, tạo lưu lượng lớn nhất, vào thời điểm 16:00 góc nghiêng 45o, tạo lưu lượng lớn Để tập trung khảo sát ảnh hưởng góc nghiêng lên đặc tính làm việc ống khói nhiệt, nhiệt xạ mặt trời mơ nguồn nhiệt phân bố bề mặt ống khói nhiệt Với phương pháp này, Chen cộng [4] Adam cộng [5] nghiên cứu ảnh hưởng cường độ nhiệt góc nghiêng thực nghiệm Thí nghiệm Chen cộng [4] thực ống khói nhiệt nghiêng phịng thí nghiệm Kết họ cho thấy, lưu lượng qua ống khói nhiệt với chiều dài 1,5m chiều rộng kênh dẫn khí 0,2m đạt giá trị lớn góc nghiêng 45o lớn 45% so với trường hợp thẳng đứng Thí nghiệm Adam cộng [5] tiến hành tương tự Chen cộng [4], với ống khói nhiệt có chiều dài 2,0m chiều rộng thay đổi từ 0,1m đến 0,3m Kết Adam cộng [5] cho thấy, ống khói nhiệt với góc nghiêng 45o cho lưu lượng khí lớn Do hai nghiên cứu Chen cộng [4] Adam cộng [5] cho thấy mối liên hệ lưu lượng góc nghiêng cho trường hợp kích thước (chiều dài, chiều rộng bề dày) ống khói nhiệt, chúng tơi muốn tìm hiểu thêm ảnh hưởng kích thước lên lưu lượng góc nghiêng ống khói nhiệt thay đổi Trong nghiên cứu này, dùng mô số để khảo sát ảnh hưởng của: chiều dài, chiều rộng, bề dày, góc nghiêng cường độ nguồn nhiệt lên lưu lượng khí 138 Nguyễn Quốc Ý Phương pháp mơ số CFD Mơ hình ống khói nhiệt dùng để mơ số giống mơ hình thực nghiệm Chen cộng [4] Adam cộng [5] mơ tả Hình Các kích thước ống khói nhiệt bao gồm: chiều dài L, chiều rộng W, bề dày kênh dẫn khí G góc nghiêng so với phương ngang θ Nguồn nhiệt phân bố bề mặt ống khói nhiệt với cường độ I L tử 1,5.10 đủ để mô trường hợp Chúng tơi giữ ngun kích thước lưới tương ứng với trường hợp cho trường hợp mơ khác Điều kiện biên cho khối khơng khí xung quanh ống khói nhiệt áp suất khí Do đó, đầu vào đầu ống khói nhiệt có điều kiện áp suất khí Gia tốc trọng trường tác dụng theo phương đứng hướng xuống Bề mặt kênh dẫn khí bên ống khói nhiệt áp đặt cường độ nhiệt phân bố Các trình truyền nhiệt: xạ từ bề mặt đến bề mặt cịn lại mơi trường, dẫn nhiệt qua mặt bên ống khói nhiệt đối lưu nhiệt bên bên kênh dẫn khí mơ Lớp cách nhiệt bên ống khói nhiệt chọn đủ dày để khơng có thất nhiệt từ bên ống khói nhiệt qua mặt W θ G Hình Mơ hình ống khói nhiệt nghiêng Mơ số thực phương pháp Động lực học lưu chất tính tốn CFD phần mềm Cradle ScStream Hệ phương trình bảo tồn mơ bao gồm: - Phương trình liên tục; - Phương trình động lượng, hay phương trình Navier-Stokes; - Phương trình lượng Các giả thiết tốn bao gồm: - Dịng khí khơng nén được; - Các trình động lực học lưu chất truyền nhiệt ổn định; - Quan hệ khối lượng riêng khơng khí nhiệt độ mơ tả theo quan hệ Boussinesq Mơ hình RANS k-epsilon sử dụng cho mơ hình tính tốn rối Hệ phương trình xấp xỉ phương pháp thể tích hữu hạn Chi tiết phương trình mơ hình tìm thấy tài liệu CFD, ví dụ Ferziger Peric [6] Miền tính tốn, mơ tả Hình 2, bao gồm ống khói nhiệt khối khơng khí bao quanh Bề mặt ống khói nhiệt mơ tả ba loại vật liệu: kính cho mặt hai mặt bên, nhơm cho bề mặt phát nhiệt bên vật liệu cách nhiệt bên bề mặt phát nhiệt Miền tính tốn kiểm tra lựa chọn để không ảnh hưởng đến kết tính với kích thước: bốn mặt bên bên cách ống khói nhiệt 1,0L mặt cách ống khói nhiệt 2,0L Hình cho thấy phân bố nút lưới tính tốn Lưới vng sử dụng cho miền tính tốn với ba kích thước lưới khác Kích thước lưới nhỏ bao quanh ống khói nhiệt 1,0mm Để khảo sát ảnh hưởng kích thước lưới tính tốn, chúng tơi mơ lại thí nghiệm Chen cộng [4] với trường hợp góc nghiêng 75o với nhiều mật độ lưới Hình cho thấy lưới với tổng số phần Hình Miền tính tốn phân bố lưới tốn Hình Khảo sát ảnh hưởng mật độ lưới tính tốn: Lưu lượng khí lưu thơng qua ống khói nhiệt thay đổi theo tổng số phần tử nút lưới Nghiên cứu Bansal cộng [7] cho thấy ảnh hưởng nhiệt độ môi trường lên lưu lượng qua ống khói nhiệt bỏ qua khoảng nhiệt độ từ 20oC đến 40oC Do đó, chúng tơi chọn nhiệt độ mơi trường 30oC tất mô Kết bàn luận Trong phần này, kết mơ so sánh với thí nghiệm Chen cộng [4] để kiểm tra độ xác mơ hình mơ số Sau đó, mơ hình mô số thực với nhiều giá trị của: - Cường độ nhiệt I bề mặt ống khói nhiệt, - Góc nghiêng θ, - Các kích thước W, G, L ISSN 1859-1531 - TẠP CHÍ KHOA HỌC VÀ CÔNG NGHỆ ĐẠI HỌC ĐÀ NẴNG, SỐ 5(102).2016 để đánh giá ảnh hưởng thông số lên lưu lượng dịng khí sinh hiệu ứng nhiệt bên kênh dẫn Lưu lượng tính theo vận tốc trung bình V dịng khí kênh tiết diện mặt cắt ngang A kênh: = × (1) với = × Vận tốc bình V dịng khí cảm ứng kênh tính cơng sụ có sẵn phần hậu xử lý (PostProcessor) phần mềm Cradle ScStream diện tích mặt cắt ngang A (cũng diện tích mặt cắt ướt dịng khí kênh): nghiêng chiều rộng W thay đổi Ống khói nhiệt có hai kích thước L=0,5m G=0,05m giữ cố định W thay đổi với giá trị 0,1m, 0,15m, 0,25m 0,5m (2) Hình So sánh lưu lượng tính tốn lưu lượng thí nghiệm Chen cộng [4] 1.6 1.4 1.2 Q/Q90 đó: dA: diện tích phân tố A, ⃗: vector pháp tuyến dA, ⃗: vector vận tốc dịng khí qua diện tích dA, dấu ⋅ tích vơ hướng hai vector Các diện tích phân tố dA xác định dựa vào phần tử lưới tính tốn nằm diện tích A Do vậy, ⃗ vận tốc dịng khí mơ nút lưới tính tốn Theo cơng thức (2), vận tốc V tính theo thành phần vận tốc vng góc với diện tích dA, tức theo chiều dọc theo kênh Nhờ đó, lưu lượng khí tính theo công thức (1) Để thống cho tất trường hợp tính, chúng tơi chọn mặt cắt ướt để tính V theo cơng thức (2) lối vào (đầu dưới) ống khói nhiệt 3.1 So sánh với kết thí nghiệm Chen cộng [4] Mơ số thực cho trường hợp thí nghiệm Chen cộng [4] Các thông số ống khói nhiệt bao gồm: L=1,5m, G=0,2m, W=0,61m, I=400W/m2, góc θ thay đổi Lưu lượng khí mơ so sánh với lưu lượng khí thực nghiệm Hình Kết mơ mơ tả xu hướng thay đổi lưu lượng góc nghiêng θ thay đổi Mức khác biệt lớn hai kết 7,0% góc nghiêng 60o 3.2 Ảnh hưởng cường độ nhiệt lên lưu lượng Mơ hình mơ số sử dụng để khảo sát ảnh hưởng cường độ nhiệt lên lưu lượng khí Ống khói nhiệt khảo sát có kích thước L=1,0m, W=0,5m, G=0,2m cường độ nhiệt I với ba mức 300W/m2, 600W/m2 900W/m2 Hình cho thấy thay đổi lưu lượng góc nghiêng thay đổi Để tiện cho việc so sánh, lưu lượng khí Q góc nghiêng thể tỉ số Q/Q90, Q90 lưu lượng khí góc nghiêng 90o (tương ứng trường hợp ống khói nhiệt thẳng đứng) Hình cho thấy tỉ số Q/Q90 cho ba trường hợp cường độ nhiệt I gần trùng với lưu lượng cực đại góc nghiêng 45o Điều cho thấy cường độ nhiệt khơng ảnh hưởng nhiều đến đường quan hệ Q θ, trường hợp khảo sát Trong mô tiếp theo, giữ cường độ nhiệt 600W/m2 khảo sát ảnh hưởng thông số khác 3.3 Ảnh hưởng chiều rộng W lên lưu lượng Hình cho thấy mối quan hệ lưu lượng góc 0.8 300W/m2 0.6 600W/m2 0.4 900W/m2 0.2 0 20 40 60 Góc nghiêng (Độ) 80 100 Hình Ảnh hưởng cường độ nhiệt lên lưu lượng khí cho trường hợp L=1,0m, W=0,5m, G=0,2m 1.8 1.6 1.4 1.2 Q/Q90 = ∫ ⃗ ⋅ ⃗ 139 0.8 W=0,5m 0.6 W=0,25m 0.4 W=0,15m 0.2 W=0,1m 0 20 40 60 Góc nghiêng (Độ) 80 100 Hình Ảnh hưởng chiều rộng W lên lưu lượng khí cho trường hợp L=0,5m, G=0,05m I=600W/m2 Kết Hình cho thấy dạng đường quan hệ Q θ tương tự W thay đổi, tỉ số Q/Q90 có xu hướng giảm W tăng góc nghiêng θ Lưu lượng cực đại trường hợp đạt góc nghiêng 60o 3.4 Ảnh hưởng bề dày G chiều dài L lên lưu lượng Để khảo sát ảnh hưởng bề dày G, chiều dài L giữ cố định 1,0m để khảo sát ảnh hưởng chiều dài L, bề dày G giữ cố định 0,05m Trong hai trường hợp, W I 0,5m 600W/m2 Để tiện cho việc so sánh, dùng tỉ số L/G để 140 Nguyễn Quốc Ý khảo sát Kết thể Hình Hình Xu hướng thay đổi tỉ số Q/Q90 hai Hình tương tự Thứ nhất, góc nghiêng cho lưu lượng cực đại tăng dần theo tỉ số L/G, tương ứng 45o, 60o, 75o tỉ số L/G tăng lên từ 5, 10, 20 Thứ hai, tỉ số L/G lớn 20, lưu lượng có xu hướng tăng dần theo góc nghiêng ống khói nhiệt, nên lưu lượng đạt cực đại θ =90o hay ống khói nhiệt thẳng đứng 1.6 1.4 o khí so với góc nghiêng 45 nên lưu lượng 30o nhỏ lưu lượng 45o khoảng 5,0% Khi góc nghiêng tăng lên, thể tích lớp khơng khí nhận nhiệt thu hẹp dần sát bề mặt ống khói nhiệt Do đó, lớp khơng khí có vận tốc lớn mỏng dần nên làm giảm lưu lượng Ở góc nghiêng nhỏ hơn, vùng khơng khí nhận nhiệt từ bề mặt rộng trường hợp 45o, chênh lệch độ cao đầu vào đầu ống khói nhiệt nhỏ hơn, nên tạo hiệu ứng nhiệt lưu lượng khí tương ứng 1.2 Q/Q90 0.8 L/G=20 0.6 L/G=5 L/G=10 0.4 L/G=14,3 0.2 0 20 40 60 Góc nghiêng (Độ) 80 100 Hình Ảnh hưởng bề dày G lên lưu lượng khí với L=1,0m, W=0,5m, I=600W/m2 1.6 1.4 1.2 Q/Q90 0.8 L/G=20 0.6 L/G=5 0.4 L/G=10 0.2 L/G=25 0 20 40 60 Góc nghiêng (Độ) 80 100 Hình Ảnh hưởng chiều dài L lên lưu lượng khí với G=0,05m, W=0,5m, I=600W/m2 Như vậy, thông số ảnh hưởng nhiều đến đường đặc tính lưu lượng dịng khí cảm ứng qua ống khói nhiệt tỉ số L/G chiều dài bề dày kênh dẫn khí Để thấy rõ ảnh hưởng tỉ số lên dịng khí cảm ứng q trình truyền nhiệt bên ống khói nhiệt, khảo sát phân bố nhiệt độ vận tốc cho hai trường hợp có tỉ số L/G nhỏ LG/ lớn Hình cho thấy phân bố nhiệt độ vận tốc dịng khí trường hợp có tỉ số L/G với ba góc nghiêng khác Các hình thể với thang màu Để tiện cho việc so sánh, vùng có nhiệt độ lớn 1oC so với nhiệt độ môi trường thể màu đỏ Phân bố nhiệt độ cho thấy nhiệt độ khơng khí bên ống khói nhiệt cao trường hợp 30o Tuy nhiên, lưu lượng cực đại đạt 45o, thể Hình 5, Điều góc nghiêng 30o, vùng tách dịng (vùng khoanh trịn Hình 9) đầu vào đầu ống khói nhiệt làm giảm tiết diện hiệu dụng A dịng Hình Phân bố nhiệt độ (bên trái) vận tốc dịng khí (bên phải) cho trường hợp L/G=5 (G=0,2m, W=0,5m, L=1,0m, I=600W/m)2 với ba góc nghiêng 30o, 45o 75o (từ xuống) Hình 10 Phân bố nhiệt độ cho trường hợp L/G 25 (G=0,05m, W=0,5m, L=1,25m, I=600W/m2) với ba góc nghiêng 45o, 75o 90o (từ trái qua phải) ISSN 1859-1531 - TẠP CHÍ KHOA HỌC VÀ CƠNG NGHỆ ĐẠI HỌC ĐÀ NẴNG, SỐ 5(102).2016 Hình 10 cho thấy phân bố nhiệt độ trường hợp có tỉ số L/G 25 Vùng màu đỏ thể lớp không khí có nhiệt độ cao 1oC so với mơi trường Phân bố nhiệt độ cho thấy tỉ số L/G này, thất nhiệt từ bên ống khói nhiệt vùng khơng khí phía kính (bên ngồi ống khói nhiệt) tăng lên góc nghiêng nhỏ Với công suất nguồn nhiệt mặt dưới, lượng nhiệt thất thoát nhiều dẫn đến lượng nhiệt mà khơng khí bên kênh dẫn nhận giảm làm giảm hiệu ứng nhiệt Do lưu lượng khí giảm góc nghiêng giảm đạt cực đại góc nghiêng lớn nhất, hay ống khói nhiệt thẳng đứng, thể Hình 3.5 Về khả áp dụng ống khói nhiệt để thơng gió tự nhiên cho cơng trình Việt Nam Với cường độ xạ mặt trời cao đến 1100W/m2 Tp Hồ Chí Minh [1], ống khói nhiệt giải pháp giúp thơng gió tự nhiên hiệu cho cơng trình địa phương để góp phần tiết kiệm lượng đảm bảo chất lượng khơng khí bên cơng trình Ví dụ, lưu lượng mà ống khói nhiệt Hình với kích thước (chiều dài x chiều rộng x bề dày) 1,0m x 0,5m x 0,2m (tương đương diện tích nửa cửa sổ thơng thường) tạo cường độ nhiệt 900W/m2 góc nghiêng 45o 27 lít/giây hay 97,2 m3/giờ Lưu lượng đủ để thơng gió tự nhiên cho phịng có kích thước 3m x 3m x 3m với bội số tuần hoàn 3,6 đáp ứng u cầu thơng gió theo tiêu chuẩn quốc gia TCVN 5687:2000 [8] Mathur cộng [2] đề xuất góc nghiêng ống khói nhiệt để nhận nhiệt xạ mặt trời nhiều vĩ độ Tp Hồ Chí Minh 50o Tuy nhiên, theo kết chúng tơi Hình 8, góc nghiêng hiệu tỉ số L/G nằm khoảng từ - 10 Do vậy, để tối ưu lưu lượng khí tạo được, ống khói nhiệt cho cơng trình Tp Hồ Chí Minh nên thiết kế với tỉ số L/G nằm khoảng giá trị với góc nghiêng lân cận 50o so với phương ngang 141 Kết luận Chúng sử dụng mô số để khảo sát đặc tính lưu lượng khí lưu thơng qua ống khói nhiệt loại nghiêng Kết mơ cho thấy, thơng số ảnh hưởng nhiều lên đường đặc tính lưu lượng tỉ số chiều dài L bề dày kênh dẫn khí G Cụ thể, thơng số khác nhau, góc nghiêng ống khói nhiệt ứng với mức lưu lượng cực đại tăng dần theo tỉ số L/G Khi tỉ số L/G lớn 20, lưu lượng đạt cực đại ống khói nhiệt thẳng đứng Các thơng số khác bao gồm chiều rộng W kênh dẫn khí cường độ nhiệt khơng ảnh hưởng đến đường đặc tính lưu lượng theo góc nghiêng Lời cảm ơn: Nghiên cứu tài trợ Trường Đại học Bách khoa - ĐHGQ Tp HCM mã số T-KTXD2016-86 Công ty Cradle CFD Acknowledgement: The author thanks Cradle CFD Co Ltd for providing a ScStream licence for this study TÀI LIỆU THAM KHẢO [1] Nguyễn Quốc Ý, “Thử nghiệm giải pháp làm mát khơng khí tiết kiệm lượng: Ống khói nhiệt kết hợp hiệu ứng bay nước – Làm mát gián tiếp điều kiện thời tiết thực Tp Hồ Chí Minh”, Tạp chí KHCN Đại học Đà Nẵng số (94), 2015, 69-73 [2] J Mathur, S Mathur, Anupma, “Summer-performance of inclined roof solar chimney for natural ventilation”, Energy and Buildings, Elsevier, 38, 2006, 1156-1163 [3] D.J Harris N Helwig, “Solar chimney and building ventilation”, Applied Energy, Elsevier, 84, 2007, 135-146 [4] Z.D Chen, P Bandopadhayay, J Halldorsson, C Byrjalen, P Heiselberg, Y Li, “An experimental investagation of a solar chimney model with uniform wall heat flux”, Building and Environment, Elsevier, 38, 2003, 893-906 [5] J Adam, T Yamanaka, H Kotani., “Mathematical model and experimental study of airflow in solar chimneys”, Proceedings of 8th International Conference on Air Distribution in Rooms (ROOMVENT 2002), 2002, 621-624 [6] J.H Ferziger M Peric, Computational methods for fluid dynamics, Springer, 2002 [7] N.K Bansal, R Mathur, M.S Bhandari, “Solar chimney for enhanced stack ventilation”, Building and Environment, Elsevier, 28 (3), 1993, 373-377 [8] Bộ KHCN Việt Nam, Tiêu chuẩn quốc gia TCVN 5687:2010 Thông gió – Điều hồ khơng khí – Tiêu chuẩn thiết kế, 2010 (BBT nhận bài: 10/03/2016, phản biện xong: 07/04/2016) ... [5] J Adam, T Yamanaka, H Kotani., “Mathematical model and experimental study of airflow in solar chimneys”, Proceedings of 8th International Conference on Air Distribution in Rooms (ROOMVENT 2002),... Chen, P Bandopadhayay, J Halldorsson, C Byrjalen, P Heiselberg, Y Li, “An experimental investagation of a solar chimney model with uniform wall heat flux”, Building and Environment, Elsevier,... 6 21- 624 [6] J.H Ferziger M Peric, Computational methods for fluid dynamics, Springer, 2002 [7] N. K Bansal, R Mathur, M.S Bhandari, “Solar chimney for enhanced stack ventilation”, Building and