[논문 리뷰] Integrated modeling and validation for phase change with natural convection

K. Schüller, Benjamin Berkels외 1명|arXiv (Cornell University)|2018. 01. 11.

Arctic and Antarctic ice dynamics참고 문헌 17인용 수 3

한 줄 요약

이 논문은 OpenFOAM에서 固定 격자 잠열 다공도 방법과 고해상도 광학 영상 및 이진 Mumford-Shah 분할을 사용하여 물-ice 시스템에서 자연대류를 동반한 상변화를 검증하기 위한 통합 시뮬레이션-실험 프레임워크를 제시한다. 이 방법은 시뮬레이션된 상계면과 측정된 상계면 간의 연속적이고 정량적인 비교를 가능하게 하여 정성적으로는 양호한 일치를 보이고 있으나, 단순화된 경계 조건과 실험적 열손실로 인해 작은 정량적 차이가 존재하는 것으로 나타났다.

ABSTRACT

Water-ice systems undergoing melting develop complex spatio-temporal interface dynamics and a non-trivial temperature field. In this contribution, we present computational aspects of a recently conducted validation study that aims at investigating the role of natural convection for cryo-interface dynamics of water-ice. We will present a fixed grid model known as the enthalpy porosity method. It is based on introducing a phase field and employs mixture theory. The resulting PDEs are solved using a finite volume discretization. The second part is devoted to experiments that have been conducted for model validation. The evolving water-ice interface is tracked based on optical images that shows both the water and the ice phase. To segment the phases, we use a binary Mumford Shah method, which yields a piece-wise constant approximation of the imaging data. Its jump set is the reconstruction of the measured phase interface. Our combined simulation and segmentation effort finally enables us to compare the modeled and measured phase interfaces continuously. We conclude with a discussion of our findings.

연구 동기 및 목표

• 물-ice 시스템에서 자연대류를 동반한 상변화를 시뮬레이션하기 위한 계산 효율성이 높은 고정 격자 모델을 개발하기 위해.
• 수축성과 시간 해상도가 높은 실험을 수행하여 변화하는 물-ice 계면의 역학을 캡처하기 위해.
• 고도의 영상 분할 기법을 통해 시뮬레이션 및 실험 데이터를 통합하여 직접적인 모델 검증을 수행하기 위해.
• 시뮬레이션과 실험 간의 격차를 규명하고 분석하여 모델의 정밀도를 향상시키기 위해.
• 복잡한 다물리적 시스템을 위한 융통성 있고 통합적인 검증 전략을 수립하기 위해.

제안 방법

• 상변화를 모델링하기 위해 혼합 이론과 단계장 공식화에 기반한 고정 격자 기반의 잠열 다공도 방법이 사용된다.
• Boussinesq 근사와 함께 운동량, 에너지, 연속 방정식을 포함하는 PDE 시스템은 OpenFOAM에서 유한체적 이산화 방법으로 해석된다.
• 고체와 액체 상의 전이를 모델링하기 위해 마찰대역 상수를 도입하였으며, 민감도 분석 결과 물-ice 시스템에서는 높은 값이 필요함을 확인하였다.
• 융해 과정의 광학적 영상은 고공간-시간 해상도로 촬영되어 실험적으로 상계면을 추적한다.
• 영상 데이터에 대해 이진 Mumford-Shah 분할을 적용하여 조각별 일정한 근사의 점프 집합으로서의 상계면을 재구성한다.
• 분할된 상계면을 여러 시간대에서 시뮬레이션된 상계면과 직접 비교함으로써 정량적 검증을 가능하게 한다.

실험 결과

연구 질문

• RQ1 고정 격자 기반으로 잠열 다공도 방법이 물-ice 상변화에서 자연대류 효과를 얼마나 정확하게 시뮬레이션할 수 있는가?
• RQ2 시뮬레이션된 상계면 형태가 실험적으로 관측된 계면 진화와 어느 정도 일치하는가?
• RQ3 상계면 위치에서 시뮬레이션과 실험 간의 주요 격차 원인은 무엇인가?
• RQ4 마찰대역 상수가 물-ice 시스템에서 물리적 거동을 포착하는 데 모델에 미치는 영향은 어떠한가?
• RQ5 고해상도 영상 분할이 상변화 시뮬레이션의 강력하고 연속적인 검증을 가능하게 하는가?

주요 결과

• 시뮬레이션과 실험은 상계면의 형태와 진화에 대해 정성적으로 양호한 일치를 보였으며, 부력에 의해 유도된 대류로 인해 상단 부근에서 상계면이 가장 빠르게 진행됨을 확인하였다.
• 시뮬레이션과 실험 간의 상계면 절대 위치에 작은 그러나 측정 가능한 격차가 존재하며, 특히 600초와 900초 시점에서 두드러짐을 확인하였다.
• 얼음의 초깃emperature는 상계면 위치에 대해 낮은 민감도를 보였으며, -25 °C에서 0 °C 사이의 변동은 모델 차이를 최소화하였다.
• 이 격차는 시뮬레이션에서 비이deal한 경계 조건으로 인한 것으로 보이며, 예를 들어 완전한 단열 벽과 일정한 히터 온도를 가정한 것과 달리 실제 히터는 28 °C에서 33 °C 사이로 진동하기 때문이다.
• 플렉시글라스 벽을 통한 실험적 열손실이 열수확체로 작용하여 격차를 더욱 심화시킬 수 있다.
• 본 연구는 고해상도 영상 촬영, 고도의 분할 기법 및 시뮬레이션을 융합함으로써 다물리적 시스템에서의 통합 모델 검증의 가능성과 가치를 입증하였다.