DBPIA-NURIMEDIA

Similar documents
14.531~539(08-037).fm

<35335FBCDBC7D1C1A42DB8E2B8AEBDBAC5CDC0C720C0FCB1E2C0FB20C6AFBCBA20BAD0BCAE2E687770>

저작자표시 - 비영리 - 변경금지 2.0 대한민국 이용자는아래의조건을따르는경우에한하여자유롭게 이저작물을복제, 배포, 전송, 전시, 공연및방송할수있습니다. 다음과같은조건을따라야합니다 : 저작자표시. 귀하는원저작자를표시하여야합니다. 비영리. 귀하는이저작물을영리목적으로이용할

Microsoft Word - KSR2012A021.doc

THE JOURNAL OF KOREAN INSTITUTE OF ELECTROMAGNETIC ENGINEERING AND SCIENCE. vol. 29, no. 10, Oct ,,. 0.5 %.., cm mm FR4 (ε r =4.4)

<31325FB1E8B0E6BCBA2E687770>

부문별 에너지원 수요의 변동특성 및 공통변동에 미치는 거시적 요인들의 영향력 분석

272 石 堂 論 叢 49집 기꾼이 많이 확인된 결과라 할 수 있다. 그리고 이야기의 유형이 가족 담, 도깨비담, 동물담, 지명유래담 등으로 한정되어 있음도 확인하였 다. 전국적인 광포성을 보이는 이인담이나 저승담, 지혜담 등이 많이 조사되지 않은 점도 특징이다. 아울

09È«¼®¿µ 5~152s

THE JOURNAL OF KOREAN INSTITUTE OF ELECTROMAGNETIC ENGINEERING AND SCIENCE Jun.; 27(6),

DBPIA-NURIMEDIA

DBPIA-NURIMEDIA

THE JOURNAL OF KOREAN INSTITUTE OF ELECTROMAGNETIC ENGINEERING AND SCIENCE Dec.; 27(12),

04-다시_고속철도61~80p

ePapyrus PDF Document


12(4) 10.fm

Analysis of objective and error source of ski technical championship Jin Su Seok 1, Seoung ki Kang 1 *, Jae Hyung Lee 1, & Won Il Son 2 1 yong in Univ

에너지경제연구 Korean Energy Economic Review Volume 17, Number 2, September 2018 : pp. 1~29 정책 용도별특성을고려한도시가스수요함수의 추정 :, ARDL,,, C4, Q4-1 -

에너지경제연구 제13권 제1호

00-CD..

인문사회과학기술융합학회

03-서연옥.hwp

DBPIA-NURIMEDIA

2. 수치시뮤레이션 2.1 기본방정식과수치조건 기본방정식은 Navier-Stokes 방정식이며 FEM 수치기법으로이산화하여구조격자를만들어계산을수행하였다. k- 을사용한수송방정식은 t (ρε)+ (ρεu x i )= i x j [( μ+ μ t σ ε ) ε + C 1ε

09권오설_ok.hwp

°í¼®ÁÖ Ãâ·Â

PJTROHMPCJPS.hwp

DBPIA-NURIMEDIA


책임연구기관

232 도시행정학보 제25집 제4호 I. 서 론 1. 연구의 배경 및 목적 사회가 다원화될수록 다양성과 복합성의 요소는 증가하게 된다. 도시의 발달은 사회의 다원 화와 밀접하게 관련되어 있기 때문에 현대화된 도시는 경제, 사회, 정치 등이 복합적으로 연 계되어 있어 특

878 Yu Kim, Dongjae Kim 지막 용량수준까지도 멈춤 규칙이 만족되지 않아 시행이 종료되지 않는 경우에는 MTD의 추정이 불가 능하다는 단점이 있다. 최근 이 SM방법의 단점을 보완하기 위해 O Quigley 등 (1990)이 제안한 CRM(Continu


10(3)-10.fm

6.24-9년 6월

???? 1

M&S. EDISON. e-science EDISON_CFD[2] EDISON. (geometry kernel software),,, EDISON_CFD.,.. 4. e-science EDISON M&S.., EDISON , b f. Find to min

ePapyrus PDF Document

THE JOURNAL OF KOREAN INSTITUTE OF ELECTROMAGNETIC ENGINEERING AND SCIENCE Nov.; 26(11),

DBPIA-NURIMEDIA

82 제 1 발표장 (2 일금 ) CFD 응용 [V] 이남훈 1*, 류태광 2 NUMERICAL VERIFICATION OF SHAKE TABLE TEST FOR THE LIQUID STORAGE TANK N. Lee and T. Yoo 1.,.,,,.,. Baek e

DBPIA-NURIMEDIA

THE JOURNAL OF KOREAN INSTITUTE OF ELECTROMAGNETIC ENGINEERING AND SCIENCE Mar.; 28(3),

High Resolution Disparity Map Generation Using TOF Depth Camera In this paper, we propose a high-resolution disparity map generation method using a lo

116 제 2 발표장 (2 일금 ) 특별세션 - 고성능컴퓨팅 박근태 1, 최해천 1*, 최석호 2, 사용철 2, 권오경 3 ANALYSIS AND CONTROL OF TIP-LEAKAGE FLOW IN AN AXIAL FLOW FAN USING LARGE EDDY SI

<303320B9DABCB1C8A328BFCF292E687770>

untitled

Journal of Educational Innovation Research 2017, Vol. 27, No. 1, pp DOI: * The

창의-열유체-3

<353420B1C7B9CCB6F52DC1F5B0ADC7F6BDC7C0BB20C0CCBFEBC7D120BEC6B5BFB1B3C0B0C7C1B7CEB1D7B7A52E687770>

< B3E2BFF8BAB828C8AFB0E629312E687770>

[ 물리 ] 과학고 R&E 결과보고서 유체내에서물체의마찰력에미치는 표면무늬에대한연구 연구기간 : ~ 연구책임자 : 홍순철 ( 울산대학교 ) 지도교사 : 김영미 ( 울산과학고 ) 참여학생 : 김형규 ( 울산과학고 ) 노준영 (

< C6AFC1FD28B1C7C7F5C1DF292E687770>

THE JOURNAL OF KOREAN INSTITUTE OF ELECTROMAGNETIC ENGINEERING AND SCIENCE Feb.; 28(2),

Journal of Educational Innovation Research 2017, Vol. 27, No. 3, pp DOI: (NCS) Method of Con

08김현휘_ok.hwp

<303120C1A6B5CEC8A32D F696C20C1D6BAAFBFA1BCADC0C720C3FEB7F92E687770>

PowerPoint 프레젠테이션

THE JOURNAL OF KOREAN INSTITUTE OF ELECTROMAGNETIC ENGINEERING AND SCIENCE. vol. 29, no. 6, Jun Rate). STAP(Space-Time Adaptive Processing)., -

에너지경제연구 Korean Energy Economic Review Volume 11, Number 2, September 2012 : pp. 1~26 실물옵션을이용한해상풍력실증단지 사업의경제성평가 1

(1 일목 ) 제 3 발표장 47 수치기법 [I] 이은택 1, 안형택 2* SIMULATION ON FLOW PAST A CIRCULAR CYLINDER USING UNSTRUCTURED MESH BASED INCOMPRESSIBLE FLUID SOLVER(ULSAN3

Journal of Educational Innovation Research 2018, Vol. 28, No. 1, pp DOI: A study on Characte

Slide 1

WOMA Pumps - Z Line

박선영무선충전-내지

Introduction Capillarity( ) (flow ceased) Capillary effect ( ) surface and colloid science, coalescence process,

THE JOURNAL OF KOREAN INSTITUTE OF ELECTROMAGNETIC ENGINEERING AND SCIENCE Sep.; 26(10),

<353920C0B1B1E2BFEB2DB0E6B0F1C0DCB1B320BBF3BACEB1B8C1B6C0C720C8DA2E687770>

DBPIA-NURIMEDIA

Microsoft Word - kai1.doc

DBPIA-NURIMEDIA

exp

¼º¿øÁø Ãâ·Â-1


<C7D1B1B9B1B3C0B0B0B3B9DFBFF85FC7D1B1B9B1B3C0B05F3430B1C733C8A35FC5EBC7D5BABB28C3D6C1BE292DC7A5C1F6C6F7C7D42E687770>

Journal of Educational Innovation Research 2017, Vol. 27, No. 4, pp DOI: * A Study on Teache

example code are examined in this stage The low pressure pressurizer reactor trip module of the Plant Protection System was programmed as subject for

Journal of Educational Innovation Research 2018, Vol. 28, No. 4, pp DOI: * A S

Journal of the Korea Academia-Industrial cooperation Society Vol. 15, No. 3 pp , ISSN 197

Microsoft Word - KSR2013A320


08.hwp

Æ÷Àå82š

(5차 편집).hwp

<31372DB9DABAB4C8A32E687770>

Microsoft PowerPoint - dev6_TCAD.ppt [호환 모드]

00내지1번2번

04김호걸(39~50)ok

±èÇö¿í Ãâ·Â

04_이근원_21~27.hwp

PowerPoint 프레젠테이션

Æ÷Àå½Ã¼³94š

Microsoft PowerPoint - 26.pptx

., (, 2000;, 1993;,,, 1994), () 65, 4 51, (,, ). 33, 4 30, 23 3 (, ) () () 25, (),,,, (,,, 2015b). 1 5,

Journal of Educational Innovation Research 2018, Vol. 28, No. 3, pp DOI: * Strenghening the Cap

<313920C0CCB1E2BFF82E687770>

03.hwp

에너지경제연구 Korean Energy Economic Review Volume 14, Number 2, September 2015 : pp. 99~126 산유국의재생에너지정책결정요인분석 1) 99


Transcription:

한국해양환경. 에너지학회지 Journal of the Korean Society for Marine Environment and Energy Vol. 16, No. 2. pp. 61-70, May 2013 http://dx.doi.org/10.7846/jkosmee.2013.16.2.61 조류발전용수직축터빈의유체동력학적 3 차원효과에관한연구 현범수 이정기한국해양대학교조선해양시스템공학부 Hydrodynamic Aspects on Three-dimensional Effects of Vertical-axis Tidal Stream Turbine B. S. Hyun and J. K. Lee Div. of Naval Architecture & Ocean Systems Eng., Korea Maritime Univ., Busan 606-792, Korea 요 약 조류발전용수직축터빈 (VAT) 의 3차원적인유체동력학적성능을효과적으로예측할수있는연구를진행하였다. 수치해석은 2차원과 3차원으로수행하였으며이를통해물리적유동현상의차이를파악하였다. 3차원효과는주로날개끝단에서발생하는날개끝보오텍스가주된원인으로서, 이로인해터빈날개가내어주는양력이손실되고회전하는터빈은토크가감소하였다. 이러한현상은본연구에서채택한통상적인수직축터빈의스팬-직경비범위에서상당한수준으로나타남을확인하였다. 본연구에서는대상터빈을선정하고 2차원으로성능해석후 3차원효과를보정하는비교적간단하고효과적인방법을제안하였다. Abstract Hydrodynamic aspects on three-dimensional effects were investigated in this study for simple and convenient conversion of tidal stream energy using a Vertical-Axis Turbine (VAT). Numerical approach was made to reveal the differences of flow physics between 2-D estimation and rigorous 3-D simulation. It was shown that the 3-D effects were dominant mainly due to the variation of tip vortices around the tip region of rotor blade, causing the loss of lift for steadily translating hydrofoil and the reduction of torque for rotating turbine blade. The 3-D effect was found to be rather prominent for the typical VATs considered in this paper. Simple and yet efficient 2-D approach with the correction of its three-dimensionality was also proposed for practical design and analysis of VAT. Keywords: Three-dimensional effects(3 차원효과 ), Vertical-axis turbine( 수직축터빈 ), Tip flow( 날개끝유동 ), Fixed hydrofoil( 고정수중익 ), Rotating Blade( 회전날개 ) 1. 서론 최근들어다양한전산유체역학기법을이용한조류발전용수직축터빈에대한연구가활발히진행되고있다. 수평축터빈과는달리다리우스형수직축터빈은형상이단순하여회전하는터빈의 2 차원수치계산으로도개략적인성능을파악할순있지만, 그결과를풍동실험이나 3차원수치계산결과와비교할때잘일치하지않는경우가자주발생하고있다. 이는터빈의회전속도에따라날개의중앙부와끝단에서다양한유동패턴이발생함에기인한것으로서, 수직축터빈의 3차원효과라고통칭한다. Corresponding author: bshyun@hhu.ac.kr 수직축터빈의 3차원효과에대한선행연구가 Table 1과같이여러연구자에의하여시도되었다. 먼저 Ye et al.[2010] 은터빈의후류와회전거동시뮬레이션에적합한 Vortex 모델을개선시킨 DVM- UBC Vortex 모델을이용하여 3차원비정상해석을수행하였는데, 난류가반영되지않은해석으로서 3차원효과가크게얻어지지않았다. Castelli et al.[2011] 은풍동실험과 RANS 모델, 공기역학모델을적용하여터빈의성능을파악하였는데 2차원 RANS 모델의경우실험결과보다 2배가량크게예측되었지만그원인에대해명확하게언급하진않았다. Howell et al.[2010] 도풍동실험과 RANS 모델을적용한수치계산을수행하였는데, TSR에따라계산및실험결과에큰차이가발생하는것은날개끝단과중축, 스트럿등의복합적인영향으로인함을추론한바있다. 하지만날개끝와 61

62 현범수 이정기 Table 1. Summary of previous studies Key research topics Open test Paraschivoiu et al.[1988] Development of BEMT code Wind tunnel test Howell et al.[2010] Application of RANS model Application of DVM-UBC Vortex model Ye et al.[2010] No bottom, no turbulence Wind tunnel test Castelli et al.[2011] Application of RANS/Aerodynamic model Fig. 3. Comparison of specific power with respect to aspect ratio (Paraschivoiu et al.[1988]). Fig. 1. Comparison of power coefficient between numerical and experimental results (Howell et al.[2010]). 류가성능에미치는영향에대한명확한단서를제공하진못하였으며, Fig. 1에서보는것과같이 TSR이높은구간에서효율이특히크게감소하는이유도밝히지못하였다. Paraschivoiu et al.[1988] 은포물선형수직축터빈의실해역실험과수치해석용 DMS기법을개발하는등의다양한방법으로연구를수행하여날개의종횡비 (AR) 에따라출력손실이달라지는것을보여주었다. Fig. 2는날개스팬방향으로속도간섭계수의변화를나타낸것으로서터빈으로유입되는속도성분이날개중앙부와끝단에서다른정도를의미한다. 이는날개끝주위에서발생하는유동은터빈중앙부와는다를것이라는것을충분히예상할수있고이로인해터빈의성능에도영향을줄것이다. Fig. 3은날개종횡비에따라달라지는비출력을나타낸것인데여기서비출력은날개의단위면적당출력을나타낸다. 일반적으로소형수직축터 빈의날개종횡비가 30 정도인점을감안하면비출력손실은 10% 정도임을알수있다. 이상의연구결과로부터 3차원효과는날개끝단의영향이주된요인으로서터빈날개의종횡비에따라그영향이달라짐을예상할수있다. 다만기존의연구들은다양한가정들에입각하여수직축터빈의 3차원적인특성을살펴보았지만불균일유동장에서비정상적인거동을보이는복잡한유동장의물리적인규명에는미치지못하였다. 본연구에서는수직축터빈의성능을효과적으로수치해석하여설계에반영하기위한방편으로서, 터빈의 2차원해석결과를바탕으로 3차원효과를보정하는방법에대하여고찰하고자한다. 수직축터빈의특성상 3차원비정상해석이필수적이지만이를직접구현하기위해서는막대한량의연산이요구되는어려움이있다. 예를들어 Intel 3.1GHz 30core 클러스터컴퓨터를사용하더라도 1 회계산에수일의계산시간이요구되는반면, 2차원해석은수시간이내의짧은시간으로해석이가능하기때문에 2차원결과를 3 차원으로보정하는방법을채택하는것이현실적이기때문이다 (Hyun et al.[2012]). 3차원효과를체계적으로관찰하기위한연구방법으로먼저고정된수중익에대하여 2차원및 3차원수치해석을수행하여날개의 3차원효과를분석해보았으며, 이를바탕으로회전하는 2차원및 3차원터빈의수치해석을수행하여회전하는터빈날개의종횡비에따른 3차원영향에대하여분석하였다. 2. 터빈선정및수치해석 Fig. 2. Upwind and downwind interference factors (Paraschivoiu et al.[1988]). 2.1 대상터빈 Han et al.[2010, 2011] 의연구결과에의하면수직축터빈설계

조류발전용수직축터빈의유체동력학적 3 차원효과에관한연구 63 Fig. 4. Schematics of 100 kw-class VAT. Table 2. Principal particulars of 100 kw-class VAT Designation Symbol Diameter(m) D 5 Number of blade Z 3 Chord length(m) c 0.44 Solidity σ 0.085 Span length to diameter ratio H/D 1.2 Blade section NACA 653-018 를위한주요변수는날개수, 솔리디티, 스팬-직경비 (H/D) 가있는데이러한주요변수를정하기전에먼저목표발전용량을설정하는것이일반적이다. 해양에너지이용기술에대한연구가활발한유럽의조류발전용터빈의발전용량추세가 100, 300, 500 kw로설정되고있는데, 수직축터빈이수평축에비하여상대적으로소형에적용되는점을감안하여 100 kw급대상터빈을선정하였다. 이는 Lee et al.[2012] 의연구에서사용하였던터빈으로서체계적인설계법에의해제시된우수한성능의터빈이다. 터빈의주요치수와형상은 Fig. 4, Table 2와같다. 2.2 수치해석 터빈의수치해석은 2차원과 3차원으로나누어수행하였다. 3차원터빈해석시격자생성의간편함을위하여스트럿은제외하였다. 난류모델은 k-ω SST model을사용하였는데이는 k-ε model에비해아주작은 y+ 가요구되기때문에그에따라격자수가증가하지만회전하는물체해석에적합한장점을가지고있다. Fig. 5는 3차원터빈격자계를나타낸것으로계산영역은터빈으로부터전방으로 5D, 후방으로 10D, 터빈높이방향으로는 4D만큼의여유를두어구성하였고터빈이회전하는영역 (Rotating zone) 과회전하지않는영역 (Non-rotating zone) 으로나누었으며총격자수는약 5,000,000개이다. 대칭이되는경계면에는 Symmetry 조건을부가하였고, 입구에 Velocity-inlet 조건, 출구에는 Pressure-outlet 조건을부가하였다. 수치계산조건은 Table 3과같으며, 로터의회전은 Jung et al.[2009] 과 Han et al.[2010] 이적용하였던슬라이딩격자기법 (Sliding mesh) 을사용하였다. 슬라이딩격자기법을적용하기위해회전영역과회전하지않는영역의경계면에인터페이스조건을사용하였는데, 이는터빈로터와그주위격자를이어주면서격 Fig. 5. Grid system around a three-dimensional turbine. Table 3. Calculation condition Conditions Turbulence model k-ω SST model (2 eq. model) Turbine rotation Sliding mesh y+ 5 Rotational speed (rad/s) 2.4 TSR 1.5~4.5 자수를감소시킬수있는강점이있다. 터빈의회전은시간간격당 2.75도 (0.02초 ) 가되도록하였고이는통상추진기성능해석시의시간간격인 5도를회전시키는것에비해더세밀한것이다. 계산의수렴여부는 TSR에따라다르지만통상 5회전이상이되면실질적으로수렴하는것으로판단되었다. 해석도구로는 ANSYS사의 FLUENT Ver.14.5 를사용하였으며 Intel 3.1Ghz 8 core 클러스터컴퓨터에서해석하였다. 수치해석의격자의존도와해석결과의검증등은 Han et al.[2010] 및 Lee et al.[2012] 가수행하였던검증과정을동일하게적용하여서로잘일치하는결과를얻었기때문에본논문에서더이상의설명은생략하기로한다. 보다구체적인내용은참고문헌의결과를참고하면된다. Fig. 6은 TSR=1.5, 3.0, 4.5에서의위상각에따른동력계수의변화를나타낸것이다. TSR에따라최대동력계수가변화하는것을잘보여주고있는데, 효율이가장우수한 TSR=3을설계점으로정하였다. 각 TSR에서 1회전시의동력계수를평균하면 Fig. 7에서와같이평균동력계수를구하게된다. Fig. 7은대상터빈에대하여 2차원및 3차원해석을수행한수치해석결과를서로비교한것이다. 3차원해석으로도출한동력계수값들이 2차원해석결과에비해작게얻어지고있다. 특히, TSR이커짐에따라그차이가커짐을알수있으며이에대한근거로는 Lee et al.[2012] 이언급한바, 높은 TSR에서동력계수가터빈의부가물 ( 중축, 스트럿등 ) 등의영향을크게받는구간 (Secondary effect region) 에해당한다고설명한바있다. 그러나본계산이스트럿을

64 현범수 이정기 Fig. 6. Comparison of power coefficient as azimuth angle (H/D=1.2). Fig. 8. Grid system for fixed hydrofoil. 위가변하게된다. α = tan 1 -------------------------- sinθ cosθ + TSR (1) Fig. 7. Comparison of power coefficient between 2D and 3D results. Table 4. Comparison of typical calculation time Number of grid Hardware Time(hr) 2-D 270,000 Intel 3.1 GHz 8core 4 3-D 5,000,000 Intel 3.1 GHz 8core 170 제외한계산임을감안한다면전체동력계수에미치는터빈날개의 3차원성에기인한다고해석하는것이보다타당할것으로여겨지는데, 이에대하여는추후정밀한조사가필요하다. 한편, Table 4에서보는것과같이 3차원터빈수치계산은 2차원에비해아주많은시간이소요된다. Ye et al.[2010] 도터빈의 3차원계산은 2차원에비해계산시간이약 40배정도더소요된다는연구과정을소개한바있다. 물론터빈로터의 z축방향으로발생하는비정상적유동현상과터빈이회전함에따른후류의영향등을포함하고있어그것을파악하는것이중요하지만, 결과를해석하는과정이쉽지않기때문에 2차원계산결과를활용해서 3차원효과를보정할수만있다면매우효과적일것이다. 3. 직진하는수중익의 3 차원효과 회전하는터빈은유입속도와로터회전속도의벡터합으로받음각을받는데, 위상각 (θ) 과 TSR에따라식 (1) 과같이받음각의범 이와같이회전하는터빈의날개는받음각변화에따라양항력계수가달라진다. 물론날개의후류영향으로실제날개가경험하는받음각은식 (1) 에서얻어지는각도와는다르며비정상유동특성을보이므로매우복잡한유동현상이된다. 하지만회전하는터빈의 3차원효과에대한연구에앞서기초연구로서, 식 (1) 에주어진다양한받음각범위에서직진하는날개에대한 2, 3차원날개의특성을먼저살펴보는것도의미있는일이다. 이는수치해석의정도검증은물론, 터빈로터의회전으로인한영향과후류의영향을무시하고순수하게날개끝유동의차이에의한양항력의손실만을살펴볼수있게한다. Fig. 8의우측상단에보여진것처럼터빈이회전함에따라순시적으로변화하는받음각들을사용하여직진하는날개에대해수치해석을수행하였고결과를 Abbott et al.[1959] 과 Jung et al.[2009] 과비교하였다. 계산은실험치가주어져있는 Re = 3 10 6 에서정상상태로수행하였다. 격자계는 Fig. 8과같이 2차원격자계 (3차원격자계의 xy면 ) 를 z축방향으로확장시켜구성하였다. 날개위격자수는 240개, 전체격자수는약 900,000개로계산영역은날개중심에서전후방, 측방으로코드길이의 30배크기로구성하였다. 계산시간을고려하여날개스팬방향으로중간위치에대칭조건을주었으며, 날개끝의영향을고려하여 z축방향으로날개끝지점으 Table 5. Calculation conditions for hydrofoil Hydrofoil Turbulence model k-ω SST model (2 Eqn.) Velocity (m/s) 6.82 Chord length (m) 0.44 Re 3.0 10 6 y+ 5 Span length(m) 6 Angle of attack(deg) ± 22

조류발전용수직축터빈의유체동력학적 3 차원효과에관한연구 65 Fig. 9. Comparison of lift and drag coefficient. 로부터 1스팬길이만큼여유를두었다. 터빈수치해석시와동일한난류모델을적용하였으며세부조건은 Table 5로정리하였다. Fig. 9는받음각에따른양항력계수를비교한것으로 Re=6 10 6 에서계산을수행한 Jung et al.[2009] 과같이본연구에서도잘일치하는결과를줌을볼수있다. 실속각인 16도정도에근접할수록수치계산결과가실험과차이를보이고있으며, 항력계수가다소크게예측되는등미소한차이를보임에도불구하고전반적으로는수중익의양항력을잘모사하고있다고판단되었다. Fig. 10와 Fig. 11은날개단면의압력계수및한계유선분포를받음각 10도와 22도에서비교한결과이며, Fig. 12는동일한받음각에서 3차원날개표면의압력분포및주위유동을나타낸것이다. 실속이일어나지않는 10도에서날개끝단에형성된날개끝보오텍스패턴만이눈에띄는반면, 22도에서는실속의영향으로날개끝이아닌영역에서도 3차원적인흐름의변형을볼수있다. 이로부터실속이발생하게되면그결과로서날개의 3차원성이증가되 는경향이있음을짐작할수있다. Fig. 13(a) 와 (b) 는각각받음각 10도와 22도에서의 2차원및 3 차원날개표면압력계수를비교한것이다. 3차원날개의압력분포는날개끝단으로부터떨어진거리 t를코드길이로무차원하여표기하였다. 즉, 스팬방향으로날개대칭면지점인 z=0 m부터 1m, 2m, 2.9m, 3m인지점은각각 t/c=6.82, 4.55, 2.27, 0.23, 0에해당한다. 비교적날개중앙부에해당하는 t/c=4.55와 2.27에서의압력계수는 2차원결과와잘일치하고있는반면 t/c=0인날개끝에서는압력면과흡입면의압력차가거의발생하지않으며, t/c=0.23 에서도 2차원에비하여그차이가작아양력이상당부분손실된다는것을알수있다. t/c=2.27은날개끝으로부터전체스팬의 1/4 지점에해당하는데날개중앙에해당하는 t/c=4.55와거의동일한결과를보여주고있다. 이러한경향은받음각 10도와 22도에서공히나타나는현상이다. Fig. 10. Pressure contours (2-D).

66 Fig. 11. Streamlines (2-D). Fig. 12. Pressure contours and streamlines (3-D). Fig. 13. Pressure distributions on blade surface. 현범수 이정기

조류발전용수직축터빈의유체동력학적 3 차원효과에관한연구 67 4. 회전하는터빈날개의 3 차원효과 회전하는터빈날개는직진하는수중익과는달리날개의원주방향위치에따라서매우복잡한유동장이형성된다. Fig. 14는터빈의위상각과속도분포및압력계수를예측한지점을나타낸것이며 t/c=0.114, 6.818에서의결과를 2차원결과와비교하였다. 그림에서원주방향으로 90도위치는교란되지않은유동장에서구동되는반면, 180도위치에서는상류에서전개되어오는터빈후류의영향이지대한영역이기때문에앞장에서살펴본직진하는수중익과그조건이매우상이할것이다. Fig. 15는 TSR=3에서터빈날개가위상각 20, 140, 260도에위치하였을때 2차원터빈과 3차원터빈주위의속도분포를나타낸것이다. Fig. 15(a) 와 Fig. 15(b) 에서는터빈및축주위의속도분포가상하대칭성을보여주지만 Fig. 15(c) 는매우비대칭적인분포를보이고있다. 이와같이터빈날개의스팬방향으로중앙부와끝부분에서의유동이다르기때문에각날개가받는압력분포역시변화가클것이다. 따라서 H/D=1.5 터빈의날개중앙 (t/c=6.818) 과날개끝단 (t/c=0.114) 에서의압력분포를 Fig. 16과같이비교해보았다. 먼저 2차원해석결과를보면, 위상각 140도의압력분포가가장큰양력을생성하며, 20도, 260도순으로양력이감소하는경향이눈에띈다. 위상각 20도는받음각이 0에접근하여양력면으로서의기능을다하지못하고있으며, 위상각 260도에위치한날개는후류의영향으로실질적으로양력을생성하는능력이현저히떨어짐을유추할수있다. 즉, Fig. 16의압력분포들로부터동력계수는위상각 140도날개가주로토크, 즉동력의생성에기여한다고볼수있겠다. 날개의중앙부 (Mid-span) 와끝단의압력분포차이를살펴보자면, 위상각 20도와 140도에서는터빈중앙부와끝단의압력계수가압력면에서상대적으로유사한경향을보이고있는반면, 260도에서는압력면과흡입면공히상당한차이를보이고있음을알수있다. Fig. 15. Velocity contours for 2-D and 3-D turbines. Fig. 14. Coordinates and notations for a turbine blade. 이러한차이는위상각에따른터빈날개의동력계수결과에서명확히확인할수있으며 Fig. 17에나타내었다. Fig. 17(a) 는터빈의한날개의중앙부와끝단의동력계수를비교한것으로 0도에서약 40도까지는동력계수에기여하지않는것을알수있고동력계수가최대가되는위상각은 100도근처임을알수있다. 또한, 후류의영향을받는영역인 180도이후부터도동력에거의전혀기여하지않는것을알수있다. Fig. 17(b) 는 3차원터빈의모든날개

68 현범수 이정기 Fig. 16. Pressure distributions on blade surface. Fig. 18. Comparison of power coefficients between 2-D and 3-D. Fig. 17. Variation of power coefficients with respect to azimuth angle. 에대한중앙부와끝단의동력계수와 2차원터빈동력계수를비교한것으로 2차원과 3차원터빈의날개중앙부동력계수는잘일치하지만끝단의경우최저동력계수차이는없고최고동력계수는 2차원보다약 16% 정도작게예측됨을알수있었다. Fig. 18은 2차원과 3차원터빈의위상각에따른동력계수를비교한것으로서터빈날개 3개의동력계수를합한결과이다. 2차원결 과와비교하였을때동력계수의위상차는없고앞선결과와마찬가지로 3차원해석결과가 2차원결과보다동력계수의최대값부근에서주로차이가발생하고있으며다른곳에서의차이는상대적으로적음을알수있다. 이로부터터빈날개의 3차원성은동력계수의피크부근의성능감소를통하여전체동력계수의저하를야기함을확인할수있었다. 다만이러한경향은터빈의스팬-직경- H/D에따라또다른경향을보일수있기때문에추후보다다양한관점에서재검토해야할것으로판단된다. 참고로 Fig. 19는 TSR=3에서스팬길이만을변화시키면서 2차원터빈의동력계수와스팬-직경비에따른 3차원터빈의평균동력계수의비를스팬-직경비와의관계로나타낸것이다 (Hyun et al. [2012]) 여기서 β=(c P ) 3-D /(C P ) 2-D 이다. 스팬-직경비가 1일때 2차원동력계수대비약 78% 에해당하는결과를보여주며, 스팬-직경비가 6 (AR=34) 일때 91% 정도로수렴하고있다. 한편 Ye 등 (2010) 은스팬-직경비가 3이상인경우자신들의터빈동력계수가 2차원터빈동력계수의 95% 까지근접한다는결론을얻는등 3차원효과를본연구보다적게예측하였는데, 그들의해석이이상유체가정을통한결과였음을감안한다면실제터빈의 3차원효과는본연구결과와같이이상유체해석결과보다다소더크게고려되어야할것

조류발전용수직축터빈의유체동력학적 3 차원효과에관한연구 69 NOMENCLATURE Fig. 19. Relationship between H/D and power coefficient (Hyun et al.[2012]). 으로사료된다. 본연구에서검토한 TSR=3인경우는수직축터빈의효율이우수한구간으로서설계 TSR로삼기에충분하므로이부근에서얻어진결과인 Fig. 19는향후유사한연구를진행함에있어서상당한계산시간이소요되는 3차원터빈의수치계산대신 2차원결과를보정하여사용하는데유용한정보로활용될수있을것이다. 향후다양한 H/D와 TSR에대하여동력계수에미치는 3차원성의해석을추가하여범용 3차원보정계수 (β) 를도출한다면대단히실용적인것으로기대된다. 5. 결론 선행연구를조사 분석하여조류발전용수직축터빈의 3차원효과에대한연구현황을파악하고, 이를검증하여설계에활용하기위하여수치해석을수행하였으며다음의결과들을도출하였다. (1) 3차원효과를파악하기위한기초연구로고정날개수치계산을통해날개중앙부와끝단에서의압력분포차이를확인하였다. (2) 터빈수치계산을통하여 3차원터빈날개의스팬방향으로터빈날개끝단과중앙부에서의속도분포및압력계수차이를살펴본후, 터빈날개끝단에서발생하는양력손실에기인하여 3차원터빈의성능이 2차원터빈성능보다저하됨을확인하였다. (3) 본연구의설계속도인 TSR=3에서는터빈의스팬-직경비가 3이상으로증가함에따라 2차원터빈의효율로수렴해감을알수있었다. 하지만공진회피설계를위해서는스팬-직경비가작은것이유리하기때문에실질적인수직축터빈의 3차원효과는 2차원해석결과대비 15% 이상의차이를감수해야할것이다. (4) 보다다양한작동조건에대한수치해석결과가보완된다면복잡한 3차원해석대신간편한 2차원해석을보정하는방식으로도수직축터빈의성능을도출하는실용적인해석이가능하다고사료된다. AR = aspect ratio, AR=H/c c = chord length C P = power coefficient (C P =Tω/0.5ρV 3 DH) D = turbine diameter (D=2R) H = span length Re = Reynolds number T = torque TSR = tip speed ratio, TSR=Rω/V V = free stream velocity Z = number of blade α = angle of attack β = Ratio between 3-D and 2-D power coefficients (β=(c P ) 3-D /(C P ) 2-D ) ρ = fluid density σ = solidity, σ=zc/πd ω = rotational speed 참고문헌 [1]Abbott, H. I et al., 1959, Theory of Wing Sections including a summary of airfoil data, 632-633. [2]Calcagno, G. and Moroso, A., 2006, Experiment and Numerical investigation of an Innovative Technology for Marine Current Exploitation: the Kobold Turbine, Proc. of 16th International Offshore and Polar Engineering Conference, San Francisco, California, USA, May 28-June. [3]Castelli, M. R. et al., 2011, The Darrieus wind turbine : Proposal for a new performance prediction model based on CFD, Energy, Vol.36, 4919-4934. [4]Han, J.S., Choi, D.H., Hyun, B.S., Kim, M.C., Rhee, S.H. and Song, M.S., 2011, Parametric Numerical Study on the Performance of Helical Tidal Stream Turbines, J. of the Korean Society for Marine Environmental Engineering, Vol.14, No.2, 114-120. [5]Han, J.S., Hyun, B.S., Choi, D.H., Mo, J.O., Kim, M.C. and Rhee, S.H., 2010, Study on Design of Darrieus-type Tidal Stream Turbine Using Parametric Study, J. of the Korean Society for Marine Environmental Engineering, Vol.13, No.4, 241-248. [6]Howell, R. et al., 2010, Wind tunnel and numerical study of a small vertical axis wind turbine, Renewable Energy, Vol.35, 412-422. [7]Hyun, BS et al., 2012, Development of Standardized S/W System for the Integrated Design of Tidal Current Turbines, Final Report, Tidal Stream Energy Research Center, Busan, Korea. [8]Jung, H.J., Lee, J.H., Rhee, S.H., Song, M.S. and Hyun, B.S., 2009, Unsteady Flow Around a Two-dimensional Section of a Vertical Axis Turbine for Tidal Stream Energy Conversion, Int l J. Naval Architecture & Ocean Engineering, Vol.1, 64-69. [9]Kim, B.S., Kim, M.E. and Lee, Y.H., 2008, Basic Configura-

70 현범수 이정기 tion Design and Performance Analysis of a 100kW Wind Turbine Blade using Blade Element Momentum Theory, J. of the Korean Society of Marine Engineering, v.32, no.6, 827-833 [10] Laneville, A. and Vittecoq, P., 1986, Dynamic Stall: The Case of the Vertical Axis Wind Turbine, Journal of Solar Energy Engineering, Vol.108, 140-145. [11] Lee, D.H., Hyun, B.S., Lee, J.K., Kim, M.C. and Rhee, S.H., 2012, Development of Hydrodynamic Capacity Evaluation Method for a Vertical-Axis Tidal Stream Turbine, J. of the Korean Society for Marine Environmental Engineering, Vol.15, No.2, 142-149. [12] Li, Y. and Calisal, S. M., 2010, Three-dimensional effects and arm effects on modeling a vertical axis tidal current turbine, Renewable Energy, Vol.35, 2325-2334. [13] Mitsuhiro, S. et al., 2000, An Experimental Study of the Characteristics of a Darrieus Turbine for Tidal Power Generation, J. of Electrical Engineering in Japan, Vol.132 (Translated from Denki Gakkai Ronbunshi, Vol.118-B, 781-787, 1988). [14] Paraschivoiu, I. et al., 1988, Blade Tip, Finite Aspect Ratio, and Dynamic Stall Effects on the Darrieus Rotor, J. of Propulsion, Vol.4, No.1, 73-80. [15] Paraschivoiu, I., 2002, Wind Turbine Design with Emphasis on Darrieus Concept, 177-188. [16] Sheldahl, R et al., 1981, Aerodynamic characteristics of seven symmetrical airfoil sections through 180-degree angle of attack for use in aerodynamic analysis of vertical axis wind turbines, Technical Report SAND 80-2114, Sandia National Labs., Alberquerque, USA. [17] Winchester, J. D et al., 2009, Torque ripple and variable blade force : A comparison of Darrieus and Gorlov-type turbines for tidal stream energy conversion, Proc. of 8th European Wave and Tidal Energy Conference, Uppsala, Sweden. 2012 년 10 월 15 일원고접수 2013 년 2 월 4 일 (1 차 ), 5 월 7 일 (2 차 ) 심사수정일자 2013 년 5 월 13 일게재확정일자