다이아몬드 필러가 적용된 통풍식 디스크 브레이크 로터의 공기 흐름 및 열전달 분석

Analysis of Air Flow and Heat Transfer in Ventilated Disc Brake Rotor with Diamond Pillars

브레이크 작동 중 발생하는 열은 디스크에 저장되며, 이 열은 주로 강제 대류를 통해 주변으로 방출됩니다. 통풍식 브레이크 로터의 효과적인 설계는 제동 시스템의 적절한 냉각을 달성하는 데 필수적입니다. 기존의 방사형 베인 로터는 베인 통로 내 온도 분포가 균일하지 않아 높은 열 응력을 유발하고, 이는 결국 열적 파손으로 이어질 수 있습니다. 본 연구는 방사형 베인의 대안으로 온도 분포가 더 균일한 다이아몬드 필러형 로터의 공기 흐름 및 열전달 특성을 이해하는 데 중점을 둡니다. CFD 코드를 사용하여 다이아몬드 필러 로터 통로를 통한 대류 열 방산을 분석하고 개선하였습니다. 특히 필러의 경사각 변화가 재순환 영역을 줄이고 질량 유량 및 열전달을 향상시키는 효과를 조사했습니다. 이 연구는 고성능 차량의 제동 시스템 설계 및 열 관리 최적화에 중요한 기여를 합니다.

메타데이터 및 키워드

Fig.2 Diamond Pillared Brake Rotor
Fig.2 Diamond Pillared Brake Rotor

논문 메타데이터

  • Industry: 자동차 (Automotive)
  • Material: 명시되지 않음 (필러형 로터의 균일한 재료 분포 언급)
  • Process: 전산유체역학(CFD) 시뮬레이션 및 실험적 검증
  • System: 다이아몬드 필러가 적용된 통풍식 디스크 브레이크 로터
  • Objective: 필러 방향 수정을 통해 다이아몬드 필러 브레이크 로터의 대류 열전달 계수 분석 및 향상

핵심 키워드

  • 브레이크 디스크
  • 통풍식
  • 열적 파손
  • 열 응력
  • 다이아몬드 필러 로터
  • CFD

핵심 요약

연구 구조

본 연구는 다이아몬드 필러가 있는 통풍식 브레이크 로터의 3D 모델링을 수행하고, ICEM-CFD와 ANSYS Fluent를 사용하여 공기 흐름 및 열전달 특성을 시뮬레이션했습니다.

방법 개요

800 rpm의 회전 속도와 700 K의 로터 벽면 온도 조건에서 시뮬레이션을 수행했으며, 테이퍼형 방사형 베인(TRV) 로터에 대한 실험 데이터와 비교하여 CFD 모델의 타당성을 검증했습니다.

주요 결과

다이아몬드 필러를 회전 방향으로 20도 기울였을 때, 대류 열전달 계수가 기본 설계 대비 11% 향상(42.16 W/m²K에서 46.80 W/m²K로 증가)되었으며, 공기 질량 유량은 16.11% 증가했습니다.

산업적 활용 가능성

고성능 자동차 제동 시스템 설계, 경주용 차량의 로터 최적화, 대형 차량의 열 관리 시스템 등에 적용 가능합니다.

한계와 유의점

CFD 시뮬레이션에서 복사 및 부력 효과는 무시되었으며, 분석은 800 rpm의 일정한 회전 속도로 제한되었습니다.

논문 상세 정보

1. 개요

  • Title: Analysis of Air Flow and Heat Transfer in Ventilated Disc Brake Rotor with Diamond Pillars
  • Author: Gorakh B. Kudal and Mahesh R. Chopade
  • Year: 2016
  • Journal: International Journal of Current Engineering and Technology
  • DOI/Link: 논문에 명시되지 않음

2. 초록

제동 작동 중 발생하는 열은 브레이크 디스크에 저장되며, 이 열은 주로 강제 대류에 의해 주변으로 방출됩니다.

제동 시스템의 적절한 냉각을 달성하기 위해서는 통풍식 브레이크 로터의 효과적인 설계가 필수적입니다.

브레이크 디스크의 열적 파손은 고온뿐만 아니라 로터 통로 내에서 발생하는 높은 열 응력으로 인해서도 발생한다는 것이 밝혀졌습니다.

필러형 로터는 로터 통로 내에서 더 균일한 온도 분포를 가지므로 방사형 베인 로터의 대안이 될 수 있습니다.

CFD 코드를 사용하여 다이아몬드 필러 브레이크 로터 통로를 통한 대류 열 방산을 분석하고 개선하였습니다.

수정된 다이아몬드 필러 브레이크 로터 구성은 기본 설계와 비교하여 대류 열전달 계수에서 11%의 개선을 보여줍니다.

3. 방법론

전산 모델링 및 시뮬레이션: 격자 생성을 위해 ICEM-CFD를 사용하고 시뮬레이션을 위해 ANSYS Fluent를 사용했습니다. 회전 대칭성을 고려하여 36개 통로 중 20도 세그먼트(2개 통로)를 분석 대상으로 삼았으며, IGS 형식의 3D 모델을 임포트하여 분석을 진행했습니다.

검증 및 격자 독립성: 테이퍼형 방사형 베인(TRV) 로터에 대해 실험적 검증을 수행했습니다. 로터를 100°C로 가열하고 500 rpm으로 회전시키면서 온도 강하를 기록했습니다. 실험값(37.41 W/m²K)과 CFD 결과(35.80 W/m²K) 사이의 편차는 약 4.30%로 양호한 일치를 보였습니다.

격자 독립성 연구: 결과가 메시 밀도에 의존하지 않도록 224,000개에서 581,000개 사이의 격자 크기를 테스트했습니다. 432,000개 이상의 격자에서 대류 열전달 계수의 변화가 미미함을 확인하고 이를 최종 시뮬레이션 격자로 결정했습니다.

4. 결과 및 분석

기본 설계 분석: 다이아몬드 필러의 상단부에서 공기 질량 유량과 열전달을 감소시키는 재순환 영역을 확인했습니다. 기본 설계의 대류 열전달 계수는 42.16 W/m²K, 질량 유량은 1.08 g/s로 측정되었으며, 재순환 영역은 다른 영역에 비해 매우 낮은 열전달율을 보였습니다.

필러 경사각의 효과: 모든 다이아몬드 필러를 회전 방향(시계 방향)으로 기울이면 재순환 영역이 줄어들고 성능이 향상되었습니다. 5도(4.65%), 10도(6.78%), 15도(7.61%), 20도(11.00%) 순으로 열전달 계수가 개선되었으며, 20도 경사 시 질량 유량은 1.254 g/s로 16.11% 증가했습니다.

Fig.5 Convective heat transfer coefficients against grid
Fig.5 Convective heat transfer coefficients against grid
Figure 7: 수정된 다이아몬드 필러 구성 결과. 경사각(5°~20°) 증가에 따른 유동 재정렬 및 재순환 영역의 감소를 시각적으로 보여줍니다.
Figure 7: 수정된 다이아몬드 필러 구성 결과. 경사각(5°~20°) 증가에 따른 유동 재정렬 및 재순환 영역의 감소를 시각적으로 보여줍니다.

5. 그림 및 표 목록 (Figure and Table List)

  • Figure 1: 디스크 브레이크 로터 기하학적 구조 (치수 단위: mm). 수치 모델에 사용된 물리적 치수(R145, R90, R76, R85 등)를 제공합니다.
  • Figure 4: 경계 조건이 포함된 CFD 모델. 20도 세그먼트와 개방형 경계 조건을 포함한 계산 영역을 시각화합니다.
  • Figure 6: 기본 다이아몬드 필러 브레이크 로터 결과. 온도 등고선, 유선 및 재순환 영역을 식별하는 속도 벡터를 보여줍니다.
  • Figure 7: 수정된 다이아몬드 필러 구성 결과. 경사각(5°~20°) 증가에 따른 유동 재정렬 및 재순환 영역의 감소를 시각적으로 보여줍니다.
  • Table 1: 수정된 다이아몬드 필러 브레이크 로터 설계의 800 rpm에서의 대류 열전달 계수 및 질량 유량. 각 경사각에 따른 정량적 개선 수치를 요약하여 제공합니다.

6. 참고문헌

  • Lisa Wallis et al. (2002). Air Flow and Heat Transfer in Ventilated Disc Brake Rotors with Diamond and Tear-Drop Pillars. Numerical Heat Transfer, Part A, 41:643-655.
  • Parish D., MacManus D. G. (2005). Aerodynamic Investigations of Ventilated Brake Discs. Proceedings of Institution of Mechanical Engineers, Part D: Journal of Automobile Engineering, 219, p.471-486.
  • Reddy S., Mallikarjuna J., Ganesan V. (2008). Flow and Heat Transfer Analysis of a Ventilated Disc Brake Rotor Using CFD. SAE Technical Paper 2008-01-0822.

기술 Q&A (Technical Q&A)

Q: 방사형 베인 로터 대신 필러형 로터를 선택한 주요 이유는 무엇입니까?

필러형 로터는 방사형 베인 로터에 비해 재료 분포가 더 균일하여 균열 전파에 대한 저항력이 더 높습니다. 또한, 로터 통로 내에서 온도 분포를 더 균일하게 형성하여 열 응력을 유의미하게 낮출 수 있기 때문에 열적 파손 방지에 유리합니다.

Q: 20도 필러 경사각에서 대류 열전달 계수는 얼마나 개선되었습니까?

연구 결과에 따르면, 다이아몬드 필러를 회전 방향으로 20도 기울였을 때 대류 열전달 계수는 기본 설계의 42.16 W/m²K에서 46.80 W/m²K로 약 11.00% 향상되었습니다. 이는 공기 흐름의 최적화를 통해 냉각 효율이 직접적으로 개선되었음을 의미합니다.

Q: 재순환 영역이 로터의 열 성능에 어떤 영향을 미칩니까?

재순환 영역은 로터 통로를 통과하는 공기의 질량 유량을 감소시키고 공기를 정체하게 만듭니다. 이로 인해 해당 영역의 대류 열전달율이 매우 낮아지며, 국부적인 온도 상승을 유발하여 핫스팟을 형성하고 열 응력을 증가시키는 부정적인 영향을 미칩니다.

Q: CFD 모델의 타당성 검증은 어떻게 수행되었습니까?

테이퍼형 방사형 베인(TRV) 로터를 사용하여 실험적 검증을 수행했습니다. 1.0 kW 히터로 로터를 가열한 후 500 rpm으로 회전시키며 온도 강하를 측정했으며, 실험을 통해 얻은 열전달 계수(37.41 W/m²K)와 CFD 결과(35.80 W/m²K)의 오차가 4.30%로 나타나 모델의 신뢰성을 확보했습니다.

Q: 격자 독립성 연구를 통해 결정된 최적의 격자 수는 얼마입니까?

224,000개부터 581,000개까지 다양한 격자 크기를 테스트한 결과, 432,000개 이상의 격자에서는 대류 열전달 계수의 변화가 무시할 수 있는 수준으로 나타났습니다. 따라서 계산 효율성과 정확성을 모두 고려하여 432,000개의 셀을 최종 분석용 격자로 결정했습니다.

Q: 필러 경사각 변화에 따른 질량 유량의 개선 경향은 어떠합니까?

필러의 경사각이 커질수록 공기 질량 유량은 점진적으로 증가하는 경향을 보였습니다. 구체적으로 20도 경사각에서 질량 유량은 1.254 g/s를 기록하여 기본 설계(1.08 g/s) 대비 약 16.11% 증가했으며, 이는 필러의 기울기가 공기 유입을 더 원활하게 유도함을 보여줍니다.

결론

본 연구는 다이아몬드 필러 로터 내의 재순환 영역이 열전달을 저해하는 주요 요인임을 확인했습니다. 필러를 회전 방향으로 20도 기울임으로써 이러한 재순환 영역을 효과적으로 줄였고, 결과적으로 대류 열전달 계수 11% 향상과 공기 질량 유량 16% 증가라는 성과를 거두었습니다.

이러한 결과는 브레이크 로터의 열적 파손을 방지하고 냉각 효율을 극대화하기 위한 기하학적 최적화의 중요성을 시사합니다. 본 연구에서 제시된 설계 변경은 고성능 차량의 제동 시스템 안정성을 높이는 데 기여할 수 있으며, 향후 다양한 운전 조건에서의 추가 검증이 기대됩니다.

출처 정보 (Source Information)

Citation: Gorakh B. Kudal and Mahesh R. Chopade (2016). Analysis of Air Flow and Heat Transfer in Ventilated Disc Brake Rotor with Diamond Pillars. International Journal of Current Engineering and Technology.

DOI/Link: 논문에 명시되지 않음

Technical Review Resources for Engineers:

▶ 논문에 명시되지 않음
▶ 기술 검토 및 적용 가능성 문의

This material is for informational purposes only. Unauthorized commercial use is prohibited.

Copyright © 2026 STI C&D. All rights reserved.

CFD를 이용한 두 가지 레이놀즈 수에서의 NACA0012 및 NACA0018 에어포일 난류 운동 에너지 분석

Turbulent Kinetic Energy Analysis of NACA0012 and NACA0018 Airfoils at Two Reynolds Number Using CFD Tool

에어포일의 설계와 성능을 최적화하기 위해서는 공기역학적 특성에 대한 정밀한 연구가 필수적입니다. 본 연구는 대칭형 에어포일인 NACA 0012와 NACA 0018을 대상으로, 레이놀즈 수 300,000 및 700,000 조건에서 발생하는 난류 운동 에너지(TKE) 분포를 전산유체역학(CFD) 시뮬레이션을 통해 분석하였습니다. 난류 운동 에너지는 항공기의 항력 및 소음 발생과 직결되는 요소로, 이를 제어하는 것은 공학적 설계의 핵심 과제입니다. 본 논문은 에어포일의 두께와 유속이 난류 생성에 미치는 영향을 정량적으로 비교하여 제시합니다. 특히 실속(stall) 발생 이후의 유동 박리와 난류 분포 변화를 시각화하여 분석의 신뢰도를 높였습니다. 연구 결과는 항공기 날개 설계와 풍력 터빈 블레이드 설계에 있어 에어포일 선택의 기준을 제공합니다. 이러한 분석은 저레이놀즈 수 영역에서의 유동 특성을 이해하고 시스템 효율을 극대화하는 데 중요한 기여를 합니다. 본 연구는 특정 운영 조건에 따른 최적의 에어포일 형상을 제안함으로써 실무적인 설계 가이드를 제공하는 데 목적이 있습니다.

메타데이터 및 키워드

Fig.3 NACA 0018 airfoil meshing domain
Fig.3 NACA 0018 airfoil meshing domain

논문 메타데이터

  • Industry: 항공우주, 신재생 에너지
  • Material: 공기 (Air)
  • Process: 전산유체역학 (CFD) 시뮬레이션
  • System: NACA 0012 및 NACA 0018 에어포일
  • Objective: 다양한 레이놀즈 수와 받음각에서 두 대칭 에어포일의 난류 운동 에너지(TKE) 분포 분석 및 비교

핵심 키워드

  • NACA
  • 레이놀즈 수 (Reynolds Number)
  • 난류 운동 에너지 (Turbulence Kinetic Energy)
  • 받음각 (Angle of Attack)
  • 속도 (Velocity)
  • 전산유체역학 (CFD)
  • 유동 박리 (Flow Separation)

핵심 요약

연구 구조

본 연구는 Ansys FLUENT 14.5를 사용하여 NACA 0012(두께 12%)와 NACA 0018(두께 18%) 에어포일의 2차원 유동 시뮬레이션을 수행하였습니다. SST k-omega 난류 모델을 적용하여 받음각 변화에 따른 난류 특성을 분석했습니다.

방법 개요

1m 코드 길이를 가진 에어포일 모델에 대해 약 89,000개의 셀로 구성된 비정렬 격자를 생성하였으며, 레이놀즈 수 300,000(유속 4.38 m/s)과 700,000(유속 10.23 m/s) 조건에서 정상 상태(steady-state) 해석을 진행했습니다.

주요 결과

NACA 0018 에어포일은 NACA 0012보다 유의미하게 높은 난류를 생성하며, 이는 두께 증가에 따른 공기 입자와의 마찰 증가에 기인합니다. 실속각은 NACA 0012의 경우 약 17도, NACA 0018의 경우 약 16도로 나타났습니다. 레이놀즈 수가 300,000에서 700,000으로 증가함에 따라 TKE 분포의 ‘레드 존(red zone)’이 확장되는 것이 확인되었습니다.

산업적 활용 가능성

난류 최소화가 중요한 항공기 날개 설계에는 NACA 0012가 권장되며, 두께가 중요하고 난류 영향이 상대적으로 적은 수평/수직축 풍력 터빈에는 NACA 0018이 적합합니다.

한계와 유의점

본 연구는 2차원 시뮬레이션에 국한되어 3차원 유동 효과를 완전히 반영하지 못할 수 있으며, 실속 이후의 복잡한 공기역학적 특성에 대해서는 추가적인 연구가 필요합니다.

논문 상세 정보

1. 개요

  • Title: Turbulent Kinetic Energy Analysis of NACA0012 and NACA0018 Airfoils at Two Reynolds Number Using CFD Tool
  • Author: Md. Abdus Shabur, Asif Hasan Khan, Mohammad Ali, SM Fahim Faisal
  • Year: 2022
  • Journal: ScienceOpen Preprints
  • DOI/Link: https://doi.org/10.14293/S2199-1006.1.SOR-.PPQGB0H.v1

2. 초록

에어포일의 설계와 성능을 개선하기 위해서는 공기역학적 특성을 적절히 연구해야 합니다.

본 연구에서는 두 가지 대칭 에어포일(NACA 0018 및 NACA 0012)의 난류 운동 에너지를 300,000과 700,000의 두 가지 레이놀즈 수에서 조사하였습니다.

본 논문에 제시된 전산 시뮬레이션 작업은 몇 가지 흥미로운 결과를 도출했습니다.

본 논문에서는 다양한 받음각과 두 가지 다른 기류 속도 하에서 두 에어포일에 걸친 난류 운동 에너지의 분포를 설명합니다.

유동 박리로 인해 실속각 이후 상당한 양의 난류가 관찰되었습니다.

최종 결과에 따르면 NACA 0018 에어포일은 NACA 0012 에어포일보다 훨씬 더 많은 난류를 생성합니다.

이는 에어포일의 두께가 증가함에 따라 공기 입자와의 마찰이 더 많이 발생하기 때문입니다.

레이놀즈 수 또는 기류 속도 측면에서의 두 번째 비교에서, 본 논문의 레드 존으로 표시된 바와 같이 높은 속도는 더 많은 난류를 생성합니다.

항공기에서는 난류가 항상 가능한 한 낮을 것으로 기대되므로, 항공기 응용 분야에는 NACA 0012 에어포일을 사용하고, 난류의 중요성이 덜한 수평 또는 수직축 풍력 터빈 응용 분야에는 NACA 0018 에어포일을 사용할 것을 권장합니다.

3. 방법론

CFD 소프트웨어 및 난류 모델: Ansys FLUENT 14.5를 사용하여 계산을 수행하였으며, 양력(CL) 및 항력(CD) 변동을 정밀하게 포착하기 위해 SST k-omega 난류 모델을 채택했습니다. 이 모델은 벽면 근처의 유동 해석에 강점이 있어 에어포일 표면의 난류 운동 에너지 계산에 적합합니다. 압력 기반 솔버와 절대 속도 공식을 사용하였으며, 2차원 평면 공간에서 정상 상태(steady-state) 해석을 수행했습니다. 솔버의 세부 설정 매개변수는 논문의 Table 1에 명시된 기준을 따랐습니다.

에어포일 기하학 및 격자 생성: 연구 대상은 1m 코드 길이를 가진 NACA 0012(12% 두께)와 NACA 0018(18% 두께) 대칭 에어포일입니다. 수치 해석의 정확도를 보장하기 위해 약 89,000개의 셀로 구성된 2차원 비정렬 격자를 생성했습니다. 격자 구조는 ‘C’형 도메인을 채택하였으며, 도메인 반경은 12.5m로 설정하여 경계 조건의 영향을 최소화했습니다. NACA 0012와 NACA 0018의 형상적 차이에 따라 면적을 정밀하게 계산하여 격자 밀도를 최적화했습니다.

경계 조건 설정: 입구 경계 조건(Velocity Inlet)은 레이놀즈 수에 맞춰 설정되었습니다. Re=300,000 조건에서는 4.38 m/s, Re=700,000 조건에서는 10.23 m/s의 유속을 적용했습니다. 출구(Pressure Outlet)는 0 Pa의 게이지 압력으로 설정하였으며, 에어포일 표면에는 점착 조건(no-slip wall condition)을 부여하여 벽면에서의 속도가 0이 되도록 했습니다. 이러한 설정은 실제 풍동 실험 환경을 수치적으로 모사하기 위한 것입니다.

수치 해석 기법: 압력-속도 커플링을 위해 Coupled 기법을 사용하였으며, 공간 이산화에는 구배 계산을 위한 Least Squares Cell Based 방식을 적용했습니다. 압력과 운동량 항에는 2차 상향풍(second-order upwind) 방식을, 난류 운동 에너지 및 소산율에는 1차 상향풍 방식을 적용하여 수렴성과 정확도를 동시에 확보했습니다. CFL 수는 0.9로 설정되었으며, 압력 및 운동량의 명시적 이완 계수는 0.75를 유지했습니다.

4. 결과 및 분석

수치 모델 검증: 본 연구의 시뮬레이션 결과인 양력 계수(CL)와 항력 계수(CD)를 Timmer(2008) 및 Jacobs(1937)의 실험 데이터와 비교하여 검증을 수행했습니다. Re=300,000, 받음각 10°에서 NACA 0018의 CL은 0.852로 계산되었으며, 이는 Jacob의 실험값 0.85와 매우 유사한 결과입니다. CD 값 또한 실험 데이터와 실질적으로 동일한 경향을 보여, 실속각 이전까지의 수치 해석 모델이 매우 높은 신뢰도를 가짐을 입증했습니다.

난류 운동 에너지(TKE) 분포 분석: TKE는 낮은 받음각에서 주로 에어포일의 뒷전(trailing edge) 부근에 집중되지만, 받음각이 증가함에 따라 에어포일 상부 표면 전체로 확산되는 양상을 보입니다. NACA 0012의 실속각은 약 17도, NACA 0018은 약 16도로 나타났으며, 실속 이후 유동 박리로 인해 TKE 수치가 급격히 상승했습니다. 특히 두께가 18%인 NACA 0018은 12%인 NACA 0012보다 훨씬 더 넓고 강한 난류 구역을 형성하는 것이 확인되었습니다.

레이놀즈 수에 따른 영향 비교: 레이놀즈 수(기류 속도)의 증가는 난류 운동 에너지의 직접적인 상승을 초래합니다. Re=300,000과 Re=700,000 조건을 비교 분석한 결과, 유속이 높은 Re=700,000 조건에서 TKE 분포 컨투어의 ‘레드 존’이 훨씬 더 광범위하게 나타났습니다. 이는 고속 유동에서 관성력이 지배적이 됨에 따라 난류 강도가 강화됨을 의미하며, 고속 비행체 설계 시 난류 제어의 중요성을 시사하는 결과입니다.

Fig.4 NACA0018 airfoil meshing (enlarged view)
Fig.4 NACA0018 airfoil meshing (enlarged view)
Fig.23 TKE Distribution at 17 degree AOA
Fig.23 TKE Distribution at 17 degree AOA
Fig.24 TKE Distribution at 20 degree AOA
Fig.24 TKE Distribution at 20 degree AOA

5. 그림 및 표 목록 (Figure and Table List)

  • Figure 1-2: NACA 0012 및 NACA 0018 에어포일 형상. 연구에 사용된 두 대칭 에어포일의 기하학적 프로파일을 보여줍니다.
  • Table 1: 솔버 설정 매개변수. Ansys FLUENT에서 사용된 CFD 솔버 구성을 상세히 정의합니다.
  • Table 2: Re=300,000에서 NACA 0018의 CL 및 CD 결과와 실험 결과 비교. 수치 모델의 정량적 타당성을 검증하는 데이터입니다.
  • Figure 7-12: Re=300,000에서 NACA 0012의 TKE 분포 (받음각 0~20도). 받음각 증가에 따른 난류 및 유동 박리의 진행 과정을 시각화합니다.
  • Figure 31: 받음각 20도에서의 TKE 분포 (NACA 0018, Re=700,000). 고속 및 고두께 조건에서 관찰된 최대 난류 상태를 보여주는 핵심 데이터입니다.

6. 참고문헌

  • M. S. Howe, G. C. Lauchle, and J. Wang. (2001). Aerodynamic lift and drag fluctuations of a sphere. J. Fluid Mech. 10.1017/S0022112001003925
  • A. Shabur, A. Hasan, and M. Ali. (2021). Comparison of Aerodynamic Behaviour between NACA 0018 and NACA 0012 Airfoils at Low Reynolds Number Through CFD Analysis. 10.5281/zenodo.4003677
  • W. A. Timmer. (2008). Summary of the Delft University Wind Turbine Dedicated Airfoils. J. Sol. Energy Eng. https://asmedigitalcollection.asme.org/solarenergyengineering/article-abstract/125/4/488/464752/Summary-of-the-Delft-University-Wind-Turbine

기술 Q&A (Technical Q&A)

Q: 에어포일의 두께가 난류 운동 에너지 생성에 미치는 영향은 무엇입니까?

본 연구 결과에 따르면 에어포일의 두께가 증가할수록(NACA 0012에서 0018로) 공기 입자와의 마찰이 증가하여 난류 운동 에너지(TKE) 생성이 유의미하게 높아집니다. NACA 0018은 18%의 두께를 가져 12%인 NACA 0012보다 더 넓은 난류 구역을 형성합니다. 이는 두꺼운 에어포일이 유동 박리를 더 일찍 유도하거나 강하게 만들기 때문입니다. 따라서 난류 억제가 중요한 설계에서는 얇은 에어포일이 유리합니다.

Q: 기류 속도(레이놀즈 수)와 TKE 분포 사이에는 어떤 상관관계가 있습니까?

기류 속도가 빨라질수록, 즉 레이놀즈 수가 높아질수록 난류 운동 에너지는 증가합니다. 연구에서 Re=300,000과 Re=700,000을 비교했을 때, 높은 레이놀즈 수 조건에서 TKE 컨투어의 ‘레드 존’이 더 크게 나타났습니다. 이는 유속 증가가 유동 내의 관성력을 높여 난류 강도를 강화시키기 때문입니다. 결과적으로 고속 운행 시 난류로 인한 에너지 손실이 더 커짐을 의미합니다.

Q: 항공기 응용 분야에서 NACA 0012가 권장되는 이유는 무엇입니까?

항공기 설계에서는 항력을 줄이고 연료 효율을 높이기 위해 난류 발생을 최소화하는 것이 필수적입니다. 본 연구에서 NACA 0012는 NACA 0018에 비해 상대적으로 낮은 난류 운동 에너지를 생성하는 것으로 입증되었습니다. 낮은 TKE는 더 매끄러운 유동 유지를 의미하며, 이는 항공기의 안정성과 성능 향상에 직결됩니다. 따라서 저난류 특성을 가진 NACA 0012가 항공기 날개에 더 적합합니다.

Q: 풍력 터빈 설계에서 NACA 0018이 사용될 수 있는 이유는 무엇입니까?

풍력 터빈 블레이드 설계에서는 구조적 강도를 위해 일정 수준 이상의 두께가 필요합니다. NACA 0018은 NACA 0012보다 두꺼워 구조적으로 더 견고한 블레이드 제작이 가능합니다. 비록 난류 발생량은 더 많지만, 풍력 터빈은 항공기에 비해 난류로 인한 성능 저하가 상대적으로 덜 치명적일 수 있는 환경에서 작동합니다. 따라서 구조적 이점과 공기역학적 효율 사이의 타협점으로 NACA 0018이 선택될 수 있습니다.

Q: 실속(Stall) 발생 이후 난류 분포는 어떻게 변화합니까?

받음각이 실속각(NACA 0012의 경우 약 17도)을 넘어서면 에어포일 상부 표면에서 급격한 유동 박리가 발생합니다. 이 시점부터 난류 운동 에너지는 에어포일의 뒷전뿐만 아니라 윗면 전체로 광범위하게 확산됩니다. 시뮬레이션 결과, 실속 이후의 TKE 값은 급격히 상승하며 매우 불규칙한 유동 패턴을 보입니다. 이는 양력의 급격한 감소와 항력의 증가를 초래하는 주요 원인이 됩니다.

Q: 본 연구에서 사용된 CFD 모델의 신뢰성은 어떻게 검증되었습니까?

연구팀은 시뮬레이션에서 얻은 양력 계수(CL)와 항력 계수(CD)를 기존의 신뢰도 높은 실험 데이터(Timmer, 2008 및 Jacobs, 1937)와 비교했습니다. 예를 들어, Re=300,000 조건의 받음각 10도에서 계산된 CL 값은 0.852로 실험값인 0.85와 매우 근접했습니다. 이러한 정량적 비교를 통해 실속각 이전까지의 수치 해석 결과가 실험 데이터와 실질적으로 일치함을 확인하여 모델의 타당성을 확보했습니다.

결론

본 연구는 CFD 시뮬레이션을 통해 에어포일의 두께와 레이놀즈 수가 난류 운동 에너지(TKE) 분포에 미치는 영향을 성공적으로 규명하였습니다. 분석 결과, 에어포일의 두께가 두꺼워질수록 공기 마찰 증가로 인해 난류 생성이 활발해지며, 유속이 빨라질수록 난류의 강도가 강화됨을 확인하였습니다. 특히 NACA 0012와 NACA 0018의 비교를 통해 각 에어포일의 특성에 맞는 최적의 응용 분야를 제시한 점이 주요 성과입니다.

공학적 관점에서, 난류 최소화가 최우선인 항공기 설계에는 NACA 0012와 같은 얇은 에어포일이 유리하며, 구조적 강성이 요구되는 풍력 터빈 블레이드에는 NACA 0018이 적합한 대안이 될 수 있습니다. 다만, 본 연구는 2차원 해석에 기반하고 있으므로 실제 3차원 환경에서의 유동 간섭이나 실속 이후의 복잡한 거동을 완벽히 예측하기에는 한계가 있습니다. 향후 연구에서는 3차원 효과와 다양한 레이놀즈 수 범위에 대한 추가적인 검증이 이루어져야 할 것입니다.

출처 정보 (Source Information)

Citation: Md. Abdus Shabur, Asif Hasan Khan, Mohammad Ali, SM Fahim Faisal (2022). Turbulent Kinetic Energy Analysis of NACA0012 and NACA0018 Airfoils at Two Reynolds Number Using CFD Tool. ScienceOpen Preprints.

DOI/Link: https://doi.org/10.14293/S2199-1006.1.SOR-.PPQGB0H.v1

Technical Review Resources for Engineers:

▶ 원문 논문 보기 (PDF)
▶ 기술 검토 및 적용 가능성 문의

This material is for informational purposes only. Unauthorized commercial use is prohibited.

Copyright © 2026 STI C&D. All rights reserved.

Fig. 2: Variations in Mach number along the axis of Bell nozzle for cold flow.

벨 노즐의 전산 해석

벨 노즐의 전산 해석

Computational Analysis of Bell Nozzles

본 보고서는 로켓 엔진에서 널리 사용되는 벨형 노즐의 설계 및 비설계 조건에서의 유동 특성을 전산 유체 역학(CFD)을 통해 분석한 연구 결과를 담고 있습니다. 연구의 핵심은 특성곡선법(MOC)을 이용한 노즐 최적화와 복잡한 충격파 구조의 수치적 규명에 있습니다.

Paper Metadata

  • Industry: 항공우주 (Aerospace)
  • Material: 티타늄 (Titanium), 공기 (Air), 수소 (Hydrogen)
  • Process: 전산 유체 역학 (CFD), 특성곡선법 (MOC)

Keywords

  • 벨형 노즐 (Bell type Nozzle)
  • 수치 해석 (Numerical Analysis)
  • 압축파 (Compressive Waves)
  • 충격파 (Shock Wave)
  • 특성곡선법 (Method of Characteristics)
  • 유동 박리 (Flow Separation)

Executive Summary

Research Architecture

본 연구는 벨형 노즐 내부의 2차원 축대칭 유동을 분석하기 위해 수치적 접근 방식을 채택하였습니다. C 언어로 작성된 특성곡선법(MOC) 및 유선 함수 코드를 사용하여 고효율 노즐 윤곽을 설계하였으며, GAMBIT 2.4.6을 통해 격자를 생성하고 FLUENT 6.3.26을 사용하여 시뮬레이션을 수행하였습니다. 해석은 냉각 유동(300K 공기)과 가열 유동(1000K 수소)의 두 가지 조건에서 설계 압력 및 다양한 비설계 배압(5, 7, 11, 15 bar)을 대상으로 진행되었습니다.

Fig. 1: The fluid domain with boundary conditions of the bell nozzle.
Fig. 1: The fluid domain with boundary conditions of the bell nozzle.

Key Findings

수치 해석 결과, 높은 노즐 압력비(NPR)에서 점성 유동에 의한 비대칭 람다 충격파가 포착되었습니다. 설계 조건에서 마하수 3.0을 목표로 하였으나, 실제 수치 해석에서는 냉각 유동 시 2.932, 가열 유동 시 2.897의 마하수가 도출되어 이론값과 약 2.2~3.4%의 오차를 보였습니다. 비설계 조건에서는 배압이 15 bar에서 5 bar로 감소함에 따라 충격파가 노즐 출구 방향으로 이동하며 강도가 세지는 현상이 관찰되었으며, 1D 이론 대비 약 20%의 위치 편차가 확인되었습니다.

Industrial Applications

본 연구의 결과는 로켓 엔진 노즐 설계 시 비설계 작동 조건에서 발생하는 성능 손실을 예측하고 최소화하는 데 활용될 수 있습니다. 특히 충격파와 경계층의 상호작용으로 발생하는 측면 하중(side loads) 및 유동 박리 지점을 정확히 파악함으로써 노즐의 구조적 안정성과 추진 효율을 높이는 기초 설계 데이터로 사용 가능합니다.

Theoretical Background

특성곡선법 (Method of Characteristics, MOC)

MOC는 초음속 유동 문제를 해결하기 위한 수치적 절차로, 균일하고 파동이 없는 초음속 유동을 생성하는 노즐 윤곽을 설계하는 데 필수적입니다. 이 방법은 유동장의 특성 라인을 따라 미분 방정식을 적분하여 노즐의 확산부 좌표를 계산하며, 등엔트로피, 비점성, 비회전 유동을 가정하여 최적의 팽창 효율을 제공하는 곡선을 정의합니다.

유동 박리 및 충격파 구조 (Flow Separation and Shock Structures)

노즐 내부 배압이 설계치보다 높을 경우 과팽창(Over-expanded) 상태가 되어 유동 박리가 발생합니다. 이때 자유 충격 박리(FSS)와 제한 충격 박리(RSS) 현상이 나타나며, 이는 람다 충격파(Lambda shock)와 마하 디스크(Mach disk)를 형성합니다. 이러한 복잡한 충격파 구조는 노즐 벽면의 압력 분포를 급격히 변화시키며 엔진의 추력 성능에 직접적인 영향을 미칩니다.

Results and Analysis

Experimental Setup

수치 해석을 위해 노즐 벽면 소재는 티타늄으로 설정하였으며, 25,800개의 구조적 격자(Structured grid)를 생성하여 해석의 정밀도를 높였습니다. 난류 모델로는 유동 박리 예측에 탁월한 성능을 보이는 SST k-ω(Shear Stress Transport) 모델을 적용하였습니다. 경계 조건으로 입구는 압력 입구(Pressure Inlet), 출구는 압력 출구(Pressure Outlet)를 사용하였으며, 2차 상향풍 스킴(Second order upwind scheme)을 통해 수렴성을 확보하였습니다.

Visual Data Summary

마하수 등고선 분석 결과, 노즐 목(Throat) 영역에서 반원형의 유동 특성이 관찰되었으며 확산부 중반 이후 유동 변화가 완만해지는 벨 노즐 특유의 설계 특징이 확인되었습니다. 비설계 조건인 15 bar 배압에서는 노즐 내부 깊숙이 충격파가 형성되었으나, 배압이 낮아질수록 충격파가 출구 쪽으로 밀려나며 유동장이 확장되는 시각적 변화가 뚜렷하게 나타났습니다.

Fig. 2: Variations in Mach number along the axis of Bell nozzle for cold flow.
Fig. 2: Variations in Mach number along the axis of Bell nozzle for cold flow.

Variable Correlation Analysis

배압(Back Pressure)과 충격파 위치 사이에는 강한 상관관계가 존재합니다. 배압이 낮아질수록 노즐 내부의 정압과 외부 압력의 차이가 줄어들어 충격파 발생 지점이 하류로 이동합니다. 또한, 가열 유동(수소)의 경우 냉각 유동(공기)보다 전압 손실이 약 1% 정도 더 크게 발생(5.46% vs 4.49%)하는 것으로 나타났는데, 이는 고온 유동에서의 경계층 효과와 점성 소산이 더 활발하기 때문으로 분석됩니다.

Paper Details

Computational Analysis of Bell Nozzles

1. Overview

  • Title: Computational Analysis of Bell Nozzles
  • Author: Beena D. Baloni, Sonu P. Kumar, S. A. Channiwala
  • Year: 2017
  • Journal: Proceedings of the 4th International Conference of Fluid Flow, Heat and Mass Transfer (FFHMT’17)

2. Abstract

벨형 노즐은 로켓 노즐에서 가장 흔히 사용되는 형태입니다. 이 유형의 노즐은 원추형 노즐에 비해 크기와 성능 면에서 상당한 이점을 제공할 뿐만 아니라 환형 노즐에 비해 복잡성을 줄여줍니다. 노즐은 연소실에서 생성된 정체 온도($T_o$)와 정체 압력($P_o$)을 사용하여 연소 가스를 높은 초음속 속도로 가속함으로써 추력을 생성합니다. 노즐 팽창비는 출구 속도에 의해 결정됩니다. 비행 중 제트 유동은 짧은 기간 동안만 주변 유동에 이상적으로 팽창 및 적응하며, 나머지 시간 동안 로켓 엔진은 비설계 조건에서 작동합니다. 본 연구는 GAMBIT 2.4.6 및 FLUENT 6.3.26 소프트웨어를 사용하여 설계 및 비설계 조건에서 벨형 노즐 내의 2D 축대칭 유동 해석을 포함합니다. 해석을 위한 고효율 노즐 윤곽을 정의하기 위해 특성곡선법과 유선 함수를 사용하는 컴퓨터 코드가 개발되었습니다. 시뮬레이션은 냉각 및 가열 유동 조건에 대해 별도로 수행되었습니다. 유동 해석을 위해 SST k-ω 난류 모델이 선택되었습니다. 수렴된 솔루션은 점성 유동에 대해 높은 노즐 압력비(NPR)에서 노즐 내 비대칭 람다 충격파를 포착했습니다. 또한 NPR에 따른 후속 충격파 및 유동 박리를 예측했습니다. 점성 예측에서 마하 스템의 법선 충격파 강도는 일반적으로 NPR이 증가함에 따라 증가합니다. 충격 구조, 충격 위치, 법선 충격 크기, 후속 충격 및 비대칭 람다 충격 측면에서 예측된 시뮬레이션과 분석 결과 사이에 양호한 일치가 관찰되었습니다.

3. Methodology

3.1. 노즐 윤곽 설계: C 언어 프로그램을 통해 MOC와 유선 함수를 결합하여 마하수 3.0에 최적화된 벨 노즐 좌표를 생성함.
3.2. 수치 도메인 및 격자 생성: GAMBIT 2.4.6을 사용하여 25,800개의 사각형 구조적 격자를 생성하고 노즐 목 부근에 격자를 조밀하게 배치함.
3.3. 물리 모델 설정: FLUENT 6.3.26에서 2D 정상 상태 RANS 방정식을 풀었으며, 박리 유동 해석에 적합한 SST k-ω 난류 모델을 적용함.
3.4. 경계 조건 및 수렴: 압력 입구 및 출구 조건을 설정하고, 잔차 기준을 $10^{-4}$로 설정하여 수렴될 때까지 반복 계산을 수행함.

4. Key Results

설계 조건에서 냉각 유동의 출구 마하수는 2.932, 가열 유동은 2.897로 측정되어 이론적 설계치인 3.0과 매우 근접한 결과를 얻었습니다. 비설계 조건(배압 5~15 bar) 분석에서는 배압이 높을수록 노즐 내부에서 충격파가 일찍 발생하며, 이로 인해 유동 박리가 유도됨을 확인하였습니다. 특히 15 bar 조건에서 충격파 이후 마하수가 2.66에서 0.498로 급격히 감소하는 현상이 관찰되었습니다. 수치 해석 결과는 1D 비점성 이론과 비교했을 때 충격파 위치에서 약 20%의 편차를 보였는데, 이는 실제 유동의 점성 효과와 경계층 성장을 반영한 결과입니다.

Figure List

  1. Fig. 1: 벨 노즐의 경계 조건이 포함된 유체 도메인
  2. Fig. 2: 냉각 유동에 대한 벨 노즐 축 방향 마하수 변화
  3. Fig. 3: 냉각 유동의 다양한 비설계 배압에 따른 마하수 등고선
  4. Fig. 4: 가열 유동의 다양한 비설계 배압에 따른 마하수 등고선
  5. Fig. 5: 15 bar 배압에서 노즐 축 및 벽면의 정압 변화

References

  1. G. Hagemann et al., “Advanced Rocket Nozzles,” Journal of Propulsion and Power, 1998.
  2. G. P. Sutton, “Rocket Propulsion Elements,” 7th ed, 2001.
  3. J. Ostlund, “Flow Process in Rocket Engine Nozzles with Focus on Flow Separation and Side-Loads,” 2002.

Technical Q&A

Q: SST k-ω 난류 모델을 사용한 특별한 이유가 있습니까?

SST k-ω 모델은 표준 k-ε 모델의 원거리 유동 해석 장점과 k-ω 모델의 벽면 근처 해석 장점을 결합한 모델입니다. 특히 로켓 노즐 내부에서 발생하는 강력한 역압력 구배에 의한 유동 박리 현상을 예측하는 데 있어 다른 난류 모델보다 실험 데이터와 더 높은 일치성을 보이기 때문에 선택되었습니다.

Q: 냉각 유동과 가열 유동의 해석 결과에서 나타난 주요 차이점은 무엇입니까?

가열 유동(수소, 1000K)은 냉각 유동(공기, 300K)에 비해 정체 압력 손실이 약 1% 더 높게 나타났습니다(5.46% vs 4.49%). 이는 고온 유동에서 점성 효과가 더 크게 작용하여 경계층 두께에 영향을 미치기 때문이며, 이로 인해 출구 마하수 또한 가열 유동에서 약간 더 낮게 관찰되었습니다.

Q: 비설계 조건에서 발생하는 람다 충격파(Lambda shock)의 구조적 특징은 무엇입니까?

람다 충격파는 입사 충격파, 반사 충격파, 그리고 마하 스템(Mach stem)이 만나는 삼중점(Triple point)을 형성하는 구조를 가집니다. 본 연구에서는 배압 5 bar 조건에서 이러한 분기된 구조를 명확히 포착하였으며, 이는 노즐 벽면의 경계층 박리와 밀접하게 연관되어 비대칭적인 압력 분포를 유발합니다.

Q: 1D 이론적 결과와 CFD 수치 해석 결과 사이에 편차가 발생하는 원인은 무엇입니까?

1D 이론은 유동을 비점성(Inviscid) 및 등엔트로피 과정으로 가정하지만, 실제 CFD 해석에서는 벽면 마찰에 의한 경계층 성장과 점성 소산 효과를 고려합니다. 이러한 점성 효과가 유동의 유효 단면적을 줄이고 에너지 손실을 유발하여 충격파의 위치와 강도에서 약 20%의 차이를 만들게 됩니다.

Q: 연구에서 언급된 ‘측면 하중(Side loads)’은 노즐에 어떤 영향을 미칩니까?

비설계 조건에서 유동 박리가 비대칭적으로 발생하면 노즐 벽면에 불균일한 압력 분포가 형성됩니다. 이는 노즐 구조물에 수직 방향의 힘인 측면 하중을 발생시키며, 심할 경우 노즐의 기계적 파손이나 비행 궤적의 불안정성을 초래할 수 있습니다.

Conclusion

본 연구는 벨형 로켓 노즐의 설계 및 비설계 작동 조건에 대한 2차원 축대칭 수치 해석을 성공적으로 수행하였습니다. MOC를 통해 설계된 노즐은 목표 마하수에 근접한 성능을 보였으며, SST k-ω 모델을 통해 과팽창 조건에서의 복잡한 충격파 구조와 유동 박리 현상을 정밀하게 모사하였습니다.

결론적으로, 노즐 내부의 충격파 위치와 강도는 배압에 따라 민감하게 변하며, 이는 1D 이론만으로는 예측하기 어려운 점성 효과를 포함하고 있습니다. 본 연구 데이터는 향후 다양한 고도에서 작동해야 하는 로켓 엔진의 노즐 최적화 및 구조적 안정성 검토를 위한 중요한 수치적 근거를 제공합니다.

Source Information

Citation: Beena D. Baloni, Sonu P. Kumar, S. A. Channiwala (2017). Computational Analysis of Bell Nozzles. Proceedings of the 4th International Conference of Fluid Flow, Heat and Mass Transfer (FFHMT’17).

DOI/Link: 10.11159/ffhmt17.110

Technical Review Resources for Engineers:

▶ Access the original research paper (PDF)
▶ FLOW-3D 솔루션 팀과 협의하여 기술적 타당성을 검토하시려면..

This material is for informational purposes only. Unauthorized commercial use is prohibited.
Copyright © 2025 STI C&D. All rights reserved.