Ứng dụng mô hình luồng khí vùng không gian giới hạn trong công trình kiến trúc và an toàn sức khỏe

ỦY BAN NHÂN DÂN TP.HCM ðẠI HỌC QUỐC GIA HỒ CHÍ MINH SỞ KHOA HỌC VÀ CƠNG NGHỆ TRƯỜNG ðẠI HỌC KHOA HỌC TỰ NHIÊN BÁO CÁO NGHIỆM THU (ðã chỉnh sửa theo góp ý Hội đồng nghiệm thu ngày 30 tháng 10 năm 2010 ) ỨNG DỤNG MÔ HÌNH LUỒNG KHÍ VÙNG KHƠNG GIAN GIỚI HẠN TRONG CƠNG TRÌNH KIẾN TRÚC VÀ AN TỒN SỨC KHOẺ CHỦ NHIỆM ðỀ TÀI: PGS TS NGUYỄN THỊ BẢY PGS TS LÊ QUANG TOẠI THÀNH PHỐ HỒ CHÍ MINH THÁNG 12/ 2010 ỦY BAN NHÂN DÂN TP.HCM ðẠI HỌC QUỐC GIA TP HỒ CHÍ MINH SỞ KHOA HỌC VÀ CƠNG NGHỆ TRƯỜNG ðẠI HỌC KHOA HỌC TỰ NHIÊN BÁO CÁO NGHIỆM THU (ðã chỉnh sửa theo góp ý Hội đồng nghiệm thu ngày 30 tháng 10 năm 2010) ỨNG DỤNG MÔ HÌNH LUỒNG KHÍ VÙNG KHƠNG GIAN GIỚI HẠN TRONG CƠNG TRÌNH KIẾN TRÚC VÀ AN TỒN SỨC KHOẺ Cơ quan chủ trì đề tài Chủ nhiệm đề tài PGS TS Nguyễn Thị Bảy PGS.TS Lê Quang Toại Cơ quan quản lý THÀNH PHỐ HỒ CHÍ MINH THÁNG 12/ 2010 Ứng dụng mơ hình luồng khí vùng khơng gian giới hạn cơng trình kiến trúc an tồn sức khoẻ TĨ M TẮ T Mơi trường khơng khí nhà nơi thường bị ô nhiễm nhiều môi trường bên ngồi Ơ nhiễm khơng khí nhà tác động trực tiếp ñến sức khoẻ người Vấn ñề ñã ñược quan tâm nhiều giới Việt Nam năm gần Tuy nhiên Các cơng trình có chưa quan tâm mức đến thay ñổi luồng khí nhà sao, thay ñổi thiết kế ứng với kịch khác luồng khí phân bố chkjkhất nhiễm nào, kịch tốt nhất,… để từ chọn phương án tối ưu thơng khí thiết kế cho cơng trình Một vấn đề cơng trình chưa đầu tư sâu vấn đề nhiễm khơng khí nhà, building,… Mục tiêu đề tài ứng dụng mơ hình luồng khí khơng gian giới hạn tính tốn phân bố luồng khơng khí theo khơng gian (áp suất luồng khí, thay đổi nhiệt độ, biến đổi chất ô nhiễm, ) phục vụ cho việc thiết kế tối ưu cơng trình dân dụng Nhằm thực mục tiêu trên, ñề tài thực nội dung sau: − Giới thiệu tổng quan mơ hình COMTAM phần mềm ANSYS − ðiều tra khảo sát ñối tượng nghiên cứu − Thu thập số liệu đầu vào cho mơ hình − ứng dụng mơ hình luồng khí tính cho đối tượng − ðo đạc thực tế phục vụ tính tốn hiệu chỉnh mơ hình − ðánh giá khả ứng dụng mơ hình ñề xuất phương án ñã ñược hiệu chỉnh − Thể kết tính thiết kế tối ưu thơng khí Autocad ðây đề tài khơng có ý nghĩa khoa học mà cịn có tính ứng dụng cao, thực dựa mơ hình tính tốn với số liệu đo đạc khảo sát thực tế, đáng tin cậy Chúng tơi xin cảm ơn cộng tác ñồng nghiệp trình thực đề tài i Ứng dụng mơ hình luồng khí vùng khơng gian giới hạn cơng trình kiến trúc an tồn sức khoẻ SUMMARY Inside building is usually more polluted than outside building Inside building air pollution effect to human health directly This problem has been much interested in the world as well as in Vietnam several recent years However, most of research works in Vietnam haven’t cared about how the air stream change inside building or air stream and how pollutant distribution will change when building design is changed or which scenario is the best,….Those are very important to chose the optimal ventilation project for designing building Another problem is that the research works in Vietnam haven’t cared about inside building pollution deeply The main goal of this research is applying limited space air stream model to calculate air stream distribution (air stream pressure, temperature change, pollutant change, ) These are helpful for designing civil architectural project optimally To achieve that goal, this topic was conducted with the following contents: - Introduce COMTAM model and ANSYS software generally - Conducting the survey research objects - Collecting input data for calculating by the model - Caculating for research objects by using the model - Surveying some places to serve calculating and correcting model parameters - Accessing the apply capacity of the model and proposing the optimal project - Showing the optimal ventilation result and design in autocad This topic has both scientific signification and high applicability It is conducted base on the model that is tested with theoretical and real result and the data is reliable, it was collected directly We express our deep gratitude to those who actively collaborated for this topic ii Ứng dụng mơ hình luồng khí vùng khơng gian giới hạn cơng trình kiến trúc an tồn sức khoẻ PHẦN MỞ ðẦU Tên đề tài: Ứng dụng mơ hình luồng khí vùng khơng gian giới hạn cơng trình kiến trúc an tồn sức khỏe Chủ nhiệm ñề tài: PGS.TS Nguyễn Thị Bảy PGS.TS Lê Quang Toại Cơ quan chủ trì: Trường ðại học Khoa học Tự nhiên Tp Hồ Chí Minh Thời gian thực đề tài: Từ tháng 12/2006 đến tháng 12/2007 Kinh phí duyệt: 255.000.000đ Mục tiêu: Ứng dụng mơ hình luồng khí khơng gian giới hạn tính tốn phân bố luồng khơng khí theo khơng gian (Áp suất luồng khí, thay ñổi nhiệt ñộ, biến ñổi chất ô nhiễm,… phục vụ cho việc thiết kế tối ưu cơng trình dân dụng) Nội dung: - Giới thiệu tổng quan mơ hình COMIS - ðiều tra khảo sát đối tượng nghiên cứu: Một ngơi nhà đơn vùng, Một bệnh viện, Một building - Thu thập số liệu đầu vào mơ hình - Ứng dụng mơ hình luồng khí COMIS tính tốn cho đối lượng ưu tiên nói - ðo đạc thực tế phục vụ tính tốn hiệu chỉnh mơ hình - ðánh giá khả ứng dụng mơ hình đề xuất phương án ñã ñược hiệu chỉnh - Thể kết tính thiết kế tối ưu khơng khí Autocad - Chuyển giao công nghệ iii Ứng dụng mơ hình luồng khí vùng khơng gian giới hạn cơng trình kiến trúc an tồn sức khoẻ MỤC LỤC TÓM TẮT i SUMMARY ii MỤC LỤC iii DANH MỤC HÌNH vii DANH MỤC BẢNG xi MỞ ðẦU CHƯƠNG I: TỔNG QUAN VỀ Ô NHIỄM KHÔNG KHÍ VÀ THƠNG GIĨ TRONG KHƠNG GIAN GIỚI HẠN 1.1 NGUỒN Ơ NHIỄM KHƠNG KHÍ TRONG NHÀ 1.1.1 Nguồn ô nhiễm bên 1.1.2 Nguồn ô nhiễm từ bên 1.2.CÁC CHẤT Ơ NHIỄM – GÂY BỆNH TRONG KHƠNG KHÍ Ở TRONG NHÀ, VĂN PHÒNG 1.2.1 Radon 1.2.2 Formaldehyde 1.2.3 Các chất hữu bay (VOCs) iv Ứng dụng mô hình luồng khí vùng khơng gian giới hạn cơng trình kiến trúc an tồn sức khoẻ 1.2.4 Các chất sát trùng gia dụng 1.2.5 Các chất khí độc từ nhà bếp 10 1.2.6 Các chất gây ô nhiễm sinh học nhà .11 1.3 TÌNH HÌNH Ơ NHIỄM KHƠNG KHÍ TRONG NHÀ 12 1.3.1 Ô nhiễm nhà so với nhiễm ngồi trời 12 1.3.2 Một số kết nghiên cứu tương quan nhiễm ngồi trời nhiễm nhà cơng bố 14 1.4 TÌNH TRẠNG Ô NHIỄM KHÔNG KHÍ TRONG NHÀ Ở CÁC NƯỚC ðANG PHÁT TRIỂN 17 1.4.1 Tình hình chung .17 1.4.2 Ơ nhiễm khơng khí nhà Việt Nam 19 1.5 SỰ THƠNG KHÍ TRONG NHÀ 20 1.5 CHẤT LƯỢNG KHƠNG KHÍ TRONG NHÀ ẢNH HƯỞNG ðẾN SỨC KHỎE CON NGƯỜI 20 1.5.1 Nhận biết khơng khí nhà gây nguy hiểm cho sức khoẻ 20 1.5.2 Sự phân bố khơng khí nhà phịng 22 1.5.3 Ơ nhiễm khơng khí nhà gây tác động xấu đến sức khỏe 22 1.6 CÁC BIỆN PHÁP NHẰM GIẢM THIỂU Ơ NHIỄM KHƠNG KHÍ TRONG NHÀ – VĂN PHỊNG .33 1.6.1 Các biện pháp chung 33 1.6.2 Chi tiết số biện pháp chất nhiễm – gây bệnh 35 1.6.3 Mười cách ñể giảm thiểu ô nhiễm nhà .38 1.7 LUỒNG KHÍ TRONG NHÀ LIÊN QUAN ðẾN SỨC KHỎE CON NGƯỜI 39 CHƯƠNG II: GIỚI THIỆU CÁC PHẦN MỀM TÍNH TỐN Ơ NHIỄM VÀ THƠNG KHÍ TRONG NHÀ 41 2.1 GIỚI THIỆU VỀ CHƯƠNG TRÌNH CONTAM 41 2.1.1 Giới thiệu 41 2.1.2 Một số ứng dụng mơ hình CONTAM .42 2.1.3 Các giả định mơ hình CONTAM 44 2.1.4 Ví dụ mẫu bước sử dụng chương trình CONTAM .46 2.2 GIỚI THIỆU VỀ PHẦN MỀM ANSYS 55 2.2.1 Cơ sở lí thuyết Flotran CFD 56 v Ứng dụng mơ hình luồng khí vùng khơng gian giới hạn cơng trình kiến trúc an toàn sức khoẻ 2.2.2 Lời giải số Flotran CFD .56 2.2.3 Kiểm tra mơ hình Ansys 57 2.2.4 Kết tính tốn phân tích 58 CHƯƠNG III: ỨNG DỤNG TÍNH TỐN THƠNG KHÍ VÀ Ơ NHIỄM KHƠNG KHÍ TRONG NHÀ .63 3.1 TÍNH CHO MẪU NHÀ VILA 63 3.1.1 Mô tả khu vực nhà Villa .63 3.1.2 Kết ño ñạc khí tượng quanh ngơi nhà .64 3.1.3 Thông số nguồn thải .68 3.1.4 Các kịch tính tốn thơng khí 69 3.1.5 Các kịch tính tốn nhiễm nhà Villa 75 3.2: TÍNH CHO MẪU NHÀ HÌNH HỘP 88 3.2.1 Mơ tả khu vực ngơi nhà hình hộp 88 3.2.2 ðiều kiện khí tượng quanh ngơi nhà .90 3.2.3 Các kịch tính tốn thơng khí 90 3.2.4 Các kịch tính tốn nhiễm nhà hộp .94 3.3: TÍNH CHO MẪU NHÀ BUILDING 107 3.3.1 Mô tả khu vực nhà Building 107 3.3.2 Thơng số nguồn thải tồ nhà 107 3.3.3 Tính tốn thơng khí 108 3.3.4 Các kịch tính tốn nhiễm ñối với nhà Building 112 CHƯƠNG IV: KẾT LUẬN & KIẾN NGHỊ 132 Kết luận 132 Kiến nghị 133 PHỤ LỤC HÌNH 18235 TÀI LIỆU THAM KHẢO 182 vi Ứng dụng mơ hình luồng khí vùng khơng gian giới hạn cơng trình kiến trúc an tồn sức khoẻ DANH MỤC HÌNH Hình 1.1: Các nguồn gây ô nhiễm không khí nhà .13 Hình 2.1: Giao diện người sử dụng 47 Hình 2.2: Khai báo thông số thời tiết 47 Hình 2.3: Khai báo thơng số gió 48 Hình 2.4: Khai báo thơng số vị trí tọa độ, vùng độ cao so với mực nước biển 48 Hình 2.5: Khai báo hiển thị thơng số áp suất gió 48 Hình 2.7: Khai báo thành phần luồng khí .49 Hình 2.8: Khai báo thơng số lọc lịch thời gian 50 Hình 2.9: Khai báo độ cao, hướng quy định giới hạn giá trị áp suất, luồng khí 50 Hình 2.10: Khai báo áp suất gió 51 Hình 2.12: Khai báo mơ hình thành phần nguồn 52 Hình 2.13: Khai báo thơng số cho mơ hình thành phần nguồn .52 Hình 2.14: Khai báo thơng số chạy mơ .53 Hình 2.15: Kiểm tra thông số thời tiết 53 Hình 2.16: Kết xuất tính tốn xong 54 Hình 2.17: Các phương pháp tính tham số cho luồng khí 54 Hình 2.18: Các phương pháp tính hệ số cho chất ô nhiễm .54 Hình 2.20 Lưới tính tốn dịng chảy bao quanh trụ Hình 2.21 Vận tốc biên mặt trụ Hình 2.23 Trường vận tốc dịng chảy bao quanh trụ tính từ Ansys, chảy tầng, Re=1,25; u0=0,0001m/s 59 Hình 2.22 Giá trị vận tốc đường trịn r = cm dịng chảy tầng, Re=1,25; u0=0,0001m/s ( a- tính từ lí thuyết; b- tính từ Ansys) Hình 2.24 Giá trị vận tốc ñường tròn r = cm dòng chảy tầng, Re=1,25; u0=0,0001m/s Hình 2.25 Trường vận tốc dòng bao quanh trụ Từ Ansys, mơ hình rối k-ε, Re=19; u0=0,015m/s Hình 2.26 Trường vận tốc dòng bao quanh trụ.Từ Ansys, mơ hình rối k-ε, Re=125; u0=0,1m/s Hình 2.27 Trường vận tốc dịng bao quanh trụ Từ Ansys, mơ hình rối k-ε, Re=188; u0=0,15m/s vii Ứng dụng mơ hình luồng khí vùng khơng gian giới hạn cơng trình kiến trúc an tồn sức khoẻ Hình 2.28 Trường vận tốc dòng bao quanh trụ.Từ Ansys, mơ hình rối k-ε, Re=250; u0=0,2m/s Hình 3.1: Mô tả tầng nhà CONTAM 63 Hình 3.2: Mơ tả tầng nhà CONTAM 64 Hình 3.3: Lưới tính tốn cho mẫu nhà vila trạng (tầng trệt) Hình 3.4: Lưới tính tốn cho mẫu nhà vila trạng (lầu 1) Hình 3.5 Áp suất khu vực vào buổi trưa ngày thứ 76 Hình 3.6 Áp suất khu vực vào buổi chiều ngày thứ 77 Hình 3.7 Thang phân chia áp suất Hình 3.8 Nhiệt ñộ khu vực vào buổi trưa ngày thứ 78 Hình 3.9 Nhiệt ñộ khu vực vào buổi chiều ngày thứ .78 Hình 3.10.Thang chia nhiệt ñộ .1 Hình 3.12 Nồng độ CO khu vực vào buổi chiều ngày thứ hai .80 Hình 3.14 Nồng độ SO2 khu vực vào buổi trưa ngày thứ hai 80 Hình 3.13 Thang chia nồng độ CO .1 Hình 3.15 Nồng ñộ SO2 khu vực vào buổi chiều ngày thứ hai 81 Hình 3.17 Nồng ñộ NO2 cc khu vực vào buổi trưa ngày thứ hai .81 Hình 3.16.Thang chia nồng ñộ SO2 Hình 3.43 Thang nồng độ NO2 .1 Hình 3.44 Diễn biến nồng ñộ CO 48 tầng 83 Hình 3.45 Diễn biến nồng ñộ CO 48 tầng 83 Hình 3.46 Diễn biến nồng ñộ NO2 48 tầng 84 Hình 3.47 Diễn biến nồng ñộ NO2 48 tầng 84 Hình 3.48 Diễn biến nồng ñộ NO2 48 nhà bếp 86 Hình 3.49 Diễn biến nồng độ NO2 48 nhà bếp 88 Hình 3.50: Mơ tả tầng nhà CONTAM 89 Hình 3.51: Mơ tả tầng ngơi nhà CONTAM .89 Hình 3.52: Lưới tính tốn cho mẫu nhà hộp trạng (tầng trệt) Hình 3.53: Lưới tính tốn cho mẫu nhà hộp trạng (lầu 1) .1 Hình 3.66 Áp suất khu vực vào buổi trưa ngày thứ .95 Hình 3.67 Áp suất khu vực vào buổi chiều ngày thứ .95 Hình 3.68 Thang phân chia áp suất .97 viii Ứng dụng mơ hình luồng khí vùng khơng gian giới hạn cơng trình kiến trúc an tồn sức khoẻ PHỤ LỤC HÌNH Hình 3.92a: Phân bố vận tốc tầng (kịch trạng), vận tốc gió SE = 0,5m/s Hình 3.92b: Phân bố vận tốc tầng (k.bản trạng), vt gió SE = 0,5m/s, phóng to cửa vào & nhà giữ xe Hình 3.92c: Phân bố vận tốc tầng (k trạng), vt gió SE = 0,5m/s, phóng to nhà giữ xe & hành lang 170 Ứng dụng mơ hình luồng khí vùng khơng gian giới hạn cơng trình kiến trúc an tồn sức khoẻ Hình 3.93a: Phân bố vận tốc tầng (kịch trạng), vận tốc gió SE = 1m/s Hình 3.93b: Phân bố vận tốc tầng (k trạng), vt gió SE = 1m/s, phóng to c? a vào & nhà giữ xe Hình 3.93c: Phân bố vận tốc tầng (k trạng), vt gió SE = 1m/s, phóng to nhà giữ xe & hành lang 171 Ứng dụng mơ hình luồng khí vùng khơng gian giới hạn cơng trình kiến trúc an tồn sức khoẻ Hình 3.94a: Phân bố vận tốc tầng (kịch trạng), vận tốc gió SE = 1,5m/s Hình 3.94b: Phân bố vận tốc tầng (k.bản trạng), vt gió SE = 1,5m/s, phóng to cử a vào & nhà giữ xe Hình 3.94c: Phân bố vận tốc tầng (k trạng), vt gió SE = 1,5m/s, phóng to nhà giữ xe & hành lang 172 Ứng dụng mơ hình luồng khí vùng khơng gian giới hạn cơng trình kiến trúc an tồn sức khoẻ Hình 3.95a: Phân bố vận tốc tầng (kịch trạng), vận tốc gió SW = 0,5m/s Hình 3.95b: Phân bố vận tốc tầng (k.bản trạng), vt gió SW = 0,5m/s, phóng to cửa vào & nhà giữ xe thang Hình 3.95c: Phân bố vận tốc tầng (k trạng), vt gió SW = 0,5m/s, phóng to nhà giữ xe & hành lang 173 Ứng dụng mơ hình luồng khí vùng khơng gian giới hạn cơng trình kiến trúc an tồn sức khoẻ Hình 3.96a: Phân bố vận tốc tầng (kịch trạng), vận tốc gió SW = 1m/s Hình 3.96b: Phân bố vận tốc tầng (k.bản trạng), vt gió SW = 1m/s, phóng to cửa vào & nhà giữ xe Hình 3.96c: Phân bố vận tốc tầng (k trạng), vt gió SW = 1m/s, phóng to nhà giữ xe & hành lang 174 Ứng dụng mơ hình luồng khí vùng khơng gian giới hạn cơng trình kiến trúc an tồn sức khoẻ Hình 3.97a: Phân bố vận tốc tầng (kịch trạng), vận tốc gió SW = 1,5m/s Hình 3.97b: Phân bố vận tốc tầng (k.bản trạng), vt gió SW = 1,5m/s, phóng to cửa vào & nhà giữ xe Hình 3.97c: Phân bố vận tốc tầng (k trạng), vt gió SW = 1,5m/s, phóng to nhà giữ xe & hành lang 175 Ứng dụng mơ hình luồng khí vùng khơng gian giới hạn cơng trình kiến trúc an tồn sức khoẻ Hình 3.98a: Phân bố vận tốc lầu 1-6 (kịch trạng), vận tốc gió SE = 0,5m/s Hình 3.98b: Phân bố vận tốc lầu 1-6 (k.bản trạng), vt gió SE = 0,5m/s, phóng to phịng 4&5 Hình 3.98c: Phân bố vận tốc lầu 1-6 (k trạng), vt gió SE = 0,5m/s, phóng to vị trí sân phơi đồ 176 Ứng dụng mơ hình luồng khí vùng khơng gian giới hạn cơng trình kiến trúc an tồn sức khoẻ Hình 3.99a: Phân bố vận tốc lầu 1-6 (kịch trạng), vận tốc gió SE =1m/s Hình 3.99b: Phân bố vận tốc lầu 1-6 (k.bản trạng), vt gió SE = 1m/s, phóng to hộ 4&5 Hình 3.99c: Phân bố vận tốc lầu 1-6 (k trạng), vt gió SE = 1m/s, phóng to hộ 177 Ứng dụng mơ hình luồng khí vùng khơng gian giới hạn cơng trình kiến trúc an tồn sức khoẻ Hình 3.100a: Phân bố vận tốc lầu 1-6 (kịch trạng), vận tốc gió SE = 1,5m/s Hình 3.100b: Phân bố vận tốc lầu 1-6 (k.bản trạng), vt gió SE = 1,5m/s, phóng to hộ 4&5 Hình 3.100c: Phân bố vận tốc lầu 1-6 (k trạng), vt gió SE = 1,5m/s, phóng to hộ 178 Ứng dụng mơ hình luồng khí vùng khơng gian giới hạn cơng trình kiến trúc an tồn sức khoẻ Hình 3.101a: Phân bố vận tốc lầu 1-6 (kịch trạng), vận tốc gió SW =0,5m/s Hình 3.101b: Phân bố vận tốc lầu 1-6 (k.bản trạng), vt gió SW =0,5m/s, phóng to hộ 5&6 Hình 3.101c: Phân bố vận tốc lầu 1-6 (k trạng), vt gió SW =0,5m/s, phóng to vị trí sân phơi đồ 179 Ứng dụng mơ hình luồng khí vùng khơng gian giới hạn cơng trình kiến trúc an tồn sức khoẻ Hình 3.102a: Phân bố vận tốc lầu 1-6 (kịch trạng), vận tốc gió SW =1m/s Hình 3.102b: Phân bố vận tốc lầu 1-6 (k.bản trạng), vt gió SW =1m/s, phóng to hộ 5&6 Hình 3.102c: Phân bố vận tốc lầu 1-6 (k trạng), vt gió SW =1m/s, phóng to hộ 1&3 180 Ứng dụng mơ hình luồng khí vùng khơng gian giới hạn cơng trình kiến trúc an tồn sức khoẻ Hình 3.103a: Phân bố vận tốc lầu 1-6 (kịch trạng), vận tốc gió SW =1,5m/s Hình 3.103b: Phân bố vận tốc lầu 1-6 (k.bản trạng), vt gió SW =1,5m/s, phóng to hộ 5&6 Hình 3.103c: Phân bố vận tốc lầu 1-6 (k trạng), vt gió SW =1,5m/s, phóng to hộ 1&3 181 Ứng dụng mơ hình luồng khí vùng khơng gian giới hạn cơng trình kiến trúc an toàn sức khoẻ TÀI LIỆU THAM KHẢO AKER, A.J., and KELSO, M (1990) “On validation of computational fluid dynamics procedures for room air motion prediction”, ASHRAE Transactions, 96(1), 760-774 BAUGHMAN, A.V., GADGIL, A.J., and NAZAROFF, W.W (1994) “Mixing of a point source pollutant by natural convection flow within a room”, Indoor Air, 4, 114-122 DRESCHER, A.C (1994) “Computed Tomography and Optical Remote Sensing” Development for the Study of Indoor Air Pollutant Transport and Dispersion, Ph.D Thesis, Dept of Civil and Environmental Engineering, University of California, Berkeley DRESCHER, A.C., LOBASCIO, C., GADGIL, A.J., et al (1995) “Mixing of a point source indoor pollutant by forced convection”, Indoor Air, 5, 204-214 FURTAW, J., PANDIAN, M.D., NELSON, D.R., et al (1996) “Modeling indoor air concentrations near emission sources in imperfectly mixed rooms”, Journal of Air and Waste Management Association, 46, 861-868 GADGIL, A.J., FINLAYSON, E.U., FISCHER, M.L., et al (2000) “Pollutant transport and dispersion in large indoor spaces: a status report for the large space effort of the interiors project”, Lawrence Berkeley National Laboratory Report, LBNL-44791, Berkeley, USA HEISELBERG, P (1996) “Room air and contaminant distribution in mixing ventilation”, ASHRAE Transactions: Symposia, 332-339 ISHIZU, Y (1980) “A general equation for the estimation of indoor pollution”, Environmental Science & Technology, 14, 1254-1257 LAMBERT, W.E., SAMET, J.M and SPENGLER, J.D (1993) “Environmental tobacco smoke concentrations in no-smoking sections of restaurants”, American Journal of Public Health, 83, 1339-1341 10 LOBSCHEID, C (2001) Isothermal mixing of a point source indoor pollutant: numerical predictions and comparison with experiments, Diploma Thesis, Faculty of 182 Ứng dụng mơ hình luồng khí vùng khơng gian giới hạn cơng trình kiến trúc an toàn sức khoẻ Process Engineering, Hermann-Rietschel Institute, Technical University of Berlin, Germany 11 LOBSCHEID, C., and GADGIL, A J (2002), “Mixing of a point-source indoor pollutant: numerical predictions and comparion with experiments,” Proceedings, Indoor Air 2002, IV, pp 223-228 Also Lawrence Berkeley National Laboratory Report LBNL- 49457 12 MORA L., AND GADGIL, A J., (2002), “Theoretical study of pollutant mixing in rooms induced by occupancy,” Proceedings, RoomVent 2002, Copenhagen, Denmark Also, Lawrence Berkeley National Laboratory report LBNL-49730 13 OZKAYNAK, H., RYAN, P B., ALLEN, G A., AND TURNER, W A., (1982), “Indoor Air Quality Modeling: Compartment Approach with Reactive Chemistry”, Environment International, 8, 461-471 14 PATANKAR, S V (1980) Numerical Heat Transfer and Fluid Flow Hemisphere Publishing, McGraw-Hill, New York 15 ROY, S., KELSO, R.M., and BAKER, A.J (1994) “An efficient CFD algorithm for the prediction of contaminant dispersion in room air motion”, ASHRAE Transactions, 100(2), 980-987 16 STUART BATTERMAN, GINA HATZIVASILIS AND CHUNRONG JIA “Concentrations and emissions of gasoline and other vapors from residential vehicle garages”.Environmental Health Sciences, University of Michigan, Ann Arbor, MI 48109-2029, USA.Migration of volatile organic compounds from attached garages to residences: A major exposure source 17 STUART BATTERMANa, , CHUNRONG JIAa and GINA HATZIVASILISb a Environmental Health Sciences, 1420 Washington Heights, Room 6075, University of Michigan, Ann Arbor, MI 48019-2029, USA.bEcophysics Inc., Ann Arbor, MI 481082200, USA 18 YAGHOUBI, M.A., KNAPPMILLER, K.D., and KIRKPATRICK, A.T (1995).Threedimensional numerical simulation of air contamination dispersal in a room”, ASHRAE Transactions, 101(2), 1031-1040 19 YOUNGRYEL RYUA, SEOGCHEOL KIMB, DOWON LEEA ‘The influence of wind flows on thermal comfort in the Daechung of a traditional Korean 183 Ứng dụng mơ hình luồng khí vùng khơng gian giới hạn cơng trình kiến trúc an tồn sức khoẻ house’ Department of Environmental Planning, Graduate School of Environmental Studies, Seoul National University, Republic of Korea Building and Enviroment./ 20 V.N BOGOSLOPSKI, V NOVODILOV “Cấp nhiệt thơng gió Tập II Nhà xuất xây dựng, Matxcơva 1978 (Tiếng Nga) 21 BÙI VẠN TRÂN Nhiệt kiến trúc -Giáo trình ðH XD Hà Nội,1977 22 HỒNG HẢI VĨ Các biện pháp chống nóng, chống nhiễm khơng khí ngồi nhà cơng nghiệp NXB KHKT Hà Nội, 1993 23 TERRY S BOUTET “Thơng gió tự nhiên nhà ở” NXB VHTT, HCM, 2005 (Sách dịch KTS Hà Nhật Tân) 24 Th.S ðỖ HỒNG OANH “Ơ nhiễm nhà gây hại nhiều nhiễm ngồi trời” (2006), Sở Tài Nguyên Môi Trường Thành Phố Hồ Chí Minh 25 QCVN 05:2009/BTNMT (Quy chuẩn chất lượng mơi trường khơng khí xung quanh) 26 QCVN 06:2009/BTNMT (Quy chuẩn số chất khí độc hại mơi trường khơng khí xung quanh) 27 QC EPA (Cục bảo vệ Môi Trường Hoa Kỳ Environment Protection Agency) 28 Tiêu chuẩn xây dựng Việt Nam NXB XD Hà Noäi,1997 184