[논문 리뷰] Fixed-grid sharp-interface numerical solutions to the three-phase spherical Stefan problem

Yavkreet Swami, Jacob Barajas외 1명|arXiv (Cornell University)|2026. 02. 09.

Phase Change Materials Research인용 수 0

한 줄 요약

논문은 고정 격자 예리 계면 방법을 유한 크기 입자에 대한 3상 spherical Stefan 문제로 확장하고, 밀도 및 운동에너지 점프를 포함하며, 2상 결과와의 검증을 통해 나노 규모 입자에서 운동에너지 효과가 중요한 반면 입자 크기가 커질수록 그 효과가 약해짐을 보인다.

ABSTRACT

Many metal manufacturing processes involve phase change phenomena, which include melting, boiling, and vaporization. These phenomena often occur concurrently. A prototypical 1D model for understanding the phase change phenomena is the Stefan problem. There is a large body of literature discussing the analytical solution to the two-phase Stefan problem that describes only the melting or boiling of phase change materials (PCMs) with one moving interface. Density-change effects that induce additional fluid flow during phase change are generally neglected in the literature to simplify the math of the Stefan problem. In our recent work [1], we provide analytical and numerical solutions to the three-phase Stefan problem with simultaneous occurrences of melting, solidification, boiling, and condensation in Cartesian coordinates. Our current work builds on our previous work to solve a more challenging problem: the three-phase Stefan problem in spherical coordinates for finite-sized particles. There are three moving interfaces in this system: the melt front, the boiling front, and the outer boundary which is in contact with the atmosphere. Although an analytical solution could not be found for this problem, we solved the governing equations using a fixed-grid sharp-interface method with second-order spatio-temporal accuracy. Using a small-time analytical solution, we predict a reasonably accurate estimate of temperature (in the three phases) and interface positions and velocities at the start of the simulation. Our numerical method is validated by reproducing the two-phase nanoparticle melting results of Font et al. [2]. Lastly, we solve the three-phase Stefan problems numerically to demonstrate the importance of kinetic energy terms during phase change of smaller (nano) particles. In contrast, these effects diminish for large particles (microns and larger).

연구 동기 및 목표

• 평면에서 구형 기하로의 3상 Stefan 문제 모델링을 유한 크기 입자로 확장한다.
• 모든 계면에서 밀도 및 운동에너지 점프 조건을 유지한다.
• 구형 좌표에서 고정 격자 예리 계면 immersed-boundary 방법을 개발하고 적용한다.
• 이동 계면을 정확히 시작하기 위해 작은 시간 해석해를 사용하여 초기화한다.
• 알려진 2상 결과에 대해 접근법을 검증하고 입자 크기에 따른 운동에너지 항의 역할을 분석한다.

제안 방법

• 밀도, 에너지, 운동을 포함한 계면 점프 및 Dirac 델타 항을 포함하는 고체, 액체, 증기에 대한 지배 방정식을 형식화한다.
• 지배 방정식을 무차원화하고 초기화를 위한 작은 시간 해를 도출한다.
• 구형 좌표에서 이동 계면의 immersed-boundary 추적을 갖는 2차 고정 격자 예리 계면 방법을 구현한다.
• 저밀도-비율 및 고밀도-비율 한계에서의 작은 시간 해석을 사용하여 초기 계면 위치/속도와 온도장을 얻는다.
• 유한 크기 구에서 3상 구형 Stefan 문제를 해결하여 계면 진화 및 밀도와 운동에너지 항의 영향을 연구한다.
• Font 등 의 2상 나노입자 용융 결과를 재현하여 수치적으로 검증한다.

(a) Two-phase spherical Stefan problem

실험 결과

연구 질문

• RQ1 유한 크기 입자에서 동시에 용융 및 끓음 전면을 가지는 3상 구형 Stefan 문제는 어떻게 진화하는가?
• RQ2 다양한 입자 규모에서 밀도와 운동에너지 점프 항이 계면 위치, 속도, 용융/끓음 시간에 미치는 영향은 무엇인가?
• RQ3 다상 Stefan 문제에 대해 구형 좌표의 고정 격자 예리 계면 방법의 정확도는 얼마나 되는가?
• RQ4 나노 규모 입자에 비해 마이크로 규모 입자에서 운동에너지 항이 용융 시간을 얼마나 영향을 미치는가?

주요 결과

• 운동에너지 항은 더 작은(나노) 입자에서 용융 시간에 상당한 영향을 준다.
• 더 큰 입자(마이크로 규모 이상)에서는 운동에너지 항의 영향이 줄어든다.
• 고정 격자 예리 계면 방법은 공간 및 시간에서 2차 정확도를 달성한다.
• Font 등 의 2상 나노입자 용융 결과를 재현하여 방법을 검증한다.
• 3상 시뮬레이션은 동시 용융 및 끓음 동안 각 상에서의 밀도 및 운동에너지 항의 역할을 정량화하며, 작은 입자에서의 중요성을 강조한다.

(b) Three-phase spherical Stefan problem