本論文研究是利用ANSYS CFX軟體進行遮煙試驗爐之內部流場之數值模擬，分析防火門在常溫與中溫環境下之空氣洩漏量，並探討離心式風機和電熱管對於溫升曲線的影響，本研究紊流模型採用K-Epsilon 紊流模型，網格系統選用六面體網格和四面體網格組成，以建立可與實驗比對並提供可靠的參考結果。
程式驗證方面，以三維空穴流、三維自然對流等，並與文獻之數據進行比對，驗證ANSYS CFX在流場流動與熱傳導、熱對流之物理性質上的可靠性與準確性。本模擬之實際的物理模型是以國立成功大學防火安全研究中心所提供的試驗爐來建立，主要探討遮煙試驗爐流場在艙體安裝測試門扇後常溫時的洩漏量是否與實驗結果符合，以及利用電熱管加熱維持實驗之中溫環境，並探討電熱管與離心式風機對於該流場溫升效果之影響。在洩漏量方面，數值模擬之結果與實驗模擬之結果相當一致其誤差為5%，在流場溫升方面，欲維持爐內的溫度，風機的設計和增加風機轉速的方法有助於提高及維持爐內溫度。

The goal of paper uses ANSYS CFX to investigate the fluid fields in air leakage test furnace. We investigate the relationship among temperature trends and the turbine rotational speed, electric heating pipe temperature. In the simulation, we use high resolution scheme and K-Epsilon turbulence model to solve the incompressible Navier-Stokes equations. The mixed grid system which combines Hexahedral and prism meshes is used.
First, several test problems, 3-D cavity flow; natural convection; and air leakage test for a plate with a hole, are simulated to understand the capability of the commercial program about the air leakage test. The numerical results are compared well with the experimental data performed by NCKU Fire Safety Research Center. Finally, the flow field and temperature in the furnace are investigated in detailed. We find out that the design angle and rotation speed of wind turbine are very important parameter to elevate temperature in the furnace.

INTRODUCTION
Fires often cause casualties due to smoke and heat. Heat is the most frightening factors. How to effectively control the proliferation of fire and smoke is very important issue. Fire doors are used to prevent the spread of smoke within a certain time. Basically a fire door testing may spend hundreds of thousands. Using ANSYS CFX can predict the performance of door and save. The goal of paper is to simulate the fluid field in air leakage test furnace by using ANSYS CFX. We investigate the temperature trends by changing the turbine rotational speed and electric heating pipe temperature.

RESULTS AND DISCUSSION
In the simulation, we use high resolution scheme and K-Epsilon turbulence model to solve the incompressible Navier-Stokes equations. The mixed grid system which combines Hexahedral and prism meshes is used. First, several test problems, 3-D cavity flow; natural convection; and air leakage test for a plate with a hole, are simulated to understand the capability of the commercial program. About the air leakage test, the numerical results are compared well with the experimental data performed by NCKU Fire Safety Research Center. Finally, the flow field and temperature distribution in the furnace are investigated in detailed.

CONCLUSION
The gas flow is roughly trend by buoyancy and wind turbine design. Heat will accumulate at the right part and are mainly affected by the wind turbine position. We find out that the design angle and rotation speed of wind turbine are very important parameter to elevate temperature in the furnace. In the furnace, temperature trends by wind turbine and heating pipes. Some lower temperature gas was caused by the door and the gas flowing through the electro-thermal tube. The air flow speed is so fast that heating electro-thermal tube can not effectively heating.

[1] Chinese National Standards CNS15038,Method of test for evaluating smoke control performance of doors buildings, The bureau of standards, Metrology and Inspection, M.O.E.A. ,ROC, 2006.
[2] ANSYS CFX, Inc. ANSYS 12.1,2009
[3]陳彥章, 防火安全用防火爐流場之數值模擬,國立成功大學航空太空工程所碩士論文,2008.
[4]杜博文, 循環式防火爐流場之數值分析,國立成功大學航空太空工程所碩士論文,2009.
[5]羅人豪, 耐火試驗爐流場與防火門熱傳現象之數值模擬,國立成功大學航空太空工程所碩士論文,2012.
[6]CEN, Eurcode 1:Action on structures part 1.2:General Action –
Actions on structures exposed to fire BS EN 1991 - 2:2002,Brussels:
CEN, European committee for standardissation,2002.
[7] I.B. Eck , Fans: design and operation of centrifugal, axial flow,
and cross- flow fans,Pergamon Press, 1973.
[8] R.H. Aungier ,Centrifugal compressors : a strategy for aero dynamic
design and analysis , ASME Press,2000.
[9]A.F. Mills, Basic heat and mass transfer,Irwin,UAS,1995.
[10]C. Shu, L. Wang, and Y.T. Chew, Numerical computation of three-
dimensional incompressible Navier-Stokes equation in primitive
variable from DQ method , Int. J. Numer. Meth. Fluids, Vol 43,
pp.345-368, 2003.
[11]T. Fusegi, J.M. Hyun, K. Kuwahara, and B. Farouk, A numerical study
of three-dimensional natural convection in a differentially heated cubical enclosure, Int. Journal. Heat and Mass Transfer,Vol34, 76,pp.1543-1557, 1990.
[12] National Fire Protection Association: NFPA 252, Standard methods
of fire tests of door assemblies, 2008.
[13] UL 10C,UL Standard for safety positive pressure fire tests of door
assemblies, First edition, Underwriters Laboratories Inc. Northbrook, IL USA,1998.
[14]T. Ohmura, M. Tsuboi, M. Onodera, and T. Tomimura, Specific heat
measurement of high temperature thermal insulations by drop
calorimeter method, Int. J. Thermophysics, Vol. 24, 2, 559-575,
2003.
[15]International Standard ISO 834-1, Fire resistance tests-elements of
building construction-Part 1: General requirements , 1999.
[16]British Standards Institution BS476, Fire tests on building materials and structures- part 20 : Methods for determination of the resistance of element of construction（general principles）, 1987.
[17]American Society for Testing and Materials Standard : ASTM E152 ,Methods of fire tests of door assemblies, 1995.
[18]National Fire Protection Association: NFPA 252, Standard methods of fire tests of door assemblies, 2008.
[19]JIS A1311,Method of fire protecting test of fire door for buildings, Japan ,1994.
[20]Draft for public Comment AS 1530-4, Methods for fire tests on building materials, components and structures-Part 4：Fire-resistance tests of elements of building construction, Draft for Public Comment Australian Standard , 2004.
[21]H.H Lee, Finite element simulations with ANSYS Workbench 13, SDC, Mission KS, 2011.