[논문 리뷰] The pear-shaped fate of an ice melting front

James N. Hewett, Mathieu Sellier|arXiv (Cornell University)|2017. 05. 07.

Icing and De-icing Technologies참고 문헌 9인용 수 3

한 줄 요약

이 연구는 수증기의 밀도 역전점(4.0 °C) 근처에서 가열된 수평 실린더 주위의 얼음 융해를 수치적으로 시뮬레이션하며, 동적 메시 기반의 노드 재배치 알고리즘을 사용해 변하는 고체-액체 계면을 추적한다. 결과적으로 열 플룸으로 인해 배추 모양의 융해 front 가 형성되며, T_w > 4.0 °C일 경우 상향, T_w < 4.0 °C일 경우 하향으로 발생한다. 이는 자연대류가 융해 속도와 계면 형태에 상당한 영향을 미친다는 것을 확인하며, 형상과 전체 융해 진행 상황에서 실험 데이터와 양호한 일치를 보인다.

ABSTRACT

A fluid-structure interaction problem with the melting of water around a heated horizontal circular cylinder is analysed with numerical simulations. Dynamic meshing was used for evolving the flow domain in time as the melting front extended radially outward from the cylinder; a node shuffle algorithm was used to retain mesh quality across the significant mesh deformation. We simulated one case above the density inversion point of water and one case below, yielding pear-shaped melting fronts due to thermal plumes either rising or falling from the cylinder, respectively. Results were compared with previous experimental studies and the melting front profiles matched reasonably well and melting rates were in agreement. We confirm that natural convection plays a significant role in the transport of energy as the melt zone increases, and needs to be considered for accurately modelling phase change under these conditions.

연구 동기 및 목표

• 가열된 수평 실린더 주위의 상전이 과정에서 수증기의 밀도 역전이 융해 front 의 형태와 역학에 미치는 영향을 조사하기 위해.
• 자연대류 효과를 고려한 명시적 계면 추적 기반의 스타인프 문제 수치 시뮬레이션을 검증하기 위해.
• 상위 및 상위 이하 역전점 온도 조건에서 실험 데이터와의 비교를 통해 시뮬레이션된 융해 front 형상과 융해 속도를 분석하기 위해.
• 융해 영역이 확장됨에 따라 자연대류가 열전달을 어떻게 향상시키는지 평가하기 위해.
• 큰 변형을 겪는 이동 경계 문제에 대해 노드 재배치 기반 동적 메시 기법의 정확도를 평가하기 위해.

제안 방법

• 비압축성 및 난류 흐름 모델을 2차원으로 사용하였으며, 부력에 의한 대류를 고려하기 위해 부소니크스 근사법을 적용하였다.
• 이동 경계(융해 front)는 동적 메시 기법을 통해 명시적으로 추적하였으며, 메시 품질을 유지하기 위해 노드 재배치 알고리즘을 사용하였다.
• 각 시간 단계에서, 메시가 업데이트된 상태에서 유동장을 정적인 조건으로 해석하여 정적 해법을 적용하였으며, 준정적인 조건을 가정하였다.
• 물의 열물성 특성은 밀도를 제외하고는 일정하다고 가정하였으며, 밀도는 온도에 따라 선형적으로 변화시켜 밀도 역전 효과를 반영하였다.
• 계면의 운동은 국소 열역학적 유량과 스타인프 조건에 의해 결정되었으며, 인위적인 상전이 소스 항은 사용하지 않았다.
• 실험 조건과 일치시키기 위해 두 가지 벽 온도 조건(T_w = 2.3 °C, T_w = 14.1 °C)에서 시뮬레이션을 수행하였다.

실험 결과

연구 질문

• RQ1 수증기의 밀도 역전점 근처(약 4.0 °C)에서 자연대류 효과가 가열된 수평 실린더 주위의 융해 front 형상과 진화에 미치는 영향는 어떠한가?
• RQ2 융해 영역이 확장됨에 따라 자연대류가 열전달에 얼마나 지배적인 역할을 하는가? 이는 밀도 역전점 이하 및 이상에서 어떻게 다를까?
• RQ3 노드 재배치 기반의 동적 메시 기법이 실험적으로 관측된 복잡한 배추 모양의 계면 형태를 정확히 포착할 수 있는가?
• RQ4 왜 시뮬레이션은 도전성 주도에서 대류 주도로의 전이 단계에서 실험 결과와는 다르게 융해 속도에 차이를 보이는가?
• RQ5 2차원 정적인 가정이 3차원 실험에서 관측된 일시적 난류 불안정성과 열전달 향상 현상을 예측하는 데 어떤 영향을 미치는가?

주요 결과

• T_w = 14.1 °C일 경우 상향, T_w = 2.3 °C일 경우 하향으로 열 플룸으로 인해 두 경우 모두 배추 모양의 융해 front 가 나타났으며, 실험 관측과 일치하였다.
• 레이놀즈 수는 각각 6700(역전)과 7400(정상)이었으며, 서로 다른 유동 구조에도 불구하고 유사한 대류 강도를 나타내었다.
• 융해 부피가 증가함에 따라 자연대류가 주요 열전달 메커니즘이 되었으며, 이는 계면 진화에 상당한 영향을 미쳤다.
• 전체 시뮬레이션 기간 동안 실험적 프로파일과 합리적인 일치를 보였으며, 모델의 타당성을 확인하였다.
• τ ≈ 0.3에서 도전성 주도에서 대류 주도로의 전이 단계에서 융해 속도의 불일치는 2차원 정적 모델에 존재하지 않는 비정상적 3차원 효과와 일시적 불안정성 때문으로 기인되었다.
• 정상 배추 모양에서는 밀도 역전점의 포함로 인해 실린더 상단에 반대 방향으로 회전하는 소용돌이가 더 복잡하게 나타났다.