SIMULATION OF TURBULENT FLOW AND HEAT TRANSFER OVER A BACKWARD-FACING STEP WITH RIBS TURBULATORS

본 연구는 리브 터뷸레이터가 장착된 채널 내 후향 계단 유동(backward-facing step flow) 및 열전달 특성을 수치적으로 분석한 결과를 제시합니다. 후향 계단 유동은 유동 박리와 재부착이 발생하는 복잡한 물리적 현상을 포함하고 있어, 공학적으로 매우 중요한 연구 대상입니다. 특히 가스 터빈 블레이드의 내부 냉각로나 열교환기 설계에서 열전달 효율을 극대화하기 위해 리브와 같은 난류 촉진 장치가 널리 사용됩니다. 본 논문은 최대 32,000의 레이놀즈 수 범위에서 계단 높이, 리브의 개수 및 두께가 유동장과 온도장에 미치는 영향을 체계적으로 조사하였습니다. 지배 방정식으로는 연속 방정식, Navier-Stokes 방정식 및 에너지 방정식을 사용하였으며, 유한 체적법과 SIMPLE 알고리즘을 통해 이산화하였습니다. 난류 모델링을 위해 벽 함수가 포함된 표준 k-ε 모델을 적용하여 수치적 정확도를 확보하였습니다. 연구 결과, 수축비(contraction ratio)의 증가가 재순환 영역의 강도와 크기에 결정적인 영향을 미친다는 점을 확인하였습니다. 이러한 분석은 복잡한 기하학적 구조 내에서의 열전달 제어 및 최적화 설계에 중요한 기초 자료를 제공합니다.

메타데이터 및 키워드

Figure 1: 고려된 문제의 개략도 (H=0.05m, L=0.4m, x1=0.0492m, H/w=11, P=0.1). 기하학적 구조와 주요 치수를 정의합니다.
Figure 1: 고려된 문제의 개략도 (H=0.05m, L=0.4m, x1=0.0492m, H/w=11, P=0.1). 기하학적 구조와 주요 치수를 정의합니다.

논문 메타데이터

  • Industry: 기계 공학 (Mechanical Engineering)
  • Material: 공기 (Air)
  • Process: 난류 유동 및 열전달 시뮬레이션
  • System: 리브 터뷸레이터가 포함된 후향 계단 채널
  • Objective: 계단 높이, 리브 개수 및 두께가 유동 및 열전달 특성에 미치는 영향 조사

핵심 키워드

  • 후향 계단 (backward facing)
  • 리브 (ribs)
  • 난류 덕트 유동 (turbulent duct flow)
  • 열전달 향상 (heat transfer enhancement)
  • k-ε 모델 (k-ε model)

핵심 요약

연구 구조

본 연구는 리브 터뷸레이터가 하단 벽면에 수직으로 배열된 후향 계단 채널 내의 난류 유동을 모델링하는 수치적 구조를 가집니다.

방법 개요

SIMPLE 알고리즘과 표준 k-ε 난류 모델을 결합한 유한 체적법(FVM)을 사용하여 지배 방정식을 해결하였습니다.

주요 결과

수축비(SR)가 0.25에서 0.5로 증가함에 따라 계단 뒤쪽의 재순환 영역 크기가 증가하였으며, 리브를 추가함으로써 매끄러운 후향 계단 대비 열전달 성능이 현저히 향상되었습니다. 또한, 리브의 폭이 좁을수록(H/w 값이 클수록) 국소 Nusselt 수가 증가하는 경향을 보였습니다.

산업적 활용 가능성

가스 터빈 블레이드의 내부 냉각 채널, 고효율 열교환기, 전자 기기의 냉각 시스템 설계에 직접적으로 적용 가능합니다.

한계와 유의점

본 시뮬레이션은 공기의 물성치가 일정하다는 Boussinesq 근사를 가정하였으며, 주로 2차원 유동 분석에 집중되어 있습니다.


논문 상세 정보

1. 개요

  • Title: SIMULATION OF TURBULENT FLOW AND HEAT TRANSFER OVER A BACKWARD-FACING STEP WITH RIBS TURBULATORS
  • Author: Khudheyer S. MUSHATET
  • Year: 2011
  • Journal: THERMAL SCIENCE
  • DOI/Link: 논문에 명시되지 않음

2. 초록

리브 터뷸레이터가 있는 채널 내부의 후향 계단 유동 및 열전달에 대한 시뮬레이션을 제시합니다.

이 문제는 최대 32,000의 레이놀즈 수에 대해 조사되었습니다.

계단 높이, 리브의 수 및 리브 두께가 유동 및 온도장에 미치는 영향을 조사하였습니다.

지배 방정식인 연속, 전체 Navier-Stokes 및 에너지 방정식을 이산화하기 위해 엇갈림 격자 기술을 사용하는 제어 체적법이 도입되었습니다.

난류의 영향은 벽 함수 공식과 함께 k-ε 모델을 사용하여 모델링되었습니다.

얻어진 결과는 수축비의 증가(즉, 계단 높이의 증가)에 따라 계단 뒤의 재순환 영역의 강도와 크기가 증가함을 보여줍니다.

리브 이후의 재순환 영역 크기와 재부착 길이는 수축비가 증가함에 따라 감소합니다.

3. 방법론

지배 방정식: 공기를 작동 유체로 하며, Boussinesq 근사를 적용한 전체 Navier-Stokes, 에너지 및 연속 방정식을 사용하여 유동을 정의하였습니다.

난류 모델링: Launder 등이 제안한 표준 k-ε 모델을 사용하여 난류 운동 에너지(k)와 소산율(ε)에 대한 수송 방정식을 해결하였으며, 모델 계수는 σk=1.0, σε=1.3, Cμ=0.09 등을 적용하였습니다.

수치 해석 기법: 비균일 엇갈림 격자 시스템에서 유한 체적법을 수행하였으며, 압력-속도 결합을 위해 SIMPLE 알고리즘을 사용하고 수렴 기준을 10^-5로 설정하였습니다.

4. 결과 및 분석

수축비의 영향: 수축비(SR)가 0.25, 0.35, 0.5로 변화함에 따라 계단 뒤의 재순환 영역은 강화되지만, 리브 이후의 재부착 길이는 오히려 단축되는 역동적인 변화를 확인하였습니다.

열전달 성능 분석: 모든 레이놀즈 수 범위에서 리브의 존재는 열전달율을 크게 향상시켰으며, 특히 계단 인근 영역에서 수축비가 클수록 국소 Nusselt 수가 높게 나타났습니다.

리브 기하학적 효과: 리브의 폭(w)이 좁을수록(H/w=11 vs H/w=5) 유동 교란이 효과적으로 발생하여 국소 Nusselt 수가 증가하는 결과를 얻었습니다.

Figure 2 Computed velocity vectors for 2 ribs and different values of contraction ratios,
Figure 2 Computed velocity vectors for 2 ribs and different values of contraction ratios,
Figure 3 Computed velocity vectors for 3 ribs and different values of contraction ratios,
Figure 3 Computed velocity vectors for 3 ribs and different values of contraction ratios,

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

  • Figure 1: 고려된 문제의 개략도 (H=0.05m, L=0.4m, x1=0.0492m, H/w=11, P=0.1). 기하학적 구조와 주요 치수를 정의합니다.
  • Figure 6: 3개의 리브와 SR=0.5 조건에서 레이놀즈 수에 따른 축 방향 속도 분포. 레이놀즈 수가 증가할수록 최대 속도와 재순환 강도가 증가함을 보여줍니다.
  • Figure 11: SR=0.5, Re=16000 조건에서 다양한 사례의 국소 Nu 변화 비교. 리브가 매끄러운 계단 유동보다 열전달을 크게 향상시킴을 입증합니다.
  • Figure 14: 현재 시뮬레이션과 기존 실험 데이터(Lio 등, H/B=1, Re=6000)의 비교. 수치 모델의 타당성을 검증하는 중요한 자료입니다.

6. 참고문헌

  • Lio, T., Hwang, J. (1992). Developing Heat Transfer and Friction in a Ribbed Rectangular Duct with Flow Separation at Inlet. ASME J. Heat Transfer. 114(3), pp. 546-573.
  • Launder, B. E., Spalding, D. B. (1974). The Numerical Computation of Turbulent Flows. Computer Methods in Applied Mechanics and Engineering.
  • Versteege, H. K., Malalasekera, W. (1995). An Introduction of Computational Fluid Dynamics. Hemisphere Publishing Corporation.

기술 Q&A (Technical Q&A)

Q: 수축비(SR)가 계단 뒤의 재순환 영역에 미치는 영향은 무엇입니까?

연구 결과에 따르면 수축비, 즉 계단의 높이가 증가할수록 계단 바로 뒤에 형성되는 재순환 영역의 강도와 크기가 증가합니다. 이는 유동 단면적의 급격한 변화가 유동 박리를 더욱 강력하게 유도하기 때문입니다. 반면, 리브 뒤쪽에서 발생하는 재순환 영역의 크기와 재부착 길이는 수축비가 증가함에 따라 오히려 감소하는 경향을 보입니다. 이러한 특성은 계단 높이 조절을 통해 특정 영역의 유동 구조를 제어할 수 있음을 시사합니다.

Q: 리브의 폭(width)이 열전달 효율에 어떤 영향을 미칩니까?

본 논문의 Figure 13 분석에 따르면, 리브의 폭이 좁을수록(H/w 값이 클수록) 국소 Nusselt 수가 더 높게 나타납니다. 구체적으로 H/w=11인 경우가 H/w=5인 경우보다 더 우수한 열전달 성능을 보였습니다. 이는 얇은 리브가 유동의 박리와 재부착을 더 빈번하게 유도하여 경계층을 효과적으로 파괴하기 때문으로 해석됩니다. 따라서 열전달 극대화를 위해서는 적절한 두께의 리브 설계가 필수적입니다.

Q: 수치 해석에서 압력-속도 결합을 위해 어떤 알고리즘이 사용되었습니까?

본 연구에서는 압력과 속도 항의 결합 문제를 해결하기 위해 CFD 분야에서 널리 사용되는 SIMPLE(Semi-Implicit Method for Pressure-Linked Equations) 알고리즘을 적용하였습니다. 이 알고리즘은 엇갈림 격자(staggered grid) 시스템 위에서 계산되었으며, 속도, 압력, 온도 및 난류 항에 대해 각각 0.5에서 0.8 사이의 이완 계수(relaxation factors)를 사용하여 수렴 안정성을 확보하였습니다. 수렴 판정 기준은 모든 변수에 대해 10^-5 이하로 설정되었습니다.

Q: 난류 모델링을 위해 선택된 모델과 그 이유는 무엇입니까?

연구자는 표준 k-ε(k-epsilon) 모델을 선택하여 난류 유동을 시뮬레이션하였습니다. 이 모델은 난류 운동 에너지(k)와 그 소산율(ε)을 계산하는 두 개의 수송 방정식을 기반으로 하며, 공학적 유동 해석에서 계산 효율성과 정확도 사이의 균형이 잘 잡힌 모델로 평가받습니다. 특히 벽면 근처의 복잡한 유동을 처리하기 위해 벽 함수(wall function) 공식을 함께 사용하여 벽면에서의 점성 효과를 적절히 반영하였습니다.

Q: 리브가 없는 일반적인 후향 계단 유동과 비교했을 때 리브의 효과는 어떠합니까?

Figure 11의 비교 데이터를 통해 알 수 있듯이, 리브 터뷸레이터를 추가한 경우가 리브가 없는 매끄러운 후향 계단 유동에 비해 Nusselt 수가 훨씬 높게 나타납니다. 리브는 유동에 지속적인 교란을 발생시켜 난류 강도를 높이고, 이는 벽면 근처의 열 확산을 촉진하는 역할을 합니다. 결과적으로 리브의 배치는 시스템의 압력 손실을 다소 증가시킬 수 있으나, 열전달 성능 면에서는 압도적인 이점을 제공합니다.

Q: 본 연구의 결과가 실제 산업 현장에서 어떻게 활용될 수 있습니까?

이 연구 결과는 고온에서 작동하는 가스 터빈 블레이드의 내부 냉각 통로 설계에 직접 활용될 수 있습니다. 블레이드 내부의 구불구불한 채널(serpentine channels)에 리브를 최적으로 배치함으로써 냉각 효율을 높여 부품의 수명을 연장할 수 있습니다. 또한, 콤팩트한 열교환기나 고출력 전자 기기의 방열판 설계 시, 유동 박리 영역을 제어하여 냉각 성능을 최적화하는 가이드라인으로 사용될 수 있습니다.

결론

본 연구는 리브 터뷸레이터가 장착된 후향 계단 채널에서의 난류 유동 및 열전달 특성을 수치적으로 규명하였습니다. 주요 결론으로, 리브의 추가는 열전달 성능을 획기적으로 향상시키며, 이러한 향상 효과는 레이놀즈 수와 수축비가 증가할수록 더욱 뚜렷해진다는 점을 확인하였습니다. 특히 수축비는 계단 뒤의 재순환 영역 크기를 결정하는 핵심 변수이며, 리브의 기하학적 형상(폭) 또한 국소 열전달 분포에 유의미한 영향을 미친다는 사실을 입증하였습니다.

공학적 관점에서 본 연구는 복잡한 분리 유동이 발생하는 시스템에서 리브를 활용한 열전달 제어의 가능성을 제시하였습니다. 다만, 본 연구는 2차원 시뮬레이션과 일정한 물성치 가정을 바탕으로 하고 있으므로, 향후 실제 작동 조건에서의 3차원 효과 및 가변 물성치를 고려한 추가 연구가 필요할 것으로 보입니다. 그럼에도 불구하고, 본 결과는 가스 터빈 및 열교환기 설계 최적화를 위한 중요한 수치적 근거를 제공합니다.


출처 정보 (Source Information)

Citation: Khudheyer S. MUSHATET (2011). SIMULATION OF TURBULENT FLOW AND HEAT TRANSFER OVER A BACKWARD-FACING STEP WITH RIBS TURBULATORS. THERMAL SCIENCE.

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.