다중 얇은 핀을 이용한 차등 가열 공동 내 자연 대류 흐름 불안정 및 열전달 수치 시뮬레이션

초록

본 연구는 차등 가열된 직사각형 공동에 얇은 핀을 부착하여 흐름 불안정과 열전달 특성을 수치적으로 분석한다. 유한체적법과 SIMPLE 알고리즘으로 온도·압력·속도장을 계산하고, 에너지 구배 이론의 K함수를 이용해 불안정 발생 메커니즘을 규명한다. 핀 길이·위치·수와 레이놀즈 수(Ra)의 변화가 열전달 효율에 미치는 영향을 조사했으며, Kmax와 Ra 사이의 선형 관계를 도출하였다.

상세 분석

본 논문은 자연 대류가 지배하는 차등 가열 공동 내부에 얇은 핀을 고정함으로써 발생하는 흐름 불안정 메커니즘을 에너지 구배 이론(Energy Gradient Theory, EGT) 기반의 K함수로 정량화하였다. 먼저, 2차원 직사각형 공동을 설정하고, 좌·우 벽을 각각 고온·저온으로 유지한다. 초기 온도는 전체 유체와 동일하게 설정하여 온도 구배가 시간에 따라 자연스럽게 형성되도록 하였다. 유한체적법(FVM)과 SIMPLE(Pressure‑Velocity Coupling) 알고리즘을 이용해 연속 방정식, 운동량 방정식, 에너지 방정식을 수치적으로 풀었으며, 격자 독립성 검증을 통해 2차 정확도를 확보하였다.

핀이 없는 기본 경우와 비교하여, 핀의 길이(Lf), 위치(Hf), 개수(Nf) 및 레이놀즈 수(Ra) (즉, 온도 차에 의한 부력 구동) 파라미터를 체계적으로 변형하였다. 결과는 다음과 같다. 첫째, K함수는 온도 등고선에서 급격히 변하는 영역, 즉 열경계층이 얇아지고 플라즈마가 형성되는 곳에서 최대값을 보였다. 이는 EGT가 제시하는 “에너지 구배가 클수록 흐름이 불안정해진다”는 가설을 실험적으로 확인한 것이다. 둘째, 핀 길이가 Ra가 낮은 경우(예: Ra < 10⁶)에는 열전달 저항을 증가시켜 Nusselt 수(Nu)를 감소시켰지만, Ra가 충분히 높아(예: Ra > 10⁸) 흐름이 강하게 전도되면서 핀 주변에 발생하는 와류가 열전달을 촉진해 Nu가 증가하였다. 셋째, 핀의 최적 위치는 공동 높이의 중앙(Hf ≈ 0.5·H)으로, 이때 Kmax가 가장 크게 나타나며, 플라즈마 형성 및 재순환 구조가 가장 활발히 발생한다. 넷째, 핀의 개수를 늘릴수록 초기에는 열전달 저항이 누적돼 Nu가 감소하지만, Ra가 충분히 큰 경우 다중 핀 사이에 복합적인 와류가 형성되어 전체 열전달이 향상된다. 마지막으로, Kmax와 Ra 사이에 선형 관계(Kmax ≈ α·Ra + β)를 도출함으로써, 전통적인 차원 분석(예: Nu‑Ra‑Pr 관계)과는 별개의 새로운 물리적 지표를 제공한다. 이러한 결과는 자연 대류 설계에서 핀 구조를 최적화함으로써 열전달 효율을 극대화할 수 있는 실용적인 가이드를 제시한다.