한국연소학회지, Vol.11, No.1, pp.1-10(2006) 1 유류연소발전용보일러에서공기공급계통의불균일성에관한실험적연구 고영건 * 최상민 ** 김영주 *** An Experimental Study on the Non-Uniform Flow Distribution in the Windbox of an Oil-Fired Boiler Young Gun Go, Sangmin Choi and Young-zoo Kim ABSTRACT Oil-fired power plant usually uses several burners and combustion air is supplied to each burner through the complicated duct which is called windbox. A windbox should be designed to supply combustion air to each burner uniformly but, due to the complicated duct shape, flow distribution in the windbox is unbalanced and non-uniform supplies of combustion air are induced by these unbalanced flows in the windbox. These flow patterns tend to make flame unstable, increase the formation of pollutants and lower the overall combustion efficiency. To prevent these disadvantages, flow patterns in the windbox should be investigated for the uniform flow distribution. In this study, computational simulation method was used to investigate the flow distribution in a windbox and measured the velocities at the exit of burners in a real windbox and model tests to compare with CFD results. The results show two significant flow patterns. One is that the flow rates of each burner are different from each other and this means that all burners operate in different conditions of air to fuel ratio. The other is that the flow distribution at the exit of each burner is not axi-symmetric although the burner shape is axi-symmetric. Additionally some modifications of windbox shape and installation of baffles were proposed to make the uniform flow in the windox. Key Words : Windbox, Oil-Fired Burner, Power Plant, Flow Distribution, CFD 1. 서론 국내화력발전설비중약 16 % 의설비용량을차지하는중유화력발전은 500 MW 급표준석탄화력발전이급속히성장하는과정에서상대적으로등한시되면서 1997 년이후지속적으로발전량이감소추세에있다 [1]. 그러나중유화력발전은뛰어난설비가용성과안정적인동력원확보의일환으로실시되고있는연료다변화정 * 한국과학기술원기계공학과대학원 ** 한국과학기술원기계공학과 *** 한전전력연구원 연락저자, smchoi@kaist.ac.kr 책의큰부분을담당하고있다. 정부의환경규제의강화정책으로인해발전업체및관련기관에서도발전효율증진과함께유해 NOx 저감등연소생성물저감에관한연구가활발히진행되고있다 [2-5]. 중유화력발전용보일러에는다수의버너가내부에설치되어있고연료공급은 Fig. 1 과같은계통을통해이루어진다. 연료의이송성을높이기위해이송라인중간에서 120 정도로가열한뒤버너로공급된다. 전체연료량을하나의조절밸브로총괄적으로조절하기때문에각버너로공급되는연료량은개별적으로조절할수없는구조를국내일부의중유화력발전소에서갖고있다. 이로인해각버너로공급되는연료량
2 고영건 최상민 김영주 Heating Steam Shut Off Valve Multi - Oil Burner Fuel Oil Daily Tank Oil Pump Oil Heater Oil Flow Control Valve Oil Pressure Control Valve Condensate Oil Recirculation Line Fig. 1 Schematic diagram of fuel supply system in the oil-fired burner 0.872 m (c) (b) (a) 0.953 m Fig. 2 Outline of a calculated windbox 이서로크게차이를보이고있으며이는각버너의연소환경이버너별로다르다는것을의미하게된다 [6,7]. 또한화력발전에서는다수의버너를사용하기때문에연료뿐만아니라연소용공기도각버너로적절히분배하게되며, 발전용보일러에서는 Windbox 라는구조물을통하여공급된다. 그기본구조는 Fig. 2 와같이버너가설치된보일러의외벽을공기가흐르는큰덕트형태의구조물이둘러싸여있는구조이다. 압입송풍기에서공급된공기는재생식공기예열기를통하여예열된뒤 Windbox 로투입된다. 버너로공급되는연소용공기의유동을안정화시키기위한 Air Foil 과댐퍼가 Windbox 에설치되어있으며각버너로적절한양의공기를공급하기위해다수의격막에의해내부가구분된다. Windbox 의통상적인구조는 Fig. 2 에서보듯이입구에서부터버너출구까지유로가급격히변하는부분이존재하기때문에내부유동이불균일하게형성되며, 그결과로각버너로공급되는공기유동이불안정해지게된다. 또한 Windbox 내부에는각버너로공급되는공기유량을자체적으로제어할수있는 (a) : Air register vane (b) : Atomizer & Impeller (c) : Burner tile & Burner throat Fig. 3 Outline of an oil-fired burner 방법이없다. 이로인해각버너의위치에따라공급되는공기유량이편차가발생하게되며이런요인들이보일러내연소환경을나빠지게하는요인이될수있다. 중유화력발전에사용되는버너의기본적인형태는공기유량을제어하는 Air register vane, 화염을안정적으로유지시키는 Impeller(Flame Stabilizer), 버너의외관을구성하는버너타일그리고연료유를무화시키는노즐로구성된다. 여기에 NOx 의생성을줄이기위한다단연소방식을도입하기도하며이런다단연소에관한연구는 Air register vane 등에의해형성되는선회강도와연관되어 Syred[8] 등에의해많은연구들이진행되어왔다. 중유화력발전에사용되는버너의형태는기본적으로축대칭적인형상을하고있지만앞에서언급했듯이 Windbox 내부의불균일한유동으로인해버너로공급되는공기는
유류연소발전용보일러에서공기공급계통의불균일성에관한실험적연구 3 중심축을기준으로비대칭적인형태로공급될수있다 [9]. 이에따라본연구에서는실제운전되고있는중유화력발전소에대해 Windbox 를통해각버너로공급되는연소용공기량의편차와개별버너의 Air register vane 을통해공급되는유동형태가이상적인조건과얼마나차이가있는가를전산유동해석, 모형실험그리고실제발전소에서의유속측정을통해확인하고자하며 Windbox 형상변경및배플설치를통한개선방향을제시하고자한다. 2. 공기공급계통의불균일성 앞절에서도언급했듯이다수의버너가사용되는경우에는각버너별로균일한연소용공기가공급되어야각버너에서화염모양및길이, 과잉공기비, 연소효율및공해물질생성을제어할수있게된다. 단지각버너별로공급되는공기량뿐만아니라버너의 Air register vane 별로도동일하게공기를공급해야버너축을중심으로대칭적인유동형태를나타내게되어연소실내화염이안정적으로형성되게된다. 이렇듯저공해및화염안정성관점에서 Windbox 내에서의유동분포를균일하게하는것은중요한요소로작용하게된다 [3,5]. 2.1 버너별공급공기유량의불균일성 각버너별로공급되는공기량의편차가발생하게되면공기가부족한쪽의버너에서는산화제부족으로인해 CO 등의불완전연소생성물이증가하게된다. 이를방지하기위해서연소용공기공급을증가하게되면다른쪽의버너에서는 과잉공기공급으로인한온도강하및연소불안정성등의영향이나타나게된다 [3,4,10]. 2.2 개별버너에서비대칭적인공기유입 버너내부에서의유동은버너중심축을기준으로대칭적인형태를보일때가화염이가장안정적으로형성되게된다. 그러나각버너로공급되는공기유량이균일하다고해도버너의 Air regiter vane 별로들어가는공기량의편차가크게되면버너내부에서의유동형태가버너축을중심으로비대칭적으로형성되어화염안정성에나쁜영향을주게된다. 그렇기때문에버너주변으로공급되는공기분포가균일하도록 Baffle 등의추가적인장치가필요하다. 3. Windbox 내유동실험및전산해석 실제발전소에서각버너로공급되는공기유량을파악할수있으면 Windbox 내유동분포및이로인한연소불안정의원인을유추할수있게된다. 그러나실제운전되고있는발전소에서는유가상황이나전력소요에따라운전상황을수시로변화시키기때문에 Windbox 내의유동분포및각개별버너로들어가는유량을실측을통해서확인하는것은매우어려운작업이다. 그러므로실질적으로다양한조건에서 Windbox 내의유동형태를관찰할수있는방법은모형실험을통해확인할수있다. 그러나축소모형을사용한실험에서는원형과의상사성을유지하기어려우며 1:1 모형사용은실험장치규모가너무크기때문에실험에제약이따르게된다. 이런실험의한계성을보완하기위하여본연구에서는전산해석결과와비교를하였으며, 740 mm Pitot tube 485 mm Manometer X-Y Robot Control & data acquisition PC (a) Schematic diagram (b) Photo of windbox cold flow model Fig. 4 Windbox cold flow model test
4 고영건 최상민 김영주 Table 1 실험조건 원형항목단위설계조건 모 형 실 험 조 건 4/4 Load Case 1 Case 2 Case 3 Case 4 Case 5 송풍기회전수 RPM 607 760 913 1066 1218 총유량 m3/min 30,730 107 135 163 191 219 버너유량 m3/min 1,099 4.46 5.63 6.79 7.95 9.13 온도 293 20 20 20 20 20 점도 kg/m s 2.903 10-5 1.78 10-5 1.78 10-5 1.78 10-5 1.78 10-5 1.78 10-5 밀도 kg/m 3 0.6527 1.225 1.225 1.225 1.225 1.225 버너출구직경 m 0.96 0.048 0.048 0.048 0.048 0.048 출구유속 m/s 25.3 50 63 77 90 103 Re( 버너출구 ) 5.467 10 5 1.662 10 5 2.095 10 5 2.528 10 5 2.959 10 5 3.398 10 5 가있으며버너내부로들어오는연소용공기의유량및선회도는 Air register vane 개도에따라제어가된다 [11]. 3.2 실험장치 Fig. 5 Schematic diagramof burner-cell and air register vane 실제운전되는발전소에서도 Windbox 를통해버너로공급되는유량의편차를측정하여모형실험과전산해석결과의비교기준으로설정하였다. 3.1 해석대상개요 본연구대상 Windbox 의기본형태는 Fig. 2 에서보는바와같이연소로주위에공기를공급하는덕트로둘러싼형태로되어있으며 (1floor), 하부층 (2floor), (F) 및 Rear(R) 의 4 부분이격막에의해구분된다. 또한 2 개의버너가상부버너와하부버너로하나의셀을구성하며, 총 12 개의버너셀이설치되어있다. 각각의버너셀은 Windbox 바깥쪽에서버너셀을봤을때왼쪽부터 A,B,C 의기호를갖게된다. 또한각버너셀에서위쪽의상부버너는 Upper(U), 아래쪽의하부버너는 Lower(L) 버너로불린다. 그러므로하부층의 면의왼쪽버너셀은 2FA 라고표시되며 2FA 버너셀의하부버너는 2FAL 로표시된다. 버너의형상은 Fig. 3 과같이버너중앙에 Impeller 와연료유를분사하는 Atomizer 실험에사용된모형은실물의 1/20 scale 로축소된모형을사용하였다. Fig. 4 에서보듯이장치구성은크게 Windbox 모형, 공기를공급하기위한송풍기, 압력측정장치, 측정포트의위치를이동시키는 X/Y Robot 그리고제어및데이터획득을위한 PC 로구성된다. 송풍기는인버터를통해서유량을조절할수있으며 Windbox 의입구전단에서양쪽으로유동이분기되어공급된다. 공급된공기는 Windbox 내부덕트를거쳐최종적으로버너의 Air register vane 을통해공급되며, 버너의모형형상은 Fig. 5 와같다. 공급공기유량은피토관과마노미터를이용해서측정을하며, 피토관은 PC 에의해제어되는 X/Y Robot 에의해이동되어각버너의유량을측정하여 PC 에저장된다. 본실험에서는작동유체를원형과똑같은공기를사용하기때문에모형과원형사이의레이놀즈수를상사하는것은현실적으로불가능하다. 그러나레이놀즈수가충분히발달된난류영역에속하는경우에는레이놀즈수값에큰영향을받지는않기때문에원형과같은차수의레이놀즈수를유지하면어느정도원형과실물사이의상사조건을만족한다고가정을할수있다. 3.3 실험및전산해석조건 실험조건은상온의공기공급유량을조절하여 Table 1에서와같이 107 m 3 /min ~ 219 m 3 /min 사이의유량조건에서다섯가지경우에대해서실험을하였다. 전산해석은실제 Windbox 형상에대해 4/4 (100%) 부하기준에대해서 Fluent 6.0을이용하
유류연소발전용보일러에서공기공급계통의불균일성에관한실험적연구 5 여해석을하였으며난류모델은 Standard κ-ε을사용하였다. 또한해석영역은 Windbox입구에연결된곡관덕트부분부터버너의 Air register vane 까지를해석대상으로설정하였다. 또한가운데버너를중심으로완벽한좌우대칭형태를보이고있기때문에 Symmetric 경계조건을사용하여전체 Windbox 중에 1/2 에해당하는영역에대해서만해석을하였다. 4. 전산해석및실험결과 4.1 Windbox 내유동분포 ( 전산해석 ) 유동분포를확인하기위한전산해석결과를 Fig. 6에나타내었다. 이그림은 쪽과 Rear 쪽에서의단면속도분포를나타낸것으로덕트가꺾이는부분들로인해서전체적으로덕트의안쪽면보다는바깥쪽에서빠른유동분포를보이고있다. 버너가설치된덕트로연결되는 Windbox의코너부분에서는 Rear쪽의유동분포가 쪽보다는불안정하며덕트아래면쪽의유속이윗면쪽보다빠르게나타나고있다. 이는 Rear 쪽으로연결되는덕트의길이가짧아유동이균일해질수있는공간이 에비해부족하기때문이다. 또한입구부분의덕트가아래쪽으로굽어져있기때문에전체적으로 Windbox 바닥면부근의유속이빠르게분포하고있다. 4.2 버너셀별유량편차 ( 전산해석및모형실험 ) 앞절의전산해석결과를보면 Windbox 내부단면을따라유량분포가 Windbox의형상에크게영향을받음을알수있으며이는각버너로공급되는연소용공기유량에도영향을주게된다. Fig. 7은모형실험을통해측정된각버너셀로공급되는유량편차를나타낸것이다. 여기서편차 (Deviation) 는각버너로공급되는유량과평균유량과의차를평균유량에대한백분율비로나타낸것이다. 각버너셀로공급되는유량분포는모든경우에서유사한경향을보이고있다. 전체적으로는 쪽으로더많은공기가공급되는것을확인할수있으며 쪽의의버너셀로조금더많은공기가공급되는반면에 Rear쪽에서는하부층으로더많이공급된다. 각버너셀별로는 쪽에서는과하부층모두에서 B열로공급되는유량이더크나 Rear쪽에서는에서는 B열이하부층에서는 A,C 열로더많은공기가공급되고있다. Table 2에는각버너출구에서의유량편차를전산해석을통해계산한결과와실험결과의일부를나타낸것이다. 발전소현장에서의측정은실제운전중인발전소에서측정을하기때문에측 원형전산해석 Rear 하부층 모형실험 (Case 3) Rear 하부층 발전소현장실험 Table 2 Air flow rate deviation upper -4.6 12.8-4.6 lower -4.5 17.1-4.5 upper -8.0 15.0-8.0 lower -16.8 8.4-16.8 upper 1.24 2.66 1.85 lower -1.26 3.00 1.45 upper -1.92-0.73-1.88 lower -2.31-0.27-2.32 upper 12.51-1.62 lower -5.99-4.9 정조건이매우제한되며이로인해 2개의버너셀에대한측정만실시하였다. 기본적으로전산해석결과와모형실험의결과는두경우모두 B 열의버너에공기가많이공급되며 A열과 C열의버너에는상대적으로적게공급되는것으로나타난다. 이외에도하부층의 Rear쪽 A, C 버너가가장적게공급이되는것도공통적으로나타났다. 그러나전산해석과모형실험결과의가장큰차이점은버너출구의유량편차가모형실험결과보다는전산해석결과에서훨씬크게예측이되고있다. 여기서실제운전되고있는발전소의일부버너에서측정된현장실험결과에서는 Windbox 중앙에위치한 B열버너로역시많은공기가공급되며유량편차크기는모형실험보다전산해석결과와비슷하게나타나고있다. 이는모형실험보다는전산해석결과의유량편차가실제상황을더정확히예측한다고볼수있다. 또한모형실험은단순히레이놀즈수만을기준으로원형과상사성을유지시켰기때문에모형의크기에비해과다한유량이공급되는조건이된다. 이로인해모형실험에서는각버너출구에서의압력손실이내부에서의압력손실보다상대적으로매우크기때문에각버너로공급되는공기유량편차가원형보다는작게측정되게된다. 전산해석과현장실험결과를기준으로판단하면각버너별로공급되는공기유량편차가크게
6 고영건 최상민 김영주 덕트안쪽의유속이바깥쪽에비해낮음 코너부분유동이 Rear 쪽보다균일한분포보임 (a) Side 덕트안쪽의유속이바깥쪽에비해낮음 코너밑부분의유속증가됨 (b) Rear Side Fig 6. Velocity distribution in the windbox (m/s) (Calculation)
유류연소발전용보일러에서공기공급계통의불균일성에관한실험적연구 7 2.5 2.0 Case1 Case2 Case3 Case4 Case5 Deviation (%) 1.5 1.0 0.5 0.0-0.5 0.5 0.0 1FA 1FB 1FC 2FA 2FB 2FC (a) 1RC 1RB 1RA 2RC 2RB 2RA Deviation (%) -0.5-1.0-1.5-2.0 (b) Fig. 7 Flow deviation w.r.t burners (a) side, (b) Rear side (Model Test) 발생함으로전체적인연소효율, 오염물질생성등의문제가발생할수있으며원형 Windbox 는구조적으로개선할여지가있다고할수있다. 4.3 개별버너에서비대칭적인공기유입 Figure 8 은전산해석을통해각버너의 Air register vane 에서의공기유량분포를나타낸것으로하나의버너셀을구성하는상부와하부버너를각각의그래프에나타낸것이다. 모든버너의공통적인특징으로각버너셀의상부버너와하부버너가중심축을기준으로상하로대칭적인유동형태를보이고있다. 이는버너셀의위쪽공간과아래쪽공간보다버너셀을이루는두버너사이의공간이협소하기때문에이부분으로유입되는공기유량은적고반대로버너셀의위쪽과아래쪽으로들어오는공기유량은크기때문에모든버너들에서상부버너와하부버너는버너셀이중심축을기준으로서로반대되는형태의유량분포를공통적으로보이고있다. 유입되는공기유량이큰 B 열버너의경우는버너상하의유동차이가더욱크게나고있다. A 열과 C 열의버너는버너상부와하부뿐만아니라좌우로도유동이불균일하게나타나고있다. 쪽의 A 열은버너중심에서왼쪽에서, Rear 쪽의 C 열은버너중심에서오른쪽에서유입 Fig. 8 Air flow shape at the exits of each burner (Calculation) 유량이크게나타나고있다. 즉 Windbox의코너부분쪽으로는유속이줄어들기때문에 Windbox의코너쪽보다는중심쪽에위치한 Air register vane을통해유입되는유량이크게나타나고있다. 이는 Fig. 9와같이버너출구에서의단면속도분포를모형실험을통해측정한결과에의해서도확인할수있다. Windbox 중앙에위치한 B열버너의경우는코너에위치한 A,C열버너에비해비교적좌우대칭적인형태를띄고있으며 A,C열버너는상하좌우로비대칭적인유속분포를보이고있다. 이렇게 Air register vane별로유입되는공기유량의차이가발생하게되면버너내부에서도버너축을중심으로균일한유동분포가나타나지
8 고영건 최상민 김영주 Fig. 9 Air flow contour at the exits of each burner (Model Test) for Case 3 않고한쪽으로유동이치우치게되므로버너내에서연료와공기의분포가균일하지않게되어연소가불안정해지는원인이될수있다. 5. 구조변경을통한 Windbox 내유동개선 5.1 구조변경 기존 Windbox 에서는각버너로공급되는공기유량편차가존재하고개별버너에서도비축대칭적으로공기가공급되고있는것을확인할수있다. 이를개선하기위해 Windbox 의구조를 Fig. 10 과같이변경하여개선방향을검토하였다. Case M-1 은과하부층의 Windbox 중심부덕트높이를 10 % 줄여 B 열버너부분의단면적을줄인형태이다. 이는대상 Widbox 에서 B 열버너로공기가많이공급되기때문에 B 열부분의유로를좁게하여 B 열로공급되는공기유량을줄이기위한방안이다. Case M-2 는 Windbox 의코너에폭 40 cm 의배플을설치하고 A,C 열버너와 B 열버너사이에는높이 50 cm 의 Baffle 을덕트의천정과바닥에설치한경우이다. 코너에설치된 Baffle 은덕트내유동이 Windbox 의코너부분을지나면서덕트내바깥면쪽면의속도가높아지기는것을상쇄하기위한것이며, A,C 열버너와 B 열버너사이에위치한 Baffle 은유동이 A,C 열버너를지나면서덕트의천정과바닥에서유속이빨라지기때문에이를상쇄하기위한것이며, Case M-4 는코너부분의배플의위치를조절한것이다. 5.2 구조변경해석결과 Table 3 은변경된형상의 Windbox 에서의 쪽버너들의출구에서의유량편차를나타낸것이며, Fig. 11 은각버너출구로공급되는공기유동형태를나타낸것이다. Case M-1 의경우는 Windbox 중심부의단면이축소됨으로인한압력손실증가로 B 열버너로공급되던유량이줄어드는효과를보이기는하나 Windbox 내부에서의공기유동이원형에서의유동형태를그대로유지하기때문에여전히높은수준의유량편차를보인다. 또한 Windbox 중앙에서의좁은유로로인해 B 열의상부와하부버너는원형의경우보다더상하로비대칭적인유동형태를보이고있다. 그러나 Windbox 형상은그대로두고내부에배플을설치한 Case M-2 과 M-3 의경우도 B 열버너로더많은공기가공급되지만원형과 Case M-1 에비해편차가크게줄어들다. 이는 Windbox 내부유동에서유량이큰부분에배플설치러인한압력손실증가에의해직접적으로유량을통제하기때문에다른경우들에비해각버너로비교적균일하게공기가공급된다고할수있다. 그러나 Case M-3 는코너에설치된배플이 A 와 C 열버너에근접하여위치하기때문에배플이설치된쪽의 Air register vane 으로는공기공급이오히려제한된다. 이로인해버너출구의유동형태는가장비대칭적으로나타나게된다. Windbox 의개선방향을유량편차와유동형태의대칭성기준으로판단하면덕트의형상을변
유류연소발전용보일러에서공기공급계통의불균일성에관한실험적연구 9 Table 3 Air flow rate deviation of modified windbox(calculation) Case M-1 Case M-2 upper -4.5 14.8-4.5 lower -2.8 6.5-2.8 upper 1.6 6.9 1.6 lower -3.5 8.2-3.5 Case M-3 upper -2.6 6.3-2.6 lower 0.3 7.6 0.3 Fig. 10 Shapes of Modified Windbox 형하기보다는 Windbox 내부의적절한위치에배플을 Case M-2 의경우가원형의 Windbox 를개선할대안이라할수있다. 6. 결론 다수의버너를사용하는발전용보일러에서는각버너별로공급되는공기유량의차이가발생함으로인해개별버너주변에서의유동이불안정할수있다. 이로인하여버너의연소불안정성을초래할수있으며전체적인연소효율저감과공 Fig. 11 Air flow shape at the exits of each burner (Modified windbox, Calculation) 해물질발생이증가할수있다. 본연구에서는실제운전되고있는발전용보일러의 Windbox 를통해각버너로공급되는공기유량을축소모
10 고영건 최상민 김영주 형실험과전산해석방법을사용하여예측하였다. 전체적인유동의경향은 Windbox 의 Rear 쪽보다는 쪽으로공기의공급이많으며주로 Windbox 의대칭축에위치한 B 열로공기의공급이많이되고있는경향을보여주고있다. 이런경향은공급층의유량변화와는상관없이동일한경향을보여주고있다. 또한버너출구에서형성된연소용공기유동형태는보일러에서버너의위치에따라버너중심축을기준으로상하좌우로비대칭적인유동형태를보이고있다. 이런공기공급의불균일성은 Windbox 내부에비교적간단한배플을설치하여유동의불균일성을효과적으로제어할수있다. 후 기 본연구는산업자원부전력기반기술개발연구사업으로선정한 중유화력분진배출저감에관한연구 와연소기술연구센터 (Combustion Engi -neering Research Center) 의연구비지원에의해수행되었습니다. 참고문헌 [1] 주요통계지표종합, 한국전력, 2004. [2] 신명철, 안재현, 김세원, 다단공기공급저 NOx 버너의선회유동및연소특성에관한실험적연구, Vol.8, No.1, 2003, pp. 25-35. [3] E. Charles and J. Baukal, Industrial Burners Handbook, CRC Press, 2004. [4] M. G. Caarvalho, W. A. Fiveland, F. C. Lockwood and C. Papadopoulos, "Combustion Technologies For a Clean Environment," Gordon and Breach Publishers, 1991, pp.11-130, 203-258. [5] J. G. Singer, "Combustion Fossil Power," Combustion Engineering INC, 1991. [6] 중유화력분진배출저감에관한연구 (1 차년도중간보고서 ), 산업자원부, 2003.8. [7] 중유화력분진배출저감에관한연구, 산업자원부, 2004.11. [8] N. Syred and J. M. Beer, "Combustion in Swirling Flow : A Review," Combustion and Flame, Vol. 23. 1974, pp.143-201. [9] 김병윤, 백희석, 산업용보일러윈드박스유동해석. 2002 Fluent Users's Group Meeting 2002. [10] E. Charles, J. Baukal and E. S. Robert, "The John Zink Combustion Handbook," CRC Press, 2001, pp. 547-586. [11] 평택화력 1-4 호기기본운전지침서 (I)," 한국전력공사평택화력발전처, 1998, pp.1-27.