Journal Archive

경사진 사각채널 내에서의 혼합대류는 태양열 집열기, 밀집형 열교환기 및 화학기상증착, 자동차 기관 및 금형장비의 냉각 등과 같은 많은 공학적 응용 분야에서 중요한 역할을 하기 때문에 많은 연구자들의 관심을 끌어 왔다. 특히 강제 및 자연대류의 상호 작용에 의한 혼합대류는 복잡한 유동 형태를 나타내어, 그 특성을 파악하는 것은 열 및 물질전달 장치의 설계에 있어 매우 중요하다. 따라서 많은 연구자들에 의해 Reynolds수, Rayleigh수, 경사각, 종횡비 및 Prandtl수 등의 영향에 대하여 다양하게 고찰되어져 왔다.

Kamotani 등¹⁾은 상부가 냉각되고 하부가 가열되는 수평 평행평판채널 내에서의 혼합대류 열전달 특성을 Pr=0.7인 공기를 작동유체로 하여 실험적으로 연구하였으며, Reynolds수(30≤Re≤1100), Rayleigh수(10³<Ra≤3.1×10⁴)rk 혼합대류에 미치는 영향을 검토하였고, Chiu 등^2,3)은 하부가 가열되는 수평 평행 채널 내에서의 혼합대류 열전달 특성에 대하여 작동유체를 공기로 하고 2.71×10²≤Gr/Re²≤6.79×10⁴ 범위 내에서 실험과 수치적 계산을 동시에 진행하여 그 결과를 고찰하였다. Nace 등⁴⁾은 하부가 가열되고 상부가 냉각되는 사각채널 내에서 AR=2, Ra=22,200, Re=18.75, 36, 54일 때 질소의 혼합대류 유동특성을 실험 및 수치해석적 방법으로 연구하였으며, 그 결과가 서로 잘 일치함을 보였다. Yu 등⁵⁾은 하부에서 가열되는 수평평행 채널 내에서의 혼합대류유동을 AR=12, Pr=0.7인 경우 Reynolds수 (20≤Re≤50)와 Rayleigh수 (Ra≤31000)를 변수로 하여 롤의 생성, 개수, 크기, 위치 등에 대하여 3차원 비정상 수치해석적으로 연구하였다. Chang 등⁶⁾은 수평 사각 덕트 내에서의 종방향 볼텍스 유동에 대하여 저 점성유체인 공기를 이용하여 종횡비의 영향에 대하여 체계적으로 실험하였다. Lin⁷⁾은 수평 사각 채널 내 부력에 의한 볼텍스 유동과 열전달 특성에 대하여 저 Reynolds수가 혼합대류 볼텍스 유동에 미치는 영향을 저 점성유체인 공기(Pr=0.7)를 작동유체로 하여 실험적으로 연구하였다. Ozsunar 등⁸⁾은 하부가 가열되는 수평 평행평판 내에서 Ar=6.67, Pr=0.7인 경우 Reynolds수(50<Re<1,000), Grashof수 (7×10⁵<Gr<4× 10⁶) 그리고 경사각(0≤θ≤90^o)이 혼합대류 유동특성에 미치는 영향을 3차원 수치해석인 방법으로 연구하였으며, 그 결과를 이전의 연구결과와 비교하여 서로 잘 일치함을 보여 주었다. 최근 Piao 등^9,10)은 고 점성유체를 사용하여 Reynolds수와 Rayleigh수가 수평채널 내의 혼합대류 볼텍스 유동에 미치는 영향을 3차원 수치해석적인 방법으로 연구하였다. 그러나 Piao 등^9,10)의 연구를 제외한 지금까지의 연구들은 대부분 저 점성유체인 공기, 질소나 물의 혼합대류 연구에 국한되었으며, 고 점성유체로 채워진 경사진 평행평판에서 경사각이 고 점성유체의 혼합대류에 미치는 영향에 대한 연구는 미흡하다. 따라서 기관 및 변압기의 냉각유체로 유용하게 사용되고 있는 고 점성유체에 관한 연구가 필요하다.

본 연구에서는 상부가 냉각되고 하부가 가열되는 고 점성유체로 채워진 경사진 사각채널 내에서 경사각이 저 Reynolds수에서의 혼합대류 볼텍스 유동 및 열전달에 미치는 영향을 유한체적법을 기초로 한 3차원 수치해석적인 방법으로 고찰하고자 한다.

본 연구는 상부가 냉각되고 하부가 가열되는 두 개의 등온 벽으로 구성된 사각채널 내에서 Pr=909, Ra=10^-3, AR=4.0일 때 경사각이 사각채널 내의 혼합대류 볼텍스 유동 및 열전달에 미치는 영향을 3차원 수치해석적인 방법으로 연구하였다.

Fig. 3은 Re=10^-3일 때 Y=0.5인 위치의 X-Z 평면에서 경사각 변화(0°≤θ≤45°)에 따른 채널 내의 유동장과 온도장을 나타낸다. 그림에서 보는 바와 같이 경사각 θ=0°인 경우 입구 유속이 상대적으로 매우 작기 때문에 전 유동장은 자연대류가 지배적인 유동을 하여 하류방향으로 천천히 이동하는 횡방향 롤 유동이다. 경사각이 증가하여 θ=15°인 경우, 경사각이 증가함으로써 수직방향 유속의 증가로 채널 중간부분에서의 횡방향 롤은 종방향 롤로 변화되어진다. 그리고 채널 입구에서는 부력에 의한 횡방향 롤 유동이 입구 유속에 의해 서로 나누어지는 것을 볼 수 있다. 경사각이 더욱 증가하면 채널 중간부분에서의 종방향 롤의 길이는 경사각의 증가로 인해 점차적으로 작아진다. 이는 지구 중력의 영향으로 경사각이 증가할수록 일정 위치를 정점으로 수직방향 유체유동 성분이 감소하기 때문인 것으로 생각된다. 그리고 채널 입구 및 출구에서는 횡방향 롤은 더욱 분명히 나누어져 구분되어지고, 그 크기가 증가하는 것을 볼 수 있다. 이로부터 AR=4.0, Pr=909, Re=1×10^-2, Ra=10⁴ 경사각 15°≤θ≤45° 범위에서 채널의 입구와 출구영역에서는 서로 구분된 횡방향 롤을 형성하고, 채널 중간부분에서는 종방향 롤을 형성하여 전체적인 유동 형태는 두 형태의 유동이 혼합된 혼합 롤 유동을 형성함을 알 수 있다.

Fig. 3 
Velocity magnitude and temperature distributions for various angles of inclination at AR=4.0, Pr=909, Re=1×10^-3, Ra=10⁴ and Y=0.5

Fig. 4는 경사각 θ=0°, 15°, 30°, 45°일 때 사각 채널 내의 혼합대류 볼텍스 유동장과 온도장의 속도벡터와 온도분포를 3차원적으로 나타낸 것이다. 그림에서 보는 바와 같이 θ=0°인 경우 3차원 사각채널의 X-Y평면상에 자연대류에 의한 횡방향 롤 유동이 형성됨을 잘 나타내고 있고, 경사각이 증가하면 종방향 롤 유동이 형성됨을 X-Z 및 X-Y면을 통해 잘 나타내고 있음을 확인 할 수 있다. 또한 경사각이 증가하면 가열된 하부평판으로부터 데워진 유체가 하부평판의 경사면을 따라 상부에 모여 하나의 열섬을 형성함을 알 수 있다. 이 데워진 유체는 저속유입과 저속유출 그리고 경사각의 영향으로 냉각된 경사진 상부평판을 따라 하부로 순환하면서 종방향 롤 유동에 의해 열전달이 이루어지고 있음을 알 수 있다.

Fig. 4 
Three-dimensional velocity vector and temperature distributions for various angles of inclination at AR=4.0, Pr=909, Re=1×10^-3, Ra=10⁴

Fig. 5는 Pr=909, Ra=10^-3, AR=4.0, Re=1×10^-3 일 때 경사각의 증가에 따른 채널 중심에서의 주 유동방향 속도분포를 그래프로 나타낸 것이다. 그림에서 보는 바와 같이 경사각 θ=0°인 경우 유동은 횡방향 롤 유동으로 속도분포는 주기적으로 변화하는 사인곡선을 나타낸다. 경사각이 증가함에 따라 주 유동방향 U속도는 증가하여 속도분포는 전부 양의 값을 나타내고 최대 속도크기는 채널 중심 영역에서 나타나며 경사각이 증가함에 따라 증가하였다.

Fig. 5 
U-velocity distributions at mid-heigh for various angles of inclination at AR=4.0, Pr=909, Ra=10⁴ and Re=1×10^-3

Fig. 6은 Pr=909, Ra=10^-3, AR=4.0, Re=1×10^-3 일 때 경사각에 따른 가열 벽에서의 국소 Nusselt수 분포를 그래프로 나타낸 것이다. 그림에서 보는 바와 같이 경사각 θ=0°에서는 횡방향 롤 유동으로 Nusselt수는 주기적인 사인곡선을 나타낸다. 경사각이 증가함에 따라 주 유동방향 속도 U가 증가함으로 입구영역에서는 Nusselt수는 증가하고 하류 방향으로 가면서 감소하지만 그 변화는 아주 작았다.

Fig. 6 
Local Nusselt number distributions of hot wall for various angles of inclination at AR=4.0, Pr=909, Ra=10⁴ and Re=1×10^-3

Fig. 7은 Pr=909, Ra=10^-3, AR=4.0, Re=4×10^-3 일 때 경사각에 따른 가열 벽에서의 국소 Nusselt수 분포를 나타낸 것이다. 이 경우는 Fig. 6에 비해 입구 유속이 증가하여 θ=0°인 경우 나타나는 자연대류에 의한 횡방향 롤은 사라지고 전 유동장은 안정적이고 대칭적인 분포를 나타내는 종방향 롤에 의해 지배되어지므로 더 이상 사인곡선은 나타나지 않는다. 또한 경사각이 증가함에 따라 Fig. 6과 같이 입구영역에서의 Nusselt수는 증가하고 하류 방향으로 가면서 감소하지만 그 변화는 아주 작았다.

Fig. 7 
Local Nusselt number distributions of hot wall for various angles of inclination at AR=4.0, Pr=909, Ra=10⁴ and Re=4×10^-3

상부가 냉각되고 하부가 가열되는 두 개의 등온 벽으로 구성된 사각채널 내에서 경사각이 사각 채널 내의 혼합대류 볼텍스 유동 및 열전달에 미치는 영향을 3차원 수치해석적인 방법으로 연구한 결과는 다음과 같다.

1) 유체의 유입속도가 낮고(Re=10^-3), 고 점성유체(Pr=909)로 채워진 사각채널 내에서의 볼텍스 유동은 경사각이 θ=0°인 경우 자연대류의 영향으로 횡방향 롤 유동을 형성하였다.

2) 경사각 θ를 15°에서 45°까지 변화시키면서 그 영향을 살펴본 결과는 경사각 15°≤θ≤45°범위에서 채널 입구와 출구에서는 횡방향 롤 유동이 형성되고, 중간부분에서는 종방향 롤 유동이 형성되는 혼합 롤 유동형태를 나타냈으며, 경사각이 증가함에 따라 사각채널 내에서 형성되는 종방향 롤의 길이는 작아졌다.

3) 경사각 θ=0°인 경우의 속도분포는 주기적으로 변화하는 사인곡선을 나타내었고, 경사각이 증가함에 따라 X방향 U속도는 증가하여 속도분포는 전부 양의 값을 가지며, 최대 속도는 채널 중심에서 형성되었고, 그 크기는 경사각이 증가함에 따라 증가하였다.

4) 경사각 θ=0°에서 Nusselt 수는 주기적인 사인곡선 형태를 나타내었고, 경사각 0°≤θ≤45°영역에서 Nusselt 수는 입구영역에서 최대값을 가지고 출구쪽으로 가면서 감소하였으며, 그 크기는 경사각이 증가함에 따라 증가하였다.